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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. 
They should be accompanied by asummary at the beginning of the paper conveying briefly 
the content of the paper, and drawing attention to important new information and to the 
main conclusions. The summary should be intelligible in itself, without reference to 
the paper, to a reader with some knowledge of the subject; it should not normally exceed 
200 words in length. Authors are requested to submit MSS. in duplicate. These should 
be typed, using double spacing and leaving a margin at one side. Corrections to the MSS. 
should be made in the text and not in the margin. Bw Council decision, MSS. of accepted 
papers are retained by the Society for one year after publication; unless their return 
is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested 
to follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or 
in a numbered list at the end of the paper. Each reference should give the name of the 
author cited, irrespectively of the occurrence of the name in the text (some latitude being 
permissible, however, in the case of an author referring to his own work). 


3. Notation.—Authors should conform closely with the recommendations of Com- 
mission 3 of the International Astronomical Union (Trans. I.A.U., Vol. VI, p. 345, 
1938). Council has decided to adopt the I.A.U. 4-letter abbreviations for constellations 
where contraction is desirable (Vol. IV, p. 221, 1932). 


4. Diagrams.—These should be drawn about twice the size required in print and 
prepared for direct photographic reproduction except for the lettering which should be 
inserted in pencil. Legends should be given in the manuscript indicating where in the 
text the figure should appear. Blocks are retained by the Society for 10 years; unless the 
author requires them before the end of this period they are then destroyed. 


5. Tables.—These should be arranged so that they can be printed upright on the 
page. 

6. Proofs.—Costs of alteration exceeding 5 per cent of composition must be borne 
by the author. Fellows are warned that such costs have risen sharply in recent years, 
and it is in their own and the Society’s interests to seek the maximum conciseness and 
simplification of symbols and equations consistent with clarity. 


7. Revised Manuscripts.—When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MS. 


Reading of Papers at Meetings 


8. Ifa paper is to be read at a meeting it should be received by the Assistant Secretary 
not later than the previous Friday. 


g. When submitting papers authors are requested to indicate whether they will be 
willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


10. Postcards giving the programme of each meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether for the Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1950 JANUARY 13 
Professor W. M. Smart, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 
Lt.-Col. William Arthur Garstin, C.B.E., 65 Campden Hill Court, London, 
W.8 (proposed by E. A. Milne); 
*Franz Daniel Kahn, 41 St. Margaret’s Road, Edgware, Middlesex (proposed 
by S. Chapman) ; 
Thomas David Kinman, University Observatory, Oxford (proposed by 


H. H. Plaskett) ; 

Joan Lilly, Waveney High School, Wortwell, Harleston, Norfolk (proposed 
by J. G. Read); 

Michael Selwyn Longuet-Higgins, Trinity College, Cambridge (proposed by 
H. Jeffreys) ; 

Hugh Ryan, B.Sc., 127 Princess Court, London, W.2 (proposed by F. M. 
Holborn) ; 

William Court Parsons Tapper, O.B.E., M.I.E.E., Godrevy, Whitecliff Road, 
Parkstone, Poole, Dorset (proposed by A. Hunter); 

Hugh Harold Lawrence Carmichael Vredenburgh, 10 Beauchamp Place, 
London, S.W.3 (proposed by W. H. Steavenson) ; 

Gert Claes Johan Walén, Roslagsgatan 25, Stockholm, Sweden (proposed 
by S. Chapman); and 

Kevin Charles Westfold, Queen’s College, Oxford (proposed by S. Chapman). 


The election by the Council of the following Junior Members was duly 
confirmed :— 


John Andrew Bastin, 24 The Drive, London, E.17 (proposed by H. H. Plaskett) ; 

Lewis Grenville Bentley, 13 Romeyn Road, Streatham, London, S.W.16 
(proposed by F. M. Holborn); 

Daphne May Fooks, H.M. Nautical Almanac Office, Herstmonceux Castle, 
Sussex (proposed by F. M. McBain); 


* Transferred from Junior Membership. 











Meeting of 1950 January 13 Vol. 110 
Mavis McIntyre Smith Gibson, B.Sc., H.M. Nautical Almanac Office, 
Herstmonceux Castle, Sussex (proposed by F. M. McBain); and 


Angela Margaret James, B.Sc., H.M. Nautical Almanac Office, Herstmonceux 
Castle, Sussex (proposed by F. M. McBain). . 


One hundred and eighty-one presents were announced as having been received 
since the last meeting. 


The President announced that the Council had awarded the Gold Medal of 
the Society to Professor Joel Stebbins, late Director of the Washburn Observatory, 
Madison, Wisconsin, U.S.A., for his development of physical methods in astro- 
nomical photometry and for the results obtained by the use of those methods. 





ON A PERIODIC FLUCTUATION IN THE 
LENGTH OF THE DAY 


H. F. Finch 


(Communicated by the Astronomer Royal) 
(Received 1950 January 10) 


Summary 


From a study of the performance of a number of quartz-crystal clocks 
employed in the Greenwich Time Service an annual periodic fluctuation in the 
length of the day is deduced. The variation is of the order of -+-0%-oo01 and has 
an accumulative effect in time of approximately +0%-o60. This is in very 
close agreement with results obtained from independent data by N. Stoyko 
and demonstrates the persistent character of the phenomenon. 





1. With the development of the quartz-crystal oscillator as a frequency 
standard, a demand has arisen for time signals capable of providing frequency 
checks to a precision consistent with the standard of performance of the crystals 
themselves. In undertaking to furnish these signals the Royal Greenwich 
Observatory has at its disposal a number of quartz clocks which serve to bridge 
the gap between the time determinations, the observational errors of which they 
also serve to smooth out. But to control the emission of a system of daily time 
signals that will define intervals of 24 hours to an accuracy as good as or better 
than one millisecond has proved a very intractable task. Of the difficulties 
encountered, only that of establishing the length of the current day, or alternatively, 
of determining the rates of the observatory clocks with the desired accuracy, is 
discussed in this paper. 

2. Beat-frequency measurements have shown that over short periods of time 
frequency stabilities of the order of one part in 10" can be achieved with quartz 
crystals, while the frequency difference of two crystals whose characteristic 
performances are known, can be predicted, in certain circumstances, over a period 
of months to within several parts in 108. In justice to the crystal, therefore, if the 
frequency of a standard oscillator is to be quoted in cycles per second, the length 
of the second must be known to something better than one part in 108. In other 
words, the astronomical observations must be capable of establishing the current 
rates of the observatory clocks to within 0*-oo1/day. 

3. At the Greenwich Observatory the available accuracy tends to be limited 
by two factors: (a) The time determinations have, up to the present, been made 
visually and it has not proved possible to reduce the probable error of a single 
determination much below +0*-020. If the clock rate, assumed to be constant, 
is determined from two such observations, they would need to be separated by 
an interval of months for the rate to be derived with certainty to the desired 
accuracy. The interval can be reduced by increasing the number of time 
determinations, but the gain is by no means proportional to the increased effort. 
(b) The rates of crystal clocks are not constant, though they may be very uniform. 
A quartz crystal, when first set into oscillation, shows a marked frequency (rate) 
drift. The drift gradually diminishes until, after some months, it settles down to 

r* 
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a fairly constant value. The clock correction (Z) can then be represented very 
closely by the expression E=a +bt + ct®, where a, 6 and c are constants and ¢ is 
reckoned in time from some arbitrary epoch. When t is expressed in days and E 
in seconds, the value of c is found in general to be between the values of + 0-00005. 
The daily rate (R) is given by R=b+2ct. Reliable prediction of the error of 
the clock requires c to be determined with great accuracy, because its effect increases 
with the square of the time. The existence of a frequency drift greatly increases 
the interval required to establish the current rate of the clock. Moreover, even 
the best clocks are liable to small erratic changes of 6 and c, and although 
sudden changes of rate are shown up in the daily intercomparisons, minor 
variations in the value of the frequency drift are not so readily detected. In 
order to obtain increased accuracy, and also to allow for the contingency of 
failures, a considerable number of clocks must be available for use in the service. 

4. After quartz clocks had been introduced into the Time Department and 
experience with them had established their high standard of performance, the 
policy of predicting their errors on a long-term basis was introduced. The clock 
constants a, b and c were derived from a series of time determinations extending 
over many months. The extrapolated corrections were then heavily weighted 
against those indicated by current observations, which might appear discordant 
because of large observational errors. As observations accumulated, any 
persistent departure of the clock corrections from the ephemeris was allowed 
for graphically and the constants were subsequently adjusted to secure the best 
fit over the last few months. Using the adopted current corrections to the 
standard clocks and the daily intercomparisons it was possible to deduce the 
amount by which the clock controlling the time signal should be corrected prior 
to transmission. 

5. The method appeared for a while to work fairly satisfactorily, but in the 
autumn of 1945 the observed clock corrections began rapidly to diverge from 
the predicted. The magnitude of this divergence varied somewhat from clock 
to clock, but all appeared to have undergone an acceleration of rate. A change 
in the period of rotation of the Earth was suspected. It was then noted that there 
had been a similar occurrence, though to a much less marked degree, during the 
previous autumn, which suggested the existence of a seasonal effect. There 
appeared to be little evidence, however, of a compensating retardation of the clocks 
during the spring. In subsequent years, as more clocks became available, the 
apparent autumnal acceleration was again observed. When, therefore, in the 
autumn of 1948 a discordance between the observed and predicted clock 
corrections again began to show, the weight given to the current observations 
was greatly increased. In the spring of 1949 there appeared a fairly definite 
indication of a retardation amongst a number of the clocks. The agreement 
was qualitative rather than quantitative, but it was felt that the differences between 
individual clocks might easily arise from erratics and the difficulty of establishing 
the true frequency drifts. 

6. The strong suspicion of a seasonal fluctuation aroused by the spring 
changes led to an examination being made of the final corrections to the times of 
emission of the Rugby time signals, which are based on the subsequent assess- 
ment of the performance of the clocks. Since the time signals are under the 
control of the Greenwich Time Service, the corrections are a measure of the errors 
in the predicted clock corrections irrespective of which clocks are in use as 
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standards at the time of transmission. The four years 1944-47 agree in showing 
a pronounced seasonal fluctuation in the signal trend, i.e. the mean daily increase 
in the value of the signal corrections (signal slow). It does not appear before 
1944 because the technique of long-term prediction had not then been fully 
developed; it disappears in 1948 because of the autumn change in technique. 
The mean seasonal trend of the signal for the four years 1944-47 is shown in 
Fig. 1. The form of this curve has necessarily been influenced to some extent 
by the discretion of the predictor in modifying his adopted clock corrections in 
the light of accumulating observations. 


ms. day Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. 


Nov. Dec. 


‘ ‘ 


2:0 - 
0-0 
26° 


Fic. 1.—Mean “ trend’’ of GBR time signal at 10.00 U.T. for the years 1944-1947. 


7. The true form of the fluctuation must be derived from the observed 
corrections. ‘The material suitable for this purpose is limited by various factors. 
It is possible to use only those clocks which have been free from disturbances 
occasioned by faults and have had undisturbed runs for many months. For the 
period 1943-47 the usable data were provided by four clocks, 44, 4B and 4C in 
service at the Post Office, and Ar at the Royal Observatory, Edinburgh (where 
a time service was operated during the war). Considerable trouble was 
experienced during this period with the clocks at Abinger. From the beginning 
of 1948, following the replacement of the original defective beat counters of the 
Abinger clocks, the installation at Greenwich of six crystal clocks, and further 
additions at Dollis Hill, a greatly increased amount of material became available. 
From about the middle of 1947 the observational data were strengthened by a 
considerable increase in the number of time determinations and by an improve- 
ment in the accuracy consequent upon the installation at Abinger of a better 
transit instrument. The material falls naturally into two groups ; separate 
analyses have therefore been made for the two periods 1943-47 and 1948-49. 

8. The analysis has been based on the assumption that the observed clock 
corrections can be fitted to an ephemeris of the form 


E=a+bt+ct* +f, 


where f is a periodic function of the time representing a hypothetical seasonal 
fluctuation, having therefore a yearly period and arising from some cause 
independent of the clocks. The second differences of a series of clock corrections 
determined at equal intervals are then of the form 2c +A"f and will differ from 
clock to clock only through the different values of the frequency drift. The 
third differences, A” f, correspond with those that would be obtained from a 
clock having zero values for a,b and¢. Random «: il! be present, but their 
effect is reduced by taking means of the various values obtained for A”f throughout 
set intervals of the calendar year. Nevertheless, owing to the possible small 
variations in rate and frequency drift and to residual random errors, it is unlikely 
that the sum of the mean values of A” f for the whole twelve months will be zero. 
The form of their distribution is not affected, however, by reducing each term 
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by its mean value for the year. The adjusted values can be integrated to give 
a series of mean values for A’f, the constant of integration being adjusted so that 
the mean yearly value is again zero. ‘The integration can be continued to give 
first the values of A’f and then those of f itself. 

g. The formation of the third differences was greatly eased by employing 
the residual clock corrections, which had been derived for the purpose of 
prediction, with reference to a parabolic ephemeris. These residuals had already 
been corrected for the effects of sudden discontinuities due to stoppages of the 
phonic motors and their subsequent restarting; such stoppages do not affect 
the rates of the clocks. The mean residual clock corrections with the corre- 
sponding mean dates of observation for each month were tabulated, after which 
the mean residuals were corrected for the effect of polar variation. Up to the 
end of 1948 these corrections were based on the values of the polar motion 
published by the Central Bureau of the International Latitude Service. ‘The 
1949 corrections were based upon the latitude variation data supplied by the 
U.S. Naval Observatory at Washington. The data were then adjusted to an 
adopted date for each month by correcting for the residual rate. As the correction 
was always small, the effect of departures of the mean from the standard date was 
never critical. The standard dates were selected to be separated as nearly as 
possible by alternate intervals of 30 and 31 days in order to make the data suitable 
for differencing. The weights of the final monthly values necessarily differ 
considerably. 

10. In order to reduce the effect of observational error in relation to the 
fluctuation, the months were divided into two groups, odd and even, to give a 
uniform working interval of 61 days. The two sets of figures were separately 
differenced and independently reintegrated as described above. Two sets of 
values for f were thus obtained, which were finally dovetailed together. This 
treatment of the figures required a minimum undisturbed run of seven months 
for any clock to furnish a single third difference. 

11. The values of the third differences for the two periods 1943-47 and 
1948-49 are given in Table I. The clock and the year are indicated in the first 
two columns and at the head of each succeeding column is listed the calendar 
month corresponding to the centre of each seven-month period involved in the 
derivation of the figures below. At the foot of the table are given the mean 
values for each column and the weight as defined by the number of entries above. 
Table II shows the number of time determinations obtained in each month 
during the whole period 1943-49. 

12. The reintegration of A” f through the stages A’f, A’f to f is contained 
in Table III; the left- and right-hand sides are devoted respectively to the odd 
and even months. The mean value printed at the bottom of each column is 
subtracted from each of the values given above it to obtain A’f, A’f and f as 
explained above. Fig. 2 contains plots of A”f, 1/61.A’f and f for each of the 
two periods. It will be noted how closely the results from the odd and even 
months separately are in agreement. The analysis of the f curves give: 


for 1943-47. f= +55°7sinn(T +141)° — 16-8 sin 2n(T +67)’, 
and for 1948-49 f= +58-3sinn(T+158)°— 7-1sin2n(T +48)°, 


where T is the day of the year, f is expressed in milliseconds and n is the diurnal 
motion of the mean Sun in degrees of arc. i.e. 0-986. 
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The expressions represent that part of the clock correction caused by the 
seasonal variation, which is now clearly apparent. The effect upon the observed 
clock rates is obtained by differentiating these expressions with respect to T, 
when we obtain: 

for 1943-47 f’ = —0-958 sin n(T + 49)° +0°578 sin 2n(T +22)°, 
and for 1948-49 f' = —1-003 sin n(T +67)° +0-245 sin 2n(T +2)°, 
the unit here being ms./day. 

13. The similarity of the curves for the two periods is very marked, 
the deep trough occurring in August being the most striking feature of the A” f 
curves. The second trough appearing in the A”f curve for 1943-47 is the 
enhanced effect of the second harmonic which is very apparent in f during this 
period, but much less so during the period 1948-49. It is the enhancement 
of this harmonic which causes the f’ curves to differ more markedly than the 
f curves. In both periods, however, the clocks appear to go fast in the spring 
and slow in the autumn, the minima and maxima of the f, or correction curves, 
occurring in May and October respectively. The date on which f passes through 
zero on its upward trend occurs early in August during both periods, at the time 
when the apparent rate, as seen in f’, isa maximum. The amplitudes of the two 
f curves are in close accordance and are about + 08-060, while those of the f’ curves 
indicate an annual fluctuation in the apparent rates of + 08-oo1/day. 


TABLE II 
Number of Observations Secured 
Jan. Feb. Mar. Apr. May June July Aug. .  “Oer, 


se ; ‘ a Me so 27 33 26 
20 24 31 25 29 19 26 25 35 26 27 
23 17 29 24 20 15 I4 23 20 24 
12 19 25 22 17 15 18 16 30 F 16 
16 10 10 15 16 10 14 24 24 21 
28 7 7 25 21 32 37 62 < 51 22 
29 48 65 5° 38 22 28 25 47 a 25 
32 48 47 43 63 47 31 36 ve “° : 


14. Having determined the seasonal fluctuation, the observed clock 
corrections can be corrected for its effects, and a new set of values for a, 6 and c 
for the various clocks determined. This has been done for the Post Office 
ring-crystal clock 5B, which from the daily intercomparisons was known to be 
amongst those whose performance was the most uniform. The full line in 
Fig. 3 represents the performance of the clock referred to the finally adopted 
G.M.T., after removal of a parabolic ephemeris and the polar effect; the dotted 
line represents its performance further corrected by the f term. The deviations 
from the ephemeris over the eighteen months ending 1949 June 15, when a 
failure of temperature control occurred, do not exceed +0%.020. Some of the 
irregularities may be due to observational uncertainties, though at least one 
is apparent in the intercomparisons with other clocks. Nevertheless, the 
performance is truly remarkable and illustrates the difficulties with which the 
astronomer is faced when asked to check the frequencies of crystal oscillators 
by using the Earth as a standard. 

15. The seasonal fluctuation f is one which clearly cannot be ascribed to 
the clocks, nor can it readily be associated with systematic errors in the observa- 
tions due to faulty star places, instrumental defects or to local peculiarities. 
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1943-1947 1948-1949 


Fic. 2.—Curves of A” f, f' and f. 


It is present irrespective of which clocks are employed, and persists in the face 
of changes of observers, instruments and observing sites. If there were any 
doubt about this it would be dispelled by the fact that the fluctuations are 
revealed also by the Washington time observations. Reception at Abinger of 
the WWV radio signals renders it possible to obtain a close comparison between 
U.T. as determined at the Naval Observatory, Washington, and the Royal 
Observatory. After removal of a systematic difference of +0%-016* and the 


* Travel time of the signal, if allowed for, would make this about +0*-036. 
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effect of polar variation, the monthly mean differences over the four years are, 
in milliseconds, as follows :— 


Jan. +14 May -—I0 Sept. —10 
Feb. +14 June —13 Oct. —- I 
Mar. +II July -11 Nov. + 9 
Apr. + 2 Aug. —I2 Dec. + 8 


These are given as corrections to be applied to the ‘‘Greenwich”’ f curve to 
obtain the corresponding Washington curve. There appears to be a slight 
annual variation between the two observatories, but this is small compared with 
the main fluctuation. Fig. 2(d) contains a plot of the mean f curve for the two 
periods referred, A to the American observations and B to the British observations. 


The difference shows itself as a small phase shift and a slightly increased amplitude 
in the American curve. 


Fic 4.—Curves showing the effect of an error in the adopted value of c. 


16, It may seem surprising that the sinusoidal character of the phenomenon 
had not been earlier recognized. During the period that the group 4 clocks were 
in operation, 4C was for good reasons acknowledged as being the most reliable. 
The value of c for this clock had been determined and the corresponding values 
for the others were readily derived by using their frequency drifts relative to 4C. 
It would now appear that the adopted value ( +0-00005) for the frequency drift 
of 4C was slightly in error, though it provided a good fit over a considerable 


portion of the f curve. The resulting qualitative effect can be seen on reference 
to Fig. 4 in which the curve f represents the seasonal fluctuation and e that part 
of the parabolic ephemeris resulting from the error in the adopted value of c. 
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Beneath is given the residual curve (f-e). Over the interval X,-X, the 
deviations from linearity are too small to be revealed by observations subject 
to a comparatively large scatter in the presence of minor variations in the clock 
rate. At X, the curve rapidly dips. Since no relative change has occurred 
in the various values of c all the clock curves will behave in the same way and the 
natural interpretation would be that a sudden change in the rotation period of 
the Earth had occurred. Following the initial run-off, the original value of c will 
again appear to hold, but at the end of twelve months the Earth will once more 
appear to slow up. 

17. The claim to have established for the first time an annual fluctuation 
in the Earth’s rotation was made by Pavel and Uhink *, based on an investigation 
of the performance of two quartz-crystal clocks at the Geodetic Institute, Potsdam, 
during the period 1934 May—1935 June. Referred to the Potsdam time deter- 
minations, both clocks showed similar fluctuations in the residuals from a parabolic 
ephemeris. As the time determinations did not show similar fluctuations 
relative to the mean of the times determined at nine observatories, it was 
concluded that the fluctuations were not due to errors of observation but to a 
variable rate of rotation of the Earth. ‘The material employed did not justify 
this conclusion, neither clock had an unbroken run during the period discussed,. 
and each was used to smooth over the disturbance of the other when a stoppage 
occurred. ‘The deduced variations in the rotation, with a range from the mean 
of from +0-005 to —0-004 secs./day, are too large to be real. 

On the other hand, N. Stoyko, by comparing the departures of mean monthly 
rates from a linear rate for the clocks at Paris, Washington and Berlin during the 
three years 1934-37 (which cover the period of investigations by Pavel and Uhink), 
obtained a result + in close agreement with that found in the present paper. 
This paper seems to have attracted little attention, probably because further 
confirmation of the reality of an annual fluctuation in the Earth’s rotation appeared 
to be required. 

The present investigation had been completed when Stoyko published, 
without details, the results of a study of the rates of the clocks at Paris and 
Washington from 1946 January to 1947 April {, which confirmed his previous 
conclusions. This period is covered by the Greenwich observations; in Fig. 2 (e) 
Stoyko’s results for the variations in the rate of rotation are plotted in dots 
upon the Greenwich f’ curve for the period 1943-47. Stoyko’s representation 
of his results can be expressed in the form: 

f’=—1-18sinn(T +59)° +0-43 sin 2n(T + 42)°, 
which is in very close agreement with the results obtained by the Greenwich 
clocks. Stoyko’s conclusions about the existence of an annual fluctuation in the 
period of rotation of the Earth are completely confirmed by the present investi- 
gation; moreover, the similarity in amplitude and phase found during the 
different periods show the effect to be a persistent one. 

18. In a recent paper§, van den Dungen, Cox and van Miegham, using 
observed barometric pressure over the surface of the Earth at opposite seasons, 
have shown that the effect can be accounted for, at least in part, by changes in 


* F. Pavel and W. Uhink, Ast. Nach., 257, 365, 1935. 

+ N. Stoyko, Comptes Rendus, 203, 39, 1936; and 205, 79, 1937- 

t N. Stoyko, Bulletin de la Classe de Science, Académie royale de Belgique, 5e Serie, XX XV, 1949. 

§ F. H. van den Dungen, J. F. Cox and J. van Miegham, Bulletin de la Classe de Science, Académie 
royale de Belgique, 5e Serie, XXXV, 1949. 
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the Earth’s moment of inertia caused by seasonal changes in the distribution of 
air masses. ‘The data are incomplete and further investigation is required, but 
the paper constitutes a first attempt to account quantitatively for the phenomena. 

19. Since the fluctuation occurs in the unit of time it will introduce a corre- 
sponding fluctuation into the observed longitude of the Moon. The observations 
of occultations have shown for many years a fluctuation of approximately annual 
period in the Moon’s mean longitude. Brouwer and Watts have analysed the 
results of the occultation observations and of the Washington meridian 
observations of the Moon*, deriving corrections to the elements of the orbits 
of the Moon and the Earth. They suggested that one term, not otherwise 
accounted for, might be attributed to an annual variation in the rate of rotation 
of the Earth. The variations shown by Stoyko’s and by the present investigation 
are, however, much too small to account for this term in the Moon’s longitude 
for which, therefore, some other explanation must be sought. 

20. In conclusion, the practical importance of the fluctuation and its effect 
upon the estimated frequency of standard oscillators must be emphasized, 
together with the need for a precise definition of the Mean Solar Second which 
takes into account both this and the effect of the polar variation upon time. 


The author wishes to acknowledge the assistance received from his colleagues 
in the present investigation and the close cooperation of the Post Office 
Engineers, without which the data available would have been greatly restricted. 


Royal Greenwich Observatory, 
Abinger : 
1949 December 15. 


* D. Brouwer and C. B. Watts, Astron. ., 52, 169, 1947. 


ADDENDUM 
Since going to press a further paper by W. Uhink (Ast. Nach., 278, No. 3, 
1950) dealing with the above subject has now become available. This paper 
contains an account of results obtained during recent years in Germany. 
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Summary 


Relative gradients for 166 southern stars were determined by photographic 
spectrophotometry with a slitless spectrograph attached to the .30-inch 
Reynolds reflector. The system is believed to be uniform with that of 
Greenwich. The relations between the Mount Stromlo gradients and the 
colours of other observers are examined, and the colours of a group of very 
blue B stars discussed with regard to their bearing on the spectrum-colour 
relation and on the stellar temperature scale. 





1. Introduction.—This paper describes observations of the relative spectro- 
photometric gradients of anumber of the brighter southern stars. The programme 
was designed as a southward extension of the Greenwich relative gradient pro- 
gramme (1, 2), and follows Greenwich closely, both in instrumental equipment 
and in general observing procedure. Observations were commenced in 
1948 May with a slitless spectrograph attached to the 30-inch Reynolds reflector. 
By 1949 June about 2500 spectra of 170 odd stars had been taken, yielding 
gradients with standard errors of 0-02 or 0-03, a satisfactorily small figure and a 
comparable one to that attained at Greenwich. At this stage it was decided to 
suspend the programme, with the intention of resuming it with lower dispersion 
and so reaching fainter stars. The present is therefore a convenient time for 
describing the results so far obtained. This is especially so, as the programme 
includes a number of interesting stars for which no other spectrophotometric 
data exist at all. 

The first part of the paper describes the observing procedure and concludes 
with the catalogue of gradients. In the second the Mount Stromlo gradients are 
compared in detail with those of a number of other observers; an apparent 
anomaly in the colours of « C Maj and « Cari is examined; and the colours of a 
group of very blue stars are discussed with regard to their bearing on the spectrum- 
colour relation and on the stellar temperature scale. 

2. The Spectrograph.—The spectra were taken on the 30-inch Reynolds 
reflector, on a slitless spectrograph similar to that described by Davidson (3), and 
built in the Observatory shops to a design by Dr Allen and Mr Gottlieb. The 
optical arrangement is shown in Fig. 1. The dispersion at Hy was 45 A./mm. 
The diagonal mirror D immediately before the plate was mounted on a spindle 
about which it could be given a small to-and-fro rotation by the usual heart-shaped 
cam-and-lever arrangement, to widen the spectra to }mm. Irregularities in the 
telescope drive presented difficulties for a while, by introducing streaks into the 
spectra. A careful remachining of the main worm, carried out at the Mary- 
byrnong Defence Research Laboratories, substantially reduced the trouble, 
which was finally almost eliminated by the adoption of a simple guiding device. 


* Received in original form 1949 November 28. 
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The third (plane) mirror A was mounted so that it could be rotated about an axis 
coincident with the axis of the primary mirror, the amount of rotation being 
controlled by an arm and screw at B. ‘The star under observation was centred 
on the guiding eyepiece cross-wire, and during exposure was maintained there, 
an appropriate rotation of the screw B counteracting any tendency to wander, 
induced by faulty driving. ‘The scheme worked well, and only a few spectra had 
finally to be rejected because of streaks. 


\ \\ 
\\ \\ 


. me. 
Diagonal D0 —____.,\ 27 
\ . . 
Guiding Telescope aoa: Fa a. 









































7) 
maa 
(f) 


Fic. 1.—(a) The optical train. (b)-(f{): Diaphragms inserted into the collimated beam at CC. 
B is a screw and arm which rotates A about an axis coincident with the axis of the primary mirror. 
The diagonal D is rocked to widen the image. 


3. Plate Calibration.—Plates were calibrated at the telescope, by taking a 
number of equal-interval exposures of the same star, with a series of diaphragms 
inserted in turn into the collimated beam (in the plane CC) to give varying reduc- 
tion factors. Two sets of diaphragms were constructed. The first, which 
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consisted of six half-magnitude steps, took the form of Maltese crosses 
for the 4 and 1 magnitude steps, and of plates with four equal symmetrical 
holes for the 14, 2, 24 and 3 magnitude steps. (See Fig. 1(6) and (c).) The 
calculation of the fraction of incident light transmitted by a diaphragm involved 
a calculation of the cross-sectional area of the beam, and this was made difficult 
by the complicated nature of the obstructions introduced by the secondary and 
diagonal mirrors and their supports. However, Mr Hogg and Miss Mason, who 
were working on the same spectrograph on a programme of monochromatic 
photoelectric photometry, kindly undertook a direct calibration of the diaphragms 
with their photocell. They repeated this twice, with results which agreed to about 
0-O1 mag. in the mean, and were of ample accuracy for photographic work. 

The constancy of the transmission factors of diaphragms of this sort depends 
on the constancy of cross-section of the collimated beam, and would have been 
upset by any such variations as, say, those introduced by changes in the setting 
of the Cassegrain paraboloid. Variations of this sort were apt to occur on the 
replacing and realigning of the spectrograph after it had been removed to make 
way for other programmes, and did on one occasion lead to a change (of 0-03 mag.) 
in the constants of diaphragms (Fig. 1). However, it was found possible to check 
the beam cross-sections, and indeed to get a good idea of the area of most of the 
obstructions, by setting on a bright star and exposing a half-plate at CC. With 
Canopus (— 0-8 mag.) an image of the beam cross-section was recorded in about 
a minute. 

Moreover, as reference to Fig. 1 will show, the transmission factors for the 
Maltese crosses (6), and especially the relative transmissions for the plates (c), 
should have been independent of the beam cross-section. Our limited evidence 
on this point was certainly not opposed to this conclusion. (It was also found that 
the relative transmissions for plates (c) as measured by the photocell agreed to 
within 0-01 mag. with those computed from the hole diameters.) Whenever 
possible the step I-1} mag. was avoided in deriving the actual calibration curves, 
and, this precaution having been adopted, the writer feels that errors arising from 
possible variations in the beam cross-section have been effectively eliminated. 

Although no trouble was expected from this source, it seemed desirable to 
carry out a completely independent basic calibration, and a second set of dia- 
phragms was accordingly constructed. These are illustrated in Fig. 1 (d, e and f) 
and worked on the half-aperture principle. In calibrating a plate, diaphragm (d) 
was first inserted, and three successive exposures made, the first with (d) only, 
and the remaining two with (e) added in such a way as to cut out the east and west 
halves of the mirror inturn. ‘The order of these exposures was varied from star 
to star. The three diaphragms (f) were used to reduce the light of brighter stars. 
These three proceeded at approximately }? magnitude steps and gave, with (e), a 
total range of three magnitudes. ‘This method of calibration was independent 
of the beam cross-sections. 

To verify the equivalence of these two methods of calibration there were 
available 42 programme stars for which, with a few exceptions, at least two 
determinations of gradient had been niade by each method. For the 42 stars 
the unweighted mean difference (diaphragm gradient)-(half-aperture gradient) 
was —0-01I1 +0-008(s.e.). For the 16 bluest stars of the group it was 
—0-015 +0-012, and for the 10 reddest stars it was +0°005+0-013. These 
figures are considered satisfactorily small. Not much significance is attached to 

2 











18 S. C. B. Gascoigne Vol. 110 


the slight systematic trend from blue to red, especially as one would have expected 
the effects of faulty calibration to be most marked for the red stars, for which, on 
our plates, the red end was much more heavily exposed than the blue. 

4. Plate Measurement.—We had intended to use Eastman 103-F plates but 
as they could not be obtained in time Kodak Panatomic-X were used instead. 
Apart from the drawback that their red sensitivity did not extend beyond about 
A6400, these Pan-X plates were very satisfactory, with fine grain, low contrast 
(an advantage in this work) and speed enough to record a fifth magnitude early- 
type star in three minutes. 22 spectra were normally secured on one quarter- 
plate, usually, although by no means invariably, all with the same exposure time, 
which was 20, 60 or 180 seconds. Each star was taken twice, if possible with 
different diaphragms, so that its spectra could be obtained at different density 
levels. If the calibration was performed with the first set of diaphragms four 
exposures of the calibrating star were usually made, to give four calibrating spectra 
(sometimes three or five). Withthe second set of diaphragms only three exposures 
were necessary, and two stars were then used to calibrate a plate. 

Plates were developed four at a time, in a tank with a plunger for agitating the 
developer. Kodak D-19 was used throughout. In the wave-length range used, 
46209 to A4050, three calibration curves usually represented well enough the 
variation of contrast with wave-length, and these curves were assumed common 
to all four plates of any one batch. 

We measured the plates in the projection photometer described previously 
(4), at the following eleven wave-lengths: AA 6209, 6035, 5825, 5550, 5259, 4945, 
4660, 4510, 4250, 4170, 4050. The harmonic mean is 44938. ‘These wave- 
lengths were chosen with some care, to ensure that they were as free as possible from 
absorption lines, and so that they would avoid regions where the plate sensitivity 
was changing rapidly. Stars of spectral type later than about F6 were measured 
only at the first eight wave-lengths, as over the last three they showed the deficiency 
of radiation noted by Plaskett (§) and others. For the same reason it was decided 
not to measure farther into the violet than A 4050, as beyond that point early A stars 
on our plates are weakened, presumably because of confluence of the wings of the 
hydrogen lines. 

The usual assumption was made, that the relative monochromatic magnitude 
difference Am, between two stars varies with wave-length according to the law 
a+b/A. The relative gradient is then 0-921 5 (A in microns), and was found from 
a least squares solution for a and 6, the Am, having been measured as above. 

One small but important addition to the photometer was that of a rheostat 
in the plate circuit for varying the amplification. This provided a means for 
holding the clear glass deflection at 100, and so enabled transmissions to be read 
directly as percentages, a considerable saving in time and clerical work. Plates 
were measured once only, previous experience, and a number of trials with the 
present series, having shown that there is little to be gained from a second 
measurement. 

It has been pointed out that spectra were not always taken with the same 
exposure times, and we had to find whether this necessitated a correction. For 
this purpose we had available about thirty pairs of spectra, taken with exposure 
times either of 20 and 60 sec., or of 60 and 180 sec., and hence always in the ratio 
of three to one. We measured and reduced these on the lines suggested by 
Greaves (6) but in spite of the appreciable reciprocity failure shown by 
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Panatomic-X (corresponding to a Schwarzschild exponent of about 0-8), it was 
found that with our small exposure time ratio the correction could be neglected, 
even for the reddest stars on the programme. 

5. The Standard System.—We selected as standards eighteen stars of early 
type, well distributed about the sky, which were intercompared frequently, 
and relative to which the gradients of all other stars were determined. The 
selection is set outin Table I. It will be noted that seven of the eighteen are also 
Greenwich standards. All comparisons of standard stars were made at equal 
altitude. It was often possible to compare a run of three, four or sometimes 
more stars which would pass through a common altitude at closely succeeding 
times. Such groups were picked from a most convenient graph on which sec z 
for each standard star was plotted against sidereal time. Fig. 2 summarizes the 
scheme, each line representing a series of comparisons between the two stars it 
joins. A number of comparisons of low weight were omitted from this diagram. 


Tasre I 
The Standard Stars 
Star H.D. type Ss Wt. G 

Y Pegs . B2 —25 “4 
« Erid . Bs —23°6 
B Aris ° AS 29°8 
B Orio ; B8p — o9 
y Orio P B2 —28°5 


—29°7 . —27 

29°1 

yes 2 -— I 
~23°7 ‘ —30 
« C Maj . Br —27°8 
8 Velr . Ao 6-0 
B Cari ‘ Ao 1°7 
« Leon : B8 —11°8 
y Cent . Ao — 1°6 


—26 


Op COC fF ORNS 


—17 


« Virg : B2 —33°7 ne Pe —31 
B Cent . Bi —33°2 ge see et wide bie 
« Ophi ‘ AS 30°5 , ° +32 
e Sgtr : Ao 3°6 a bis whe ee fe 
« Aqil ‘ AS 41°6 8 , . +43 


« Pavo 2‘1 B3 —24'9 cue 
a PscA . A3 10°7 12 a ne ne an +16 
« Pegs : Ao ol 4 — 3°0 31 — 8 8 — 8 
S indicates the gradient as determined at Mount Stromlo; G at Greenwich; H by Hall; 
BC by Barbier and Chalonge. The unit is o-o1. 


A least squares solution of all these data was made and the results entered in 
column 4 of Table I, with approximate weights in column 5, unit weight corre- 
sponding to a standard error of +0-04. The zero, which of course is arbitrary, 
was adjusted to make the weighted mean difference Mount Stromlo minus 
Greenwich vanish. This adjustment fixes the zero of all the gradients in this 
paper. The Greenwich values in column 6 are from page 20 of reference (r). 
The figures after the decimal point in columns 4 and 6 have no real significance. 

From the residuals of the least squares solution the standard error of one 
comparison was found to be + 0-04, appreciably lower than the Greenwich figure 
of +007 or 0-08. The difference may well be attributed to the clearer atmo- 
sphere over the New South Wales table-lands. The unweighted r.m.s. difference 
Mount Stromlo minus Greenwich (col. 7) comes out as + 0-032 (allowance being 
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made for the zero adjustment), a very satisfactory figure in view of the s.e. of 
+ 0-02 quoted for the Greenwich standard stars. ‘The unweighted mean absolute 
difference is +0-025. 

The gradients in column 8 are taken from Hall’s photoelectric work (7), the 
values in the table being the mean of the Sproul values with 0-01 subtracted, and 
of the Amherst values with 0-02 added. Hall made these adjustments to bring 
his gradients into systematic agreement with those of Greenwich. The mean 
absolute difference Mount Stromlo minus Hall is +0-041; the mean absolute 
difference Greenwich minus Hall was +0-044. According to Hall 0-01 should 
be added to our gradients. 


Fic. 2.—Equialtitude comparisons between standard stars. 


Finally Barbier and Chalonge’s gradients (8), measured photographically in 
the region 43700-4600, appear in column 10. We have tabulated ¢,,—1-06, 
¢,, being Barbier and Chalonge’s 4, corrected by their Table VII, i.e. referred 
from a mean wave-length of A 4250 to a mean of A5000, and the 1-06 a constant 
which brings the French gradients into systematic agreement with both Greenwich 
and ourselves. The mean absolute difference Mount Stromlo minus Barbier 
and Chalonge is then +0-03. 

More complete comparisons with other observers are made later in the paper. 
As far as the standard stars are concerned, the general agreement must be 
considered good. 
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6. Programme Stars.—The programme was made up of most of the brighter 
southern stars of early type, enough northern stars to ensure a good overlap with 
northern observers, a number of later-type stars, and some fainter stars (down to 
5:I mag.) of interest. We planned to measure each programme star on at least 
three nights, against at least two standard stars. Actually only five stars received 
less than three full-weight determinations, and over half received four or more. 
As it was impracticable to make all programme star comparisons at equal altitudes 
we had to correct for atmospheric reddening. ‘The usual assumption was made, 
that the lower of two stars being compared appears redder by A (sec 2, —sec 2,), 
where z, and z, are the zenith distances of the two stars and A is a coefficient 
which has to be determined. A was found as at Greenwich, by comparing standard 
stars at different altitudes, groups of standards, well distributed in azimuth, 
being compared where possible. The values for A ranged from 0-25 to 0-40; 
on the best nights A was very near 0-30, and only on two out of 42 nights did it exceed 
0-40, 

In Table II are assembled all the relative gradients determined. The 
columns show in order the number of the star in the Yale Bright Star Catalogue 
(i.e. the H.R. number), its name, its R.A. and Dec. for 1950, its galactic latitude 
and longitude and its visual magnitude. Column 8 is the H.D. spectral type, 
and column 9 the type either as given in the Yerkes “‘ Atlas”’ (g) or as determined 
by Miss M. L. Woods at this Observatory from a classification on the Yerkes 
system of a large number of objective prism spectra taken by the late Mr W. B. 
Rimmer (10). ‘‘Atlas’’ types are italicized. Finally in column 10 is the mean 
gradient as determined at Mount Stromlo. 
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Ko (G8 V) 1°43 


As (A5 V) 0°30 
Fo Fo V 0°59 
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Spectral type 


H.D. 


Oe5 
Bo 
Bo 
Bo 
Bsp 
F8 
A2 
Bo 
B3 
Br 
Fo 
Ao 
B8 
Oe5 
Ao 
A5 
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F8p 
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Oe 
Oes 
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B8 
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Bip 
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B2 
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(Al V) 
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TasLe I1—Continued 

x Galactic : 

B.S. : A. Dec Spectral type Gradient 
No. %r 1950 to =a H.D. MKK G 


4798 a Muse —68 52 269 — 6 ‘ B3 B3 IV —o:28 
4819 y Cent —48 41 269 +13 . Ao Ao III —0°02 
4825 ~=—y Virg —or 11 268 ‘ Fo (FO V) 080 
4844 B Muse —67 50 270 — , B3 B3 V —0°24 
4853 B Crux —59 25 270 + : Br Bo III —0o'36 


4898* =» Crux —56 54 271 ‘5 B3/B3p B3 V —o'18 
5020 y Hyda —22 54 281 G5 (G65 IID) 1°72 
5028 « Cent —36 27 278 A2 A2V 0-08 
5056 « Virg —10 54 286 B2 (B1 IlI-IV) —0°34 
5132 e Cent —53 13 278 Br Br IV —0°30 


5190 v Cent —41 26 283 B2 B2 V —o'29 
5193 pw Cent —42 14 282 Bap B2 Ve —o'1o 
5231 ¢ Cent —47 03 282 Bzp B2IV —o'31 
5235 7 Boot +18 39 332 Go (G0 IV) 1°05 
5248 ¢ Cent —41 52 284 B3 B2V —o'26 


5267 B Cent —60 280 Bi Bx Il —0°33 
5354 t Lupi —45 286 + B3 B3 V —0o'22 
5440 7» Cent —4I 291 + B3p/Az2p <8 —0'23 
5459/60*« Cent —60 284 Go/Ks5 G2V 1°27 
5463 « Circ —64 283 ° Fo Fo VSr 0°49 


5469 « Lupi —47 290 ‘ B2 B2 Il —0°25 
5531 a? Libr —15 310 + , A3 A3 (met) 0°39 
5571 B Lupi —42 294 . Bap B2IV-V —o-25 
5576 «x Cent —4I 295 ° B3 B2 V —0o'29 
5671 y TriA —68 283 ° Ao Ao VEu 0°03 


5685 B Libr 320 . (B8 V) —0-'09 
5695 8 Lupi 299 : B2 IV —0o'31 
5708 e Lupi 297 . B3 V —o'22 
5776 y Lupi c 301 ‘ B3 V —0'24 
5812 7 Libr 309 3° B3 V —o'18 


5897 BTrniA 290 3° F2V 0°68 
5928 p Scor 313 . B2V —o'26 
5944 a Scor 315 + . (B2IV)  —0o-23 
5948 7 Lupi 307 . B2IV-V —o-29 
5953 8 Scor 318 , (Bo IV) —o'1l 


5984 B Scor 321 + ; (B0-5 IV) o’o1 
6084 o Scor . 319 ° (B1 IIT) 0°48 
6165 7 Scor 319 - ; (BO V) —0'27 
6175 ¢ Ophi 334 - . (09-5 V) 0°22 
6241 e Scor 317 + ° Ka III-IV 2°00 


6247* yp! Scor 314 + nek —0'26 
6252 p* Scor 5 314 4 ‘ B2 IV —0°30 
6262* ¢} Scor 311 aii 1°45 
6378 7 Ophi 335 + ° (A2 V) o'13 
6380 7 Scor 312 , . F2 Ill 0°86 


6453 @ Ophi 328 + ° (B21IV) -—o-21 
6462 y Arae 302 ‘ Bi V —o'1l 
6500 6 Arae 299 , B8 V —O'14 
6508 v Scor 319 . B2 IV —0'33 
6510 « Arae 308 , B3 Ve —0'20 


6527 A Scor 320 ‘ (B2IV) —0°34 
6553 @ Scor 315 ‘ Fo I-Il 0°86 
6556 « Ophi 4 + ° (Aé III) 0°30 
6580 « Scor 319 . B2 IV —0°32 
6615 t! Scor 318 ° Fol 1°15 
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TasBLe I1—Continued 
Galactic , ‘ 
; Dec. satiate Spectral type Gradient 
“ ee eas H.D. MKK G 


6 Arae : —50 06 311 —16 . Bip Br Il 0°03 
e Sgtr 3 20° —34 25 327 —12 . Ao Bg IV 0°04 
«a Tele 3° —46 316 —16 . B3 B3 V —o'2! 
¢@ Sgtr 3 “ —27 03 336 —12 ° B8 B8 III —0'07 


A Pavo . —62 14 —25 ° B2 B2 II le Var 


«x Pavo 3 51° —67 18 —27 ‘2V Fsp sas Var 

o Sgtr ° —26 22 7 —14 ° B3 (B3 IV-V) —0-27 
¢ Setr 59° —29 57 3 —17 ° A2 A2IV O°14 
¢ Aqil 3° 13 47 5 +2 Ao Ao V 0°05 
A Aqil 3% 58 35 — 7 3° Bog B8-5 V 10 


a7 Sgtr 8 06 , F2 (F2 IT) “71 
v Sgtr 8 O4 : f Aope ne 54 
« Aqil ° o8 - . aae 26 
a Aqil 3+ 44.0 (A7 V) “42 
5 Pavo 3° 3°5 Gs IV 43 


6 Adil 
a? Capr 
a Pavo 
e Delf 


B Pavo 


Bo IV —0'07 
G7 Ill 1°85 
B3 IV —0°25 
Bs V ool 
As Ill °° 


ovr n si 


w 


¢ Capr 3° 3 E ° ; G4 Ib 
B Aqar " _ 3 ° ; (GO Ib) 
8 Capr : - - . As (met) 
y Grus 50° — ° B8 III 
« Aqar 3° : ‘ ° ; (G1 Ib) 


Bs V 
B8-5 V 
A3 V 
Az Ill 
(43 V) 


(B9 V) 


x Grus 
¢ Pegs 
e Grus 
8 Aqar 
« Pisc A 
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NOTES ON INDIVIDUAL STARS 

897 @ Erid. Visual binary (8”). Mean gradient taken. 

3207 y® Velr. Wave-lengths A 4660 and A 5825, which coincided with emission bands, were omitted 
in measuring this Wolf-Rayet. 

4133 p Leon. Gradient according to Stromlo —o-10; Greenwich +0:°04; Barbier and Chalonge 

0:06; Williams —o-20. The gradient equivalent of Williams’ colours has been used, not 
the colours themselves. 

4730/31 a Crux. The bluest star. Mean of two components. 

4898 » Crux. Bright component only. 

5459/60 a Cent. Both components. Gradient for brighter estimated as + 1-20. 

6247 w' Scor. Eclipsing binary. No evidence for periodic change of colour. 

6262 © Scor. One of the reddest of the brighter B stars. 

7074 A Pavo. Magnitude and colour both variable. See remarks in Section 11 of this paper. 

7107 « Pavo. Cepheid. Colour range (from seven observations) about +0°80 to + 1°65. 

7852 « Delf. Gradient according to Mount Stromlo +0-’o1; Greenwich —o-12; Barbier and 
Chalonge—o:'23; Williams —o-10. The gradient equivalent of Williams’ colours has been 
used, not the colours themselves. 

8322 6 Capr. Metallic line star. Gradient equivalent to F2. 

8634 ¢ Pegs. Gradient according to Mount Stromlo +-0:04; Greenwich —o-11; Barbier and 
Chalonge +003; Williams —o-15. The gradient equivalent of Williams’ colours has been 
used, not the colours themselves. 

The following stars received less than three full-weight observations : 2282 ¢C Maj; 2451 
vPupp ; 245615 Mono; 2550 «Pict ; 2845 8 C Min. 
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Consideration of the residuals from these mean gradients gives an internal 
figure on their accuracy. We find that for stars with G less than about 0-80, 
the mean s.e. of a gradient is about +0-025, and of one determination about 
+0:05. (The corresponding Greenwich figures quoted by Greaves (11) are 
+0-03 and +0-06.) For red stars, with G greater than about 0-80, the mean 
s.e. of a gradient is about +0-04, and of one determination about +0-08. The 
shorter base line (0-7 of the full base line) accounts for most of the increase in 
errors of the red stars. 

7. Comparisons with other ‘ubservers.—Three other observers have measured 
the colours of the brighter southern stars—King (12), Wallenquist (13) and 
E. G. Williams (14). King used extra-focal photometry and three plate-filter 
combinations to measure Pg, Pv and Pr (photographic, photovisual and photo-red) 
magnitudes. His well-known work covers the whole sky, except for the Pr, 
where he measured only the bright northern stars. Wallenquist made a visual 
match of the colours of his stars against the colour of a lamp, varying the latter 
with a movable blue wedge filter. Williams employed a Fabry lens and three 
plate-filter combinations to give him magnitudes at the effective wave-lengths 
A 4100, 5150 and 6200. His programme was limited for the most part to early 
B stars south of — 30°. 

Comparisons between these colours and the Mount Stromlo gradients look 
very like the comparisons found by the Greenwich observers between their 
gradients and various colour indices (15). We have little to add to their discussion 
and will not describe our results in detail. ‘There are one or two anomalies, 
however, which should be mentioned. ‘The first, and most important, arises 
when our gradients are compared with King’s colour indices (Pv-Pg). ‘The 
agreement is good, except for the bright stars « C Maj and « Cari, both of which we 
make appreciably bluer than does King. The possible systematic effect so suggested 
is discussed in detail in the next section. ‘There is a good deal more scatter when 
the gradients are plotted against Wallenquist’s “reciprocal temperatures”’ ; 
in particular, our colours for « Grus and 7 Ceti are appreciably bluer than Wallen- 
quist’s, and for @ Scor, 4 Scor and « Agar appreciably redder. The correlation 
between our measurements and Williams’ is good, especially for the red (space- 
reddened) stars. The mean curve agrees well with that calculated from the 
effective wave-lengths quoted above, and, excluding a few stars with spectral 
peculiarities (y? Velr), noticeable discrepancies exist only for p Leon and ¢ Pegs, 
for which Williams’ colours may be less accurate. Both these stars are discussed 
in the notes to Table IT. 

The very close relation between gradients and King’s red indices (Pr—Pg) is 
worth mentioning. It is illustrated in Fig. 3. Most of the stars being in the 
north, we have had to use Greenwich gradients for more than half. «C Maj 
again stands off, as does 5 Leon in the opposite sense. ‘The Greenwich gradient 
for the latter is well confirmed by other spectrophotometric, observers. 

Turning now to the northern spectrophotometric work, we have made 
extended comparisons of our colours with those of Greenwich, Hall, and Barbier 
and Chalonge, all of whom were mentioned in connection with the standard 
stars. ‘The comparisons are shown graphically in Figs. 4 and 5. 

The mean absolute difference Mount Stromlo minus Greenwich for 23 stars in 
common is +0-053, appreciably larger than the +0-025 found in Section 5 for 
standard stars alone. ‘The agreement can nevertheless be considered satisfactory. 
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These stars are far south at Greenwich; some do not rise more than 40° above the 
horizon there, and the Greenwich observers may well have had difficulty with 
extinction corrections. By the same token they are well north for us, and while 
the Mount Stromlo atmosphere should give less trouble, some of these stars, 
because of their relatively short observing season (curtailed for a number by a 
run of bad weather), were not observed as often as we would have liked. The 
worst discrepancies occur for p Leon (Stromlo minus Greenwich = — 0-14), e Delf 
(+0°13) and { Pegs (+0-15). As it happens, our gradients for these three stars 
have quite good weight. All the five Orion stars we have in common with 
Greenwich we make bluer, by 0-04 in the mean, a trend in which we are supported 
by Stebbins and Whitford (16), who find 5, « and £Orio among the bluest stars 
which they have measured. And the comparison with the 23 Greenwich stars, 


as opposed to that with the seven standard stars only, indicates a correction to our 
zero of +0-0I. 


King’s 
red index 
(Pg-Pr) 

1-5 











! 1 I 
0-0 0-5 1-0 1-5 2-0 
gradient ——> 





Fic. 3.—Gradients and King’s red indices. 


For the 17 stars we have in common with Hall we find a mean absolute 
difference +0-038. This is ‘smaller than the mean absolute difference —0-044 
for Greenwich minus Hall. According to this comparison our zero should be 
increased by 0-017. 

‘The comparison with Barbier and Chalonge shows good agreement except for 
a few rather large residuals and the comparatively big correction of —0-06 
indicated to our zero-point. The latter was unexpected after the good agreement 
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shown with Barbier and Chalonge in Table I for the standard stars. Closer 
examination shows that while the Barbier and Chalonge gradients for these are 
all of high weight (except perhaps for « Pisc A), this is not the case with many of the 
stars shown in Fig. 4, especially as less weight is given to those south of —15°. 
Perhaps over-correction by the French observers for atmospheric extinction for 
these southern stars is the reason for the discrepancy. 
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—0-8 








j l | 
—0-5 0-0 0-5 1-0 
s-—> 





Fic. 4.—Comparisons among gradients determined at Mount Stromlo (S), Greenwich (G), and by 
Hall (H), and Barbier and Chalonge (BC). Crosses are G-BC. 


Of the stars which stand off the main curve, A Orio is of very early type(MKKO8) 
and various measurements of the colours of « Delf disagree rather widely. (See 
note to Table II.) ¢Ophi is a strongly reddened star ; its position on the figure is, 
as was first suggested by Greenstein (17) (see also references (8), (11) and (16)), a 
consequence of the fact that interstellar reddening, as measured by gradient excess, 
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increases with increasing wave-length, so that a reddened star is redder at 45000 
than at 44250. xLeps is an Fo supergiant and 5Capr a well-known metallic 
line A. Beyond values of G of approximately 0-6 the correlation deviates from 
linearity, in the sense that the colour temperature then begins to increase with 
wave-length. A few Greenwich stars, indicated by crosses, have been added to 
illustrate the point. « Hyda and z* Orio, both apparently normal, are now seen 
to stand off. But for these and the other stars mentioned in this paragraph, 
except ¢ Ophi, the Barbier and Chalonge gradients are of low weight, with s.e.’s 
of the order +0-08, so that too much importance should not be attached to these 
anomalies. 


























4 z 
Magnitude 





Fic. 5.—Comparisons between Mount Stromlo (S), Greenwich (G), and Hall (H) gradients, and 
King’s colour indices (K), to test for a possible magnitude affect. 


8. The Gradients of Sirius and Canopus.—According to these comparisons our 
real errors are not much in excess of the + 0-02 or + 0-03 calculated from the internal 
agreement of the individual observations, and hence no appreciable systematic 
effects were expected. It was therefore disconcerting to find that the two 
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brightest stars on the programme, Sirius (m,—1-4; AI V; G=-—o-16) and 
Canopus (—0-7; Fo I; +0-25), were both much bluer than their spectral types 
would indicate, the more so as the comparison with King showed either that we 
were making these stars too blue or that King was making them toored. Further, 
the determination, by Dr Woolley and the writer (4), of 1-06 for the gradient of the 
Sun relative to Sirius implied that, if the gradient of Sirius is — 0-16, that of the 
Sun is 0-90. Not only is this much too blue for a G2 V, but we are then also in 
difficulties with the stellar temperature scale. 

The suggestion, that we are making the colours of the bright stars systematically 
too blue, is clearly serious. ‘To examine it further, we plotted in Fig. 5 the colour 
differences Mount Stromlo minus Greenwich, Hall and King (King’s colour 
indices were multiplied by 1-5) against magnitude. It can be seen that, apart 
from « Cari and « C Maj, there is no evidence for a magnitude effect in our colours. 

In Table III we have assembled data concerning the brightest stars on the 
programme. ‘This includes the gradient as determined at Mount Stromlo and, 
anticipating Section 9, that expected from the spectral type. Apart from the 
two brightest stars, the differences found are small. With «Cent both the 
measurements and the estimation of G(.Sp) are complicated by the presence of 
the K-type companion.  Orio must be slightly reddened, as one would expect. 
We make it bluer than does Hall (cf. Table I), but as it is a super-giant and also 
shows traces of emission, this difference is possibly a real indication of anomalies 
in its continuous spectrum. «Crux is certainly a very blue star, the bluest we 
have measured; but in view of its position as the nearest Br in the sky, situated 
moreover in a very clear region of the galaxy, this is not altogether unexpected. 


Tasie III 
The Brightest Stars 
mM» MKK type G(m) G(Sp) G(m)-G(Sp) 
—1'4 A1lV —o'16 0°05 —o°'21 
—o'8 Fo la 0°25 o°55 —0°30 
0°3, 1°7 G2 V+ KoV 1°27: 1°35 —o08 : 
3 B8 Ia —o-ol —o'lo 0°09: 
o'4 Fs IV 0°88 0-90 —o0'02 
°'6 B2 IV —0o'24 —0'25 ool 
o'7 Br II — 0°33 —0°35 0°02 
o'8 A7V 0°42 0°42 000 
o'9 Br IV —o'4I —0°35 —o'06 
G(m) is the gradient measured on this programme. 
G(Sp) is the gradient estimated from the spectral type. 


It is certainly very difficult to account for the presence of a systematic error 
of this kind, and even more difficult to see how such an error should not affect 
the fainter stars, making them progressively redder. As reference to Table IV 
will show, some of our bluest stars were fainter than 3-0. And the good corre- 
lation found between our gradients and spectral type (Fig. 6) is a strong argument 
against systematic effects. 

Two recent measurements of the colour of Sirius are now available *, although 
with Canopus the only colour we have, beside King’s, is Wallenquist’s, which 
neither disagrees with nor confirms ours. Danjon (18) has published visual 


* Barbier and Chalonge’s gradient for Sirius is omitted; they state that the gradients for their 
three brightest stars (« C Maj, 8 Orio and « Lyra) are unreliable. 
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measures of eight bright stars, in which he makes « C Maj definitely bluer than 
B Orio, and may be said to support our value. Seven of the eight stars are on 
Hall’s programme, with whose figures Danjon agrees well. Eggen (1g) has 
determined an accurate photoelectric colour of — 0-07 for Sirius, on a scale which 
is near to, though not identical with, the International photographic scale * ; 
—0-07 is near the International colour for spectral type Ar. 

To the writer, Dr Eggen’s result is conclusive so far as Sirius is concerned. 
Meanwhile we can only repeat that, apart from these two stars, the evidence is 
against the presence of systematic effects, and we shall proceed on the assumption 
that they are in fact absent. Further observations of Sirius and Canopus are 
planned. 

g. The Spectrum-Colour Relation.—The preponderance of bright B stars in the 
southern sky has been known for many years, and one of the objects of the present 
work was to find if the colours of these stars, a fair proportion of which might be 
near enough to avoid appreciable space-reddening, would help to fix more 
definitely the blue end of the spectrum-colour relation. It is the bluest stars 
which are of interest here, and a list of these is givenin Table IV. Five are stated 


TaBLe IV 
The Bluest Stars 


Distance 
(parsecs) 
3447 o Velr B3 V 3°6 ° 70 
4199 6 Cari Bo Vp 2°9 ‘ 120 
4656 5 Crux B2 IV 29 , 130 (B) 
4730/1 a Crux Br IV °'9 . 50 
4853 B Crux Bo III 1°4 . 120 (B) 
5056 a Virg Bi III-IV 1‘2 90 
5132 e Cent Br IV 2°4 -o° 180 (B) 
5231 ¢ Cent B2 IV 3°0 160 
5267 B Cent Bx II o'7 125 
5695 5 Lupi B2 IV 3°3 ° 180 (B) 
5262 pe? Scor B2 IV 3°7 200 (B) 
6508 v Scor B2 IV 29 —0°33 135 
6527 A Scor B2IV ° —0°34 go 
6580 « Scor B2IV . —0°32 115 
G 39 y Pegs B2V . —0°30 80 
G 46 7 Auri B3 V : —0'22 70 
G 113 7 U Maj B3 V ° —0'23 40 
G 221 B Ceph Br IV . —0'30 180 
For stars with numbers prefixed with a G the gradients are from the Greenwich list (2) 


the number being the Greenwich number. Distances followed by (B) are from Blaauw’s 
parallax for the Scorpio—Centaurus cluster (20). 


B.S. Star Type My 


by Blaauw (20) to be “ certain”’ or “ probable’ members of the Scorpio—Centaurus 
cluster; his distances are given for these stars and are distinguished with a B. 
The distances of the remainder are from Miss Woods’ luminosities, using the 
luminosity calibration in the Yerkes “‘ Atlas”. The possibility of a correction of 
+0°OI or +0-02 to our zero should be kept in mind, but will not affect any of our 
general conclusions. The four stars at the foot of the table are the two bluest 


* The relations between the two scales will be described by Dr Eggen in a forthcoming paper in 
the Astrophysical Journal. 
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(y Pegsand 8 Ceph) and the two nearest early B stars (7 U Maj and » Auri) observed 
at Greenwich. In Fig. 6 we have plotted gradient against Yerkes spectral type, 
using only the nearest stars ofeach type. For Fo and earlier only Mount Stromlo 
gradients were used, and for B3 and earlier only those of Table IV. Later than 
Fo we included a number of Greenwich stars; these are not distinguished in any 
way. As one would expect from Morgan and Bidelmann’s work (21), the 
correlation is extremely good. Indeed, earlier than about Go the entire scatter 
may well be attributed to observational error. The one exception is y Virg, an 
Fo V only about 100 parsecs away, in galactic latitude 61°. The Mount Stromlo 
gradient of 0-80 agrees well enough with Hall’s mean of 0-74. On its colour this 
star could well be F3, or even F5. 





1 T 








J i j 
BO AO FO KO 
Yerkes Type 





. 6.—Gradients and Yerkes spectral types. Filled circles represent dwarfs, crosses giants and 
open circles super-giants. No differentiation is made for stars earlier than FO. 


At about Go the dwarf-giant differentiation can be observed setting in, 
in the customary way, though here as elsewhere one cannot be sure 
that space-reddening has been entirely eliminated. Space-reddening is in 
fact a major difficulty throughout this subject, well illustrated by 7 U Maj and 
n Auri. Though the two nearest in a table of early B stars, they are definitely the 
reddest. (We must not forget, however, the possibility of dispersion in the 
colours of the early B stars.) It is surprising that a starin such a low galactic latitude 
as « Crux should be so blue, and groups such as y, A and « Scor must be in exceed- 
ingly transparent regions of space. (Yet o Scor, just as near—it is also in the 
Scorpio—Centaurus cluster—is a Br with a gradient of 0-48.) Stebbins and 
Whitford remark that the situation of « Virg, in galactic latitude 50°, might make 
it the bluest star inthe sky. If this is so, then all the stars in Table IV are sensibly 
unreddened, and the intrinsic gradient of a Bo is no bluer than about —0-35. 
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Atthe same time the present work makes the early B stars at least as blue as — 0-35, 
and possibly as blue as, say, —0-40. ‘This question of the intrinsic colours of 
B stars is important for space-reddening and other problems, and it is of interest 
to find the equivalent of our colours in other systems. For instance, it has been 
suggested by Oort (22) and others that the normal C1 colours of Stebbins and 
Whitford should be corrected by —0-04 (so that, for example, the Cr for an 
unreddened Br should be —0-26 and not —0-22). In Fig. 7 Cr colours (taken 
from M.W.C. 621) are plotted against Greenwich and Mount Stromlo gradients. 
In choosing the line 

G=0°76 +4:30CI 


to represent the relation, more weight has been given to the points at the blue 
end of the curve than to the intermediate ones; this seemed the better procedure 
for the present particular application. We then find that for G=—o-40, 
C1 = —0-27 and for G= —0-35, €r= —0-26. Our work clearly supports Oort. 
Not enough observational material exists at present to enable the corresponding 
figures to be derived for other colour systems, in particular for the International 
and photo-red colour indices. 


cl t T 





0-00 


—0-30 - 








i a. i i 
— 0-40 —0-20 0-00 0-20 0-40 0-60 0-80 
Gradient 
Fic. 7.—Stebbins and Whitford’s C1 colours plotted against Greenwich (dots) and Mount 


Stromlo (crosses) gradients. Early-type stars only are included; the red stars are space- 
reddened blue stars. 





In this general connection it is worth pointing out that different observers 
do not always agree on what are the bluest stars. Thus Stebbins and Whitford 
include among their half-dozen bluest V-I stars 10 Lacr (Greenwich G = —0-15), 
) Orio (G = — 0-21), 8 Ceph (G = — 0-30) and ¢ Pers (G= —0-23). Of these only 
the short-period variable 8 Ceph has a really blue gradient, while both Greenwich 
and Mount Stromlo make 6 Orio slightly but definitely reddened, as one would 
expect from the strength of its interstellar K line. Barbier and Chalonge 
make 10 Lacr very blue. These differences are outside the limits of observational 
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error, and are perhaps to be explained in terms of the deviations from black-body 
conditions noted by Stebbins and Whitford, and apparent in Fig. 3a of reference 
(16). 

10. The Stellar Temperature Scale.—According to the preceding section the 
intrinsic colour of a main sequence Bo-B3 can be represented on the Greenwich 
system by a relative gradient of —0-35. Ifthe colour temperature at Ao is 16,000, 
that at Bo will then be 38,000 deg., which may be compared with the admittedly 
approximate value of 57,000 deg. found theoretically by Unséld for + Scor (29). 

Stebbins and Whitford remark on a difficulty they experienced in adopting 
a temperature scale which would allow Ao and Go to agree with observation and 
at the same time not make Bo too hot. Thus if they assume 16,000 deg. for Ao, 
their (_V —I) values give them temperatures in the region of 100,000 deg. for Bo, 
And if, as our work suggests, and as they themselves realize was possible, even 
their bluest stars were slightly reddened, this difficulty is accentuated. But if 
gradients or (V —R) colours are used the difficulty is not so apparent, and it is 
instructive to see how it comes about. 
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Gradient (G) 
Fic. 8.—Gradients plotted against Stebbins and Whitford’s V-R and V-I colour indices. 


Open circles are space-reddened stars, and the full lines represent the relation between colours on 
each pair of systems and of a series of black bodies of differing temperatures. 


In Fig. 8 we plot G—o-67(V—J) and G—o-g0(V—R) against G for the 
fifty-odd Stebbins and Whitford stars for which gradients are known. It will 
be recalled that V, R and J are heterochromatic magnitudes measured over fairly 
narrow spectral regions of effective wave-lengths 44220, A7190 and A10,300 
respectively. In Table 5 of (16) Stebbins and Whitford have computed (inter alia) 
the colour indices (V —R) and (V —J) as functions of black-body temperature, 
the colours having a common zero at T=5500deg. The corresponding relation 
between gradient and temperature is set out in Table V of (x), These figures 


3 
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enable us to draw the full lines of Fig. 8, relating in each pair of systems the colours 
of a series of black bodies of varying temperature. The zero has been chosen 
so that T=6500deg. corresponds to G=1-15, (V—J)=—0-25 and (V—R) 
=-—0-21. This change of zero point necessitated the addition of 0-32 to the 
tabulated (V —J) and of 0-21 to the tabulated (V—R). From the plot of (V —J) 
against G it is clear that departures from black-body conditions are present *, and, 
bearing in mind the more linear relation which evidently exists between (V — R) 
and G, it seems safe to attribute the greater part of these to the I region, in particular 
to the Paschen discontinuity. Hall and Williams (28) have shown that the dis- 
tribution of energy in the spectrum of Vega (Ao) is quite closely that of a black 
body up to the Paschen limit, beyond which it shows an excess of energy of about 
0-23mag. And Hall (7) has shown that this behaviour is typical of stars in the 
region of Ao, showing, as one would expect, close analogies to the Balmer 
discontinuity. 

Somerelevantfiguresare collected in Table V, which should be self-explanatory. 
The (V —J) colour for Bo is rather arbitrary but cannot be far out. In one group 
we have adopted the common temperature of 6500deg. for Go, in the other 
30,000 deg. for Bo. The better agreement between (V—R) and G is clear. 
The lower temperature which the photoelectric colours attribute to Go must be 
due to heavy line absorption in the V region; this would not affect gradients. 

To the writer these considerations argue rather tellingly in favour of photo- 
graphic methods for absolute temperature work. In this problem the lower 
accuracy obtainable is not as important as the ability of the photographic plate 
to isolate spectral regions a few angstroms in width. This, of course, is far from 
. decrying the great accuracy and convenience brought by the photocell to many 
problems of stellar colours. 


TABLE V 


Variation of colour temperature with spectral type according to gradients and to two photoelectric 
colour indices 





G V-I V-R G V-I V-R G V-R 





(deg.) (deg.) (deg.) | (deg.) (deg.) (deg.) 
I‘I5 —0'25 —0'20 6,500 6,500 6,500 6,350 6,000 6,150 
0-00 —1°85 —1°45]| 16,000 15,000 17,500 | 15,000 12,400 14,500 

—0o'35 —2°60 —1°85| 36,000 80,000 50,000 | 30,000 30,000 30,000 




















11. Individual Stars of Interest.—We will conclude this paper with a few remarks 
on individual stars of interest. 

First y* Velr, the brightest Wolf-Rayet star, has as common proper motion 
companion the apparently normal B3 y! Velr. Not only can we determine the 
absolute magnitude of y* Velr from the spectroscopic parallax of y1, but also we can 
use the colour of y! to correct that of y? for space absorption. We find 7? rather 
bluer than y!; intrinsically it must be a very blue star (24). 

The following table contains data relative to four of our Be stars, selected 
because they appear to lie in very transparent directions of the galaxy. 
205G Cari is quite near in space to, and at about the same distance from the 


* As otherwise there would be a random scatter of points about the solid line. 
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Sun as, @ Cari, of gradient —0-32. Cent is similarly only one degree away 
from v Cent, of gradient—o-29. 6Cent is in a very transparent region, and 
a Arae is included as the bluest of the Be stars. For none of these stars would 
we expect space reddening to be appreciable, and taken as a group we consider 
that they present strong evidence in favour of the intrinsic reddening of Be 
stars (25). 


TaBle VI 
Type Gradient Distance 
(parsecs) 
4140 205 G Cari Bs IVe —o'08 
4621 8 Cent B2 Ve —or'lr 80 
5193 » Cent B2 Ve —o'ro go 
6510 a« Arae B3 Ve —0'20 60 


Finally we have A Pavo, an apparently normal B3 which, after the style of 
y Cass, developed emission characteristics at some time in the last 20 years (26). 
As the accompanying journal of measurements shows, the colour of this star 
appears to vary. No obvious variations in emission-line intensities have been 
noticed. 


Tasie VII 
Measurements on A Pavo 


Date Gradient 


27/7/48 “eer 


1/8/48 0°07 
3/8/48 0°04 
12/10/48 o'18 
24/10/48 o'2I 
29/10/48 016 
23/6/49 —o'lo 


Cousins (27) and others have suspected this star of variability. 
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IC 4406: A DOUBLE NUCLEUS PLANETARY NEBULA 
David S. Evans 
(Received 1949 October 5) 


Summary 


IC 4406 is found to consist of an elongated mass of gas under the 
gravitational and radiative influence of two stars. The results of a photo- 
metric and spectroscopic investigation are presented. ‘The system does not 
appear to form an equilibrium configuration and seems to resemble closely the 
type of model much discussed in connection with the theory of the origin of 
the solar system. It is hoped that the observations reported may be of interest 
to theoretical workers. 





In the course of a survey of bright southern planetary nebulae attention was 
attracted to the system IC 4406, which, because of its unusual astrophysical 
interest, seemed to merit a detailed study. Dreyer* describes the system as 
‘Planetary, stellar, 10 mag., extended in position angle 80°’’. Campbell and 
Moore + describe it as being of about 20” diameter, slightly flattened on the 
east and west edges. These workers obtained three spectra showing no more 
than the N, and N, lines, and deduced a velocity of —22 km./sec. A drawing by 
Wilson shows an undifferentiated area of luminosity of the size and shape described 
above. ‘There is no indication that a central star was observed. In the drawing 
a star is shown about 33” to the west of the centre of the nebula. Fleming { 
gives the magnitude of the nebula as 9-5. 

The following plates of this object are on file at this observatory :— 


TABLE I 
IC 4406: R.A. 14" 19™-3; Dec. —43° 55’ (1950.0) 
Galactic Long. 288°, Lat. + 15° 
Date Emulsion Exposure Filter Calibration Observer 
(min.) 

1949 April 6 1.A.Z. 10, 3 we ADT 
April 19 L.A.Z. 15 a> DSE 
April 19 1.A.Z. 10 wi DSE 
April 19 K 103a-O 15 wie DSE 
June 21 1.A.Z. 10, 5 ai ~ DSE 
July 21 K 103a—O 15 DSE 
July 21 K 103a—O 15 DSE 
July 22 K 103a-O 5 at DSE 
July 22 K 103a-E 15 W 25 DSE 
July 23 K 103a-E 20 W 25 DSE 
July 24 K 103a-E 10 W 25 DSE 
July 24 I. HP3 15 W 12 DSE 
July 24 I. HP3 15 W 12 DSE 
July 26 I. HP3 40 W 12 DSE 
Aug. 23. K 103a-O 40° W 18a DSE 
Aug. 23K 103a-E 20 W 25 _ DSE 

I.A.Z.=Ilford Astronomical Zenith; K=Kodak; I=Ilford; W=Wratten. 
W 25: transmission, 6000 A. and to red. 
Filters 


-++ -f 


{+ +. +f 


l+++4 


W 12: transmission, 5200 A. and to red. 
W 18a: transmission (10 per cent +) 3300-3900 A. 





* J. L. E. Dreyer, Second Index Catalogue, Mem. R.A.S., §9, 105, 1908. 
+ W. W. Campbell and J. H. Moore, Lick Observatory Publications, 13, Part 4, 1918. 
t Mrs Fleming, Harvard Circular, No. 60, 1901. 
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These plates were calibrated in the same way as those discussed in a 
previous paper.* 

In the 74-inch telescope the object appears visually as a small almost rectangular 
blue smudge extended in an approximately east-west direction, sharply bounded 
to north and south, but somewhat blurred to east and west. No central star 


P(W) 





F (E) 
Fic. 1(a).—Blue. 


can be seen, but there is a prominent star to the west, lying apparently well outside 
the system. 


On short-exposure red and blue plates, the central star (star C) is visible 
faintly. It appears to lie at the centre of a structure of oval shells, having an 
* M.N., 109, 94, 1948. 
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axis ratio of about 1-5:1, the longer axis lying north and south, and the overall 
north-south diameter of the largest shell being about 34”. A small, rather 
asymmetrical dark area surrounds star C. The nebula has east-west symmetry 
about the central star, and on short-exposure plates its eastern and western ends 
appear slightly concave. On yellow plates star C can scarcely be seen. 

P(W) 


__30 
oer ga 
Ne salt 


———_—- 








Fic. 1 (6).— Yellow. 

As soon as the exposures are increased, star C disappears and the nebula 
can be traced as extending to greater and greater distances to east and west, 
remaining closely symmetrical about the position of the central star. The 
eastern and western ends remain slightly concave, until the exposures are sufficient 
to trace the parts of the nebula surrounding the second star (star S) lying 33” 
to the west of the centre. A shell structure surrounding star S can now be seen, 
and this has the effect of making the western end of the nebula concave: at this 
stage the eastern end is somewhat irregular and slightly concave. Finally, it 
is possible to trace still greater faint extensions at both ends. To the west star S 
is surrounded, while on the east there is a rounded end. ‘The whole nebula is 
shaped like a slightly curved elongated oval approaching a bean shape. In 
addition to this general structure there are a number of details to be described. 
On the east, short dark bays run inwards, while on the west, the form of the shell 
structure surrounding star S includes an asymmetrical tongue of brightness, 
issuing from the nebula and running close round the star to the south-west of it. 
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To the north-east of star S lie two small and slightly irregular fan-shaped dark 
areas centred on the star. The structures shown on the red, yellow and blue 
plates appear to be essentially the same, and it is believed that the material 
available does reveal the full extent of the system. 

The surface of the nebula has been photometered, and the slightly smoothed 
version of the surface intensity variation which the method produces is illustrated 
in Fig. 1(a), (6) and (c). Fig. 1(a) shows half-magnitude isophotes for blue 
plates, A369* and A374, uncorrected for background sky fog. Similar plots 
for yellow plates A 378 and A 380 and for red plates A376 and A377 are shown 
in Figs. 1(5) and 1 (c). For the red and blue plates, allowance for sky fog involves 
renaming the isophotes according to the following scheme :— 

Taste II 
Blue Red 
Plates A 369*, A 374 Plates A 376, A 377 
Sky-fog level : 2°43-+0°045 Sky-fog level : 3°93+0°030 
Uncorrected Corrected Uncorrected Corrected 

2°20 4°01 3°50 4°72 

1°70 2°48 3°00 3°60 

1°20 1°63 2°50 2°84 

0°70 0°95 2°00 2°20 

0°20 0°35 1°50 1°62 
—0°30 —o°2I 1°00 1°08 
—o-80 —o-74 

It will, of course, be appreciated that both these magnitude scales are 
arbitrary as to their zeros. 

There is surprisingly little difference in the appearance of the system in light 
of various colours, except that the region near star C appears somewhat bluer 
than the periphery of the nebula. The red plates probably represent the distri- 
bution mainly in Ha light, the blue plates mainly in N, and N, and possibly 
A3727 light. Fig. 2 shows a comparison, taken along the mid-line of the nebula, 
of the corrected brightnesses in red and blue light. The positions of stars C 
and S are shown by the arrows. The interesting features are: firstly, that the 
differences are relatively small; secondly, that the plot is closely symmetrical 
about star C and seems hardly affected at all by the presence of star S; thirdly, 
that the variation, especially in the blue, is not far off an inverse square law. 
By arguments similar to those advanced in a previous paper * we may probably 
attribute the nebular light to radiation originating in star C and scattered by the 
surrounding gas. 

At my request Dr Thackeray was good enough to obtain a spectrum of IC 4406 
with the Newtonian spectrograph, which gives a dispersion of 3-7 mm. between 
Hf and He. The slit was placed longitudinally on the nebula, and the following 
lines were found :— 

Tas_e III 
Spectrogram NS 106 
1949 Aug. 25. Exp. 25™. Slit 10. Emulsion R 55 
Line N, N, HB Hy H8 He+3967(NelID) 3869 (Nell) 
Strength 40 12 6 4 I 5 5 
Length 55 oS ep 16 19 21 
(seconds of arc) 


The line A 3727 unfortunately would fall precisely at the edge of the film 
and is therefore not recorded. 


* Loc. cit. 
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The magnitude of star S could be estimated from image diameters on plate 
A317 by comparing it with plate A328 of selected area No. 132, taken by 
Dr Thackeray on April 23-24 on the same emulsion (K 103a—O), with exposure 
times of 2™, 6™ and 18™. With an error of at most two to three-tenths of a 
magnitude, star S is estimated as having a photographic magnitude of 14:7. 
Star C shows two or three magnitudes fainter than this, but it is clearly much 
subject to absorption, and is hardly brighter than the nebular surface. 

Plates taken at different wave-lengths show that while the transition from 
red to ultra-violet markedly alters the appearance of the surrounding star field, 
it hardly changes at all the appearance of the nebula and star S. It is clear that 
star S is of a very early type, and is hardly, if at all, redder than the nebula itself. 
No other information concerning the colour of star S is available at present. 

To summarize, we have the following picture of the system: it is a mass of 
gas of a fairly uniform colour, of an elongated oval shape, measuring about 
100” x 37”, having a star, presumably of type W, but much dimmed by absorption, 
at its centre of symmetry. Star S, of early type and apparent photographic 
magnitude about 14:7, lies at a distance of 33” to the west of the central star, 
also on the mid-line of the nebula and within the nebulosity. 

With the means at our disposal at present we cannot get an estimate of a 
spectroscopic parallax for star S, so that the distance and scale of 1C 4406 are 
unknown. However, we can consider the various possibilities which can be 
summarized in the following table :— 


TaBLe IV 
Star S Distance (parsecs) Distance CS (parsecs) 
Hypothetical WM 
—5'3 16 
—0'3 1 6 
+4°7 0°16 
Of these three guesses, the first makes star S intrinsically very luminous and 
puts it at a distance which, according to present ideas of galactic structure, should 
take it right outside the galaxy past its centre. There is no particular objection 
to the system IC 4406 lying at a distance greater than the galactic centre, since 
the latitude is moderately large and so is the longitude difference from the centre. 
On the other hand, if there is a window here we should expect to find extra-galactic 
nebulae appearing on the plates. Even on plates which should have a limiting 
magnitude of, say, 18, no nebulae can be found. Thus, if there is no spatial 
absorption and even if IC 4406 lies to one side and past the galactic centre, we 
must regard the first guess as giving a distance at least three times the possible 
reasonable maximum. This argument suggests that an upper limit for the 
distance C'S is about 5 parsecs. The third guess, according to the classical 
‘| Hertzsprung-Russell diagram would make star S of far too lateatype. Guessing 
at a maximum possible value for M of, say, +1-7, the distance comes out at 
2500 parsecs, giving CS=o-4 parsecs. Thus, assuming no absorption, we can 
probably place the distance C'S as lying in the tenfold range from 0-5 to 5 parsecs. 
The discussion of absorption effects falls into two parts. Either the system 
is very distant and lies clear of the central absorbing layer of the galaxy, or it is 
relatively near and lies within this layer. In the first case the absorption will be 
constant and may be taken as 0™-25 cosec 15° or almost exactly I magnitude. With 
assumed values for M of —5-3 and — 0:3, values of the distance, D, then come out 
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as 6:3 x 10* parsecs and 6-3 x 10% parsecs. The corresponding distances above 
the galactic plane, Dsin15°, are then 1-63 x 10* parsecs and 1-63 x 10% parsecs, 
with C'S being 10 parsecs and 1 parsec in the two cases. 

The first estimate is probably much too large, for it makes IC 4406 an object 
lying well outside the galaxy and, just as in the unobscured case already discussed, 
we should be able to detect other extra-galactic objects on the plates. The second 
estimate is more reasonable. In this case IC 4406 lies more or less on the boundaries 
of the galaxy. It seems that we must retain as our rough upper limit for the 
distance C'S the value of 5 parsecs already deduced for the unobscured case. 

If however we assume a smaller absolute luminosity for star S, the system 
IC 4406 becomes a member of the galaxy, and the correct way to allow for the 
absorption is by reckoning a dimming of so many magnitudes per kiloparsec. 
The following table shows values of absorption computed on various 
assumptions :— 

TABLE V 
oM-s5 per kiloparsec oM-75 per kiloparsec 


fice Total se oto Total 
D Dsin15° CS absorption DD Dsin15 CS absorption 





M=-—2°'3| 6100 1580 0°98 3™1 4800 1240 o'77 3™-6 
M=—0'3 | 4000 1030 0°64 2™-o 3250 840 0°52 2™-4 
M=-+1°'7 | 2300 590 0°37 ym-2 2000 520 0°32 1™+5 








The last entry represents probably the low limit for the value of CS. The most 
probable values would seem to be D = 4000 parsecs, CS =0-65 parsecs, but there 
is a wide range of possibilities. 

It is not the purpose of this paper to attempt a theoretical discussion, but it is 
possible to advance certain general physical arguments which may provide a 
pointer to the nature of the system. What seems of greatest importance is that 
although star S is involved in the system, the distribution is closely symmetrical 
about star C. Ifwe consider a very artificial model consisting of two stationary 
point masses, m and m’, surrounded by a tenuous cloud of gas, then eventually 
the gas will conform to the equipotential surfaces defined by these two masses. 
These are readily sketched out in bi-polar coordinates as the surfaces 


m/r +m’ /r’ = constant. 


For various values of the ratio m/m’ they give a series of dumb-bell-shaped curves 
centred about the mass-centre of the two points. Some forms can be obtained 
which are reminiscent of IC 4406 but the centre of symmetry does not coincide 
with one of the point masses. Symmetry about m can only be secured by taking 
m/m' large, but then the surfaces approximate to spheres centred about m. 

A physical possibility is that the two stars form a binary of long standing. 
In this case, neglecting radiation pressure, the equipotential surfaces will be 
somewhat modified, but still, in a general way, it is possible to maintain that 
the equilibrium form would have its centre near the mass-centre of the two stars 
and not coincide with one of them. 

General considerations of radiation pressure seem to present the same 
dilemma. ‘That is, if star C contributes all the radiation-pressure effects, it 
would be expected that the nebula would have spherical symmetry about this 
star. We know from the observational results that it is probable that star C 
does contribute most of the radiation exciting the nebula. One possibility still 
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remaining, of accounting for the state of the system as one which is of long standing 
in its present form would be to attribute all the radiation-pressure effects to star C 
while sharing out the gravitational effects between the two stars. It seems 
however that there are difficulties about this interpretation too, for the joint 
effect would be to leave the centre of symmetry lying between C and S, and, 
moreover, according to the mass-luminosity rule it is a highly artificial hypothesis 
to make star S both relatively faint and massive. 

I am inclined to think that the solution of the dilemma lies in the hypothesis 
that IC 4406 is a system which has not yet reached equilibrium and which cannot 
therefore be of long standing. In fact, IC 4406 may represent a case of a stellar 
encounter. On this view star C would originally have been the central star of 
an ordinary planetary nebula having symmetry about the central star. Star S 
then produces tidal effects, which should give at any instant ellipsoidal equi- 
potential surfaces centred about star C. This is, in fact, what we do find. Moreover, 
although there is a rough north-south symmetry about the line CS there are 
important detailed structures which show asymmetry. The whole nebula is 
slightly curved as if there were a component of transverse relative velocity in the 
plane of the sky. Moreover, the shell structure round star S shows the asym- 
metrical bright tongue already described, which suggests a modification of an 
originally spherical shell system due to relative motion. 

In view of the possibility of changes I have consulted early descriptions of 
this object. Referring to the period 1898-1901 Delisle Stewart *, using the 
Bruce telescope at Arequipa, gives the description, “‘two nebulous stars make 
dumb-bell, remarkable”. Innes, writing in 1902, describes his observations 
with the Cape 7-inch:—‘“‘Is a fine planetary, 10 mag., about 20” in diameter. 
Examined on the same night with the 18-inch it appeared dumb-bell shaped. 
With the 24-inch object-glass prism it was found to be a gaseous nebula. On 
Carte du Ciel Plate 3689 with 1 hour exposure are shown two very elongated 
spindle shaped nebulae of the same length, parallel to each other and in contact 
at their points of greatest condensation or brightness (angle of.elongation 80°)’. 

There is a good deal which is puzzling about these descriptions. Even allowing 
for differences of scale and focal ratio it seems doubtful whether the nebula as 
now seen could fit Stewart’s description. If it were not for the discrepancy of 
appearance in Innes’ observations, to which he himself calls attention, one might 
be confident in saying that the appearance had changed, but the variety of the 
descriptions given casts some doubt on their reliability. ‘The most that can be 
said is that there is a suspicion of change. 

The investigation suggests that IC 4406 may represent a stellar encounter, 
and the system is clearly one which will repay not only an expert theoretical 
discussion but also more detailed spectroscopic examination when this becomes 
possible. 

Acknowledgments.—I am much indebted to Dr Thackeray for the use of 
some of his material and for his critical comments on this paper, and to Dr Knox 
Shaw for comments on the manuscript. 


Radcliffe Observatory, 
Pretoria : 
1949 September. 
* D. Stewart, Harvard Annals, No. 60, p. 163 
+ R. T. A. Innes, M..N., 62, 469, 1902. 
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Montuty Notices of R.A.S. 


Fic. 3.—IC 4406. 


Date: 1949 April 19; Emulsion: Kodak 103a—0; Exposure: 15 minutes; Magnification: xX 10 
from original plate, giving 2:25 seconds of arc per mm. on this plate. 


D. S. Evans, IC 4406: A Double Nucleus Planetary Nebula. 





Montuty Notices of R.A.S. Vor. 110, Prate 2. 


He comparison 














(a) Proxima Centauri, 1949 Fuly 4 (NS 45). 

(b) Proxima Centauri, 1949 August 3 (NS 78). 
(c) RR Telescopti, 1949 August 8 (NS 89). 

(d) RR Telescopii, 1949 October 7 (NS 121). 

(e) 47 Tuc. Var. No. 1, near maximum (NS 97). 


A. D. Thackeray, Five Southern Stars with Emission-line Spectra. 





FIVE SOUTHERN STARS WITH EMISSION-LINE SPECTRA 


A. D. Thackeray 


(Received 1949 December 28) 


Summary 


Low-dispersion spectra showing bright lines of Proxima Centauri, 
RR Telescopii and of three long-period variables in 47 Tucanae are described. 
Proxima Centauri is of dMe type, the emission lines apparently varying in 
intensity. RR Telescopii underwent a rather abrupt transition from F-type 
absorption to emission plus continuous during the period 1949 August to 
September. 





1. The Newtonian spectrograph at present in use on the 74-inch Radcliffe 
reflector is best suited for low-dispersion studies of faint extended sources such as 
diffuse nebulae, globular clusters and extra-galactic nebulae. In the absence of 
a spectrograph with higher dispersion it has also been used, chiefly at times of 
bright moon, for obtaining miscellaneous stellar spectra likely to be of interest. 
The dispersion (240 A./mm. at Hy) is too low for accurate spectral classification, 
but the instrument has proved useful for detecting emission-line spectra with 
which the present paper deals. 

All the spectra here described were taken on R55 Orthochromatic film 
covering the range 5700 to 3750A., with the exception of two of RR Telescopii 
(at its brightest) for which a rather slower emulsion, Ansco Plenachrome, was 
used. The f/1 Schmidt camera gives a reduction factor of 4:2 and hence a wide 
slit could be tolerated. 

2. Proxima Centauri.—Table I lists nine spectroscopic observations of this 
star which, as far as the writer is aware, has not been hitherto accessible to a slit 
spectrograph. In this and subsequent tables in the “‘ Remarks” column, (a) 
denotes absorption, (e) emission. 


TABLE I 
Observations of Proxima Centauri 


Date Exp. Slit 

JD 243+ (min) (mm) Rementes 
3088-301 60 0°22 TiO+ MgH(a); no emission 
3102°297 75 20 H(e)!; Cate)! 
3130°210 25 20 TiO strong; no emission 
3132°203 30 15 Hy, Cat(e) 
3133°224 30 "15 Hye); Ca, Mn(a) 
3135°222 30 "15 Hy, Ca*(e) 
3137°205 30 13 ? Hye); some moon-fog 
3139°202 32 13 No emission 

103 3153°210 30 Ors No emission 


The spectrum is of advanced dwarf M type with bands due to Ti0. Comparison 
spectra were obtained of 20C 923 and 20C 1143 (dM5 and dM3 respectively 
according to Joy’s list of standard types recommended tothe I.A.U.). The MgH 
band at 5211A., found strong in dwarfs relative to giants by Ohman*, has also 


* Y. Ohman, Stockholm Ann., Bd. 12, No. 8, 1936. 


POSTS WAL TAP MA ag ED OCC 











46 A. D. Thackeray Vol. 110 


been measured. A rather strong diffuse structure extending to the violet from a 
measured wave-length 5567 +7. has been found on three spectra with correct 
exposure in this region; it is believed that the (1,2) band of MgH with P and Q 
heads at 5568 and 5550 may be responsible. It is perhaps surprising that it 
appears stronger than the (0, 0) band at 5211, although the latter is masked by 
the strong TiO band at 5167. Ohman has also found absorption in this region in 
M dwarfs extending farther to the red and he has attributed it to the (0, 1) MgH 
band with P and Q heads at 5621 and 5609. 

Line absorption due to Ca4226 (apparently asymmetrical as found by 
Lindblad in dwarfs), Mn 4030 and Fe5269 has been measured; Cat H and 
K absorption is suspected on a weak continuous background on NS 82. 

Emission lines due to hydrogen and Ca+ H and K have been measured on 
five of the nine spectra. This is a common feature of late M dwarfs according 
to Joy *, who finds that late type rather than duplicity of the stars appears to be 
the determining factor. The emission lines appear remarkably strong on NS 45 
(1949 July 4); see Plate 2(a). In this instance Hf, Hy, Hd, He +H, K, HZ 
and possibly two further Balmer lines appear bright, whereas on the other spectra 
Hy, H and K were usually the only bright lines detected. ‘The longer exposure 
given on this occasion compared with later spectra was probably more than 
offset by the fact that the later spectra were obtained with a freshly silvered 
mirror. The density of the continuous spectrum in the violet is much the same 
as on some of the later spectra with weak emission; see Plate 2(b). Two 
alternative conclusions suggest themselves: either (1) Proxima had undergone 
an abrupt physical change on 1949 July 4 perhaps analogous to that of the flares 
of Luyten’s dwarf L—726-8 f, or (2) the origin of the bright-line spectrum does 
not coincide with that of the M-type continuum and a relatively greater proportion 
of the former region may have been transmitted by the slit on this single occasion. 
Unfortunately, the dispersion is too low to permit adequate measurement of 
radial velocity of the emission relative to the absorption or to the known velocity 
of the «Centauri system; further observations at the Cassegrain focus are 
urgently needed to settle this interesting point. 

No accurate magnitude observations of Proxima are as yet available at the 
Radcliffe Observatory to test whether the star is liable to outbursts like those of 
L-726-8. The Cape observers have examined 58 plates of the star taken for 
parallax between 1932 and 1949 and have found no evidence of any variation in 
brightness exceeding 0-10 mag. However, the possibility of variations of about 
0-2 mag. lasting less than an hour, such as have been observed tf recently for 
B.D. +20° 2465, is not excluded. 

3. RR Telescopti.—This star was at one time suspected to be of U Geminorum 
type, but Kukarkin and Parenago’s 1947 Catalogue lists it as “‘semi-regular”’, 
varying between 12:5 and less than 14 with a period of 387 days. Dr P. Kirchhoff 
of Johannesburg discovered that the star had brightened to 6-0 in 1948 July and 
photographs at the Cape Observatory showed that the brightening had taken 
place at least as early as 1946 September.§ 


* A. H. Joy, Ap. F., 105, 96, 1947. 

t+ W. J. Luyten, P.A.S.P., 61, 179, 1949. 

¢ K. C. Gordon and G. E. Kron, P.A.S.P., 61, 210, 1949. 

§ R. P. de Kock, M.N.A.S.S.A., 7,74, 1948. Harvard records show that the star was at mag. 10 
and brightening towards the end of 1944 (M. W. Mayall, H.B., No. 919, 16, 1949). 
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The star was placed on the observing list for the Newtonian spectrograph at 
the Radcliffe Observatory and spectra obtained on 1949 June Ig and 20 showed 
pure absorption apparently of supergiant F type with H and K three times as 
intense as Hd. The only other element detected in absorption was 77+ (five 
lines), though Si+ 4130 was suspected on both films. HB could not be detected 
but the strength of the other members of the Balmer series in absorption suggested 
that HB was lost due to a balance of absorption and emission as in some Be stars. 
This explanation receives some support from the next spectrum of 1949 August 8 
which shows Hf bright, Hy undetectable, while Hd and higher members of the 
series are fairly strong in absorption; see Plate 2(c). In addition to HB, Fe+ 5018 
and 4923 appear bright. As on June 20, 7i* is well represented in absorption 
and Fe* absorption is conspicuous by its absence. The next two spectra (taken 
on 1949 September 12 and October 7) show a remarkable change to a rich bright 
line spectrum superposed on a continuous background in which no absorption 
is measurable; see Plate 2(d).* The dense September exposure shows the 
Balmer lines bright from Hf to H 3750 at the limit of the film; Ca+K is also 
bright with intensity approximately equal to H3889. Nearly all the remaining 
lines are due to Fe+; multiplets 27, 28, 32, 37, 38, 42, 48 and 49 of the “ Revised 
Multiplet Table” are all represented, lines with laboratory intensity equal to 
or greater than 4 appearing with very few exceptions. 7+ 3913 and 3900 are 
weakly present; other strong 77+ lines in the blue are lost in blends with Fer. 
Further spectra taken on November 1 and December 4 show no marked change, 
though there has perhaps been some fading of the bright lines relative to the 
continuous. 

TaBLe II 
Observations of RR Telescopii 
Date Exp. Slit 
JD 243+ (min.) (mm.) 
31 3087-480 O'1s H, Ca*, Tit(a); supergiant F spectrum 
37 3088-45 "15 H, Ca*, Ti*(a); supergiant F spectrum 
89 3137°429 13 HB, ? Fe+(e): Cat, Ti*(a) 
III 3172°365 "10 H, Fet+, Cat(e)!! 
121 3197°267 , “15 H, Fet, Cat(e)!! 
131 3222°315 ° “10 H, Fet*, Cat*(e)!! 
132 3255°266 ° O'1s H, Fe+, Cat(e)!! 


NS No. Remarks 


According to information supplied by Dr Kirchhoff the star began to fade 
in July. The transition to a pure emission spectrum is thus explicable in terms 
of a weakening of the continuous backgroynd which up till October 7 at least 
left the emission within the line cores relatively unaffected. The prominence 
of Ti+ in the absorption phase and of Fet+ in the emission phase recalls the 
differential behaviour of these two elements in the various shells of Nova Herculis. 
The spectra of RR Tele are of too low dispersion to state whether or not the emission 
lines are accompanied by violet absorption as in the usual nova, but the long 
period of maximum light and gradual fading which now seems to have set in 
suggest that the star should be classified as a “‘slow nova” with small velocity 
of expansion. The examples of Nova Pictoris and Nova Herculis show that the 
star merits the attention of double-star observers. 

* This film shows very irregular trailing of the star along the slit due to the shortness of the 


exposure for such a bright star; but the reader should have no difficulty in distinguishing the stellar 
lines from the helium comparison. 
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4. Long-Period Variables in 47 Tucanae.—This well-known globular cluster, 
NGC 104, is remarkable for its apparent lack, despite searches, of typical short- 
period variables. Bailey’s investigations revealed eight variables, apparently 
either irregular or of long period, of which numbers I, 2, 3 have been shown by 
Shapley * to have definitive periods of 211, 203 and 192 days respectively and 
maxima of between I1-3 and I1-5 mag. 


Tasie III 
Observations of 47 Tucanae Variables 


Date Exp. Slit 
JD 243+ (min.) (mm.) 
3144°463 45 ors He)! 
3152°439 25 "15 H(e); star near maximum 
3172°401 22 “15 H(e); star near maximum 
3174°474 40 ‘10 Underexposed 
3187°433 40 "15 
3197°332 30 "15 H(e) 
3197°389 60 O'1s 


NS No. Var. No. Remarks 

97 
105 
112 
116 
120 
122 
123 


~~ Ne ew We 


All three of these stars have been photographed near maximum light with 
the Newtonian spectrograph (see Table III) and have been found to be of Me 
type, with 7:0 bands of only moderate strength, considerably weaker than in 
RCent (M4e). Hy and Hé have been found bright on all spectra with adequate 
exposure. The two lines are in most cases almost equal in strength, but the first 
two spectra of variable No. 1 taken at or perhaps a little before maximum show 
H8 somewhat stronger than Hy; see Plate 2(e). NSg7, a fairly dense spectrum, 


shows four Balmer lines beyond He; He itself has not been observed on any of 
the spectra, presumably due to the usual obscuration by Cat+H. In all these 
features the three variables appear to be quite typical of long-period variables 
of Me spectra, despite the fact that according to Shapley they have absolute 
magnitudes at maximum some 5 magnitudes brighter than the average such star. 
(It is of interest to point out that if this estimate is correct Plates 2(a) and (e) 
correspond to stellar spectra, superficially not very dissimilar, but differing by 
20 magnitudes in the absolute brightness of their sources.) Attempts to follow 
variable No. 1 some months after maximum, when other bright lines might be 
expected to appear if the analogy holds for the phases of decreasing light, have 
not as yet met with success. 


Radcliffe Observatory, 
Pretoria : 
1949 December 13. 


* H. Shapley, P.N.A.S., 27, 440, 1941; Harv. Repr., No. 228, 1941. 





SECULAR CHANGES OF STELLAR STRUCTURE 
AND THE ICE AGES 
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Summary 


A sudden increase in the rate of nuclear energy generation by a certain 
amount leads to expansion and to a decrease of the luminosity by nearly an 
equal amount (if measured in stellar magnitudes). Theoretical considerations 
are confirmed by actual calculations of expanding stellar models. 

The corresponding changes of structure consist of two components :— 
(i) homologous changes, in which the release or absorption of gravitational 
energy is proportional to T/u; (ii) pure changes of structure at constant 
potential energy, in which a certain energy of transition, w, is released in the 
central regions and absorbed in the outer regions, or vice versa. An approxi- 
mate formula for w is given in terms of the variation of central density and 
temperature. 

For homologous changes the time of relaxation, to, of a disturbance is 
inversely proportional to s+3, s being the temperature exponent of nuclear 
energy generation; for the Sun fy is about 5 x 10° years with s~15; this falls 
well within the order of magnitude of the duration of the Quaternary glaciation. 

A scheme of complex expansion, produced by a sudden increase of the 
intrinsic energy generation, is worked out by numerical integration; pure 
changes of structure and the boundary condition are schematically taken into 
account. As compared with homologous expansion, the drop in luminosity is 
delayed and reduced in depth. 

Non-static changes of solar structure are suggested as a plausible cause of 
the recurrent ice ages on the Earth. An increase in the circumcentral 
hydrogen content in the Sun, produced by turbulence suddenly starting, may 
lead to an increased energy production and thus to a decreased luminosity. 

Other existing theories of Quaternary climatic changes are briefly 
reviewed; it is shown that none of these is able to account for the observed 
size of the changes. 





1. Introduction.—It has been formerly suggested by the writer* that a 
temporary expansion of the Sun may have caused a dimming of solar luminosity 
and a corresponding deterioration of the terrestrial climate known as the still 
lasting Quaternary glaciation, similar causes possibly being responsible for other 
ice ages which occurred on the Earth several times since the Archaean. After 
some disturbance in the energy balance the Sun is supposed to start re-adjusting 
itself toward a new equilibrium configuration; if this re-adjustment requires 
expansion, gravitational energy is absorbed at the expense of the radiative energy 
flux and the apparent luminosity is likely to decrease. The time intervals 
involved in such a re-adjustment agree more or less with the known duration 
of the Quaternary. 

In the present paper a closer theoretical consideration is given to the above 
suggestion. This involves the theory of non-static stellar models.¢ As the 


*E. Opik, Tartu Obs. Pub., 30, No. 3, 1938; Acta et Commentationes Univ. Tartu (Dorp.), 
Ser. A, 33, 1938. 

t F. Hoyle and R. A. Lyttleton, M.N., 107, 246, 1947, discuss non-static stellar models from 
a standpoint which, however, is hardly related to our present problem. The present theory was 
practically ready several years ago, but circumstances of the war prevented its publication. 
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time intervals involved in the changes considered here are large as compared 
with the period of free pulsations of the star (which is about a couple of hours 
for the Sun), the deviations from hydrostatic equilibrium are insignificant; 
therefore, all the ordinary equations of hydrostatic and radiative equilibrium 
remain practically valid in these non-static models. The only difference with 
respect to a static model consists in the distribution of the energy sources, 
gravitational energy sources being additionally active. 

2. The Model and the Disturbance.—We will consider a typical stellar model, 
with the intrinsic (nuclear) energy sources confined to the central core (not 
necessarily a point-source); such a model consists of three main portions*: 
a convective core of radius r,, mass M,; an intermediate radiative equilibrium 
shell between r,, M, and r,, M,; and a peripheral convective region, extending 
from r., M, outwards until the photospheric boundary is reached at the radius 
r=R and total mass M of the star. In the core, the hydrogen content is supposed 
to decrease with time, as the result of nuclear reactions; in this way models with 
non-uniform chemical composition are formed. Several scores of different 
typical cases of such models have been calculated by the writer.t Cowling’s 
model f is a particular case of uniform composition throughout and of a vanishing 
peripheral convective region. The different composite models are generally 
non-homologous with respect to one another. 

In a real case, the type of the model is determined by the photospheric 
boundary condition (p,, T,); the density p, at photospheric level depends upon 
the law of photospheric opacity; the effective temperature T, must satisfy the 
condition of radiation to space equalling the heat transfer from the interior; 
the photospheric level itself may be practically identified as the level at optical 
depth §. The solution of the equations of hydrostatic, radiative and convective 
equilibrium is made to satisfy the boundary condition. 

The boundary condition is extremely sensitive to the extent of the peripheral 
convective region, whereas for an equal extent of the latter it depends very little 
upon the conditions in the deep interior. Therefore we will symbolize the 
boundary condition by its practical equivalent, or by a certain ratio M,/M, 
which is the relative mass inside the second (subphotospheric) interface. Cowling’s 
model, with a “‘ complete radiative outer shell’’, corresponds to M,/M=1; from 
what we know of the photospheric opacity, the Sun’s case probably does not 
differ much from that. Another extreme case is M,/M = M,/ M, or when the two 
interfaces coincide, or the ‘‘complete adiabatic model’’. Most real cases of 
dwarf stars may be supposed to fall between these two extremes. 

The boundary condition cannot be violated very much; but a slight increase 
in the extent of the peripheral convective region leads to a considerable decrease 
in the photospheric temperature and radiation, which means that a considerable 
fraction of the heat generated inside cannot be radiated any more to space; the 
extra heat becomes stored inside in the form of potential energy of gravitation 
and leads to expansion; as a consequence, the central temperature decreases, 
the nuclear energy sources are slowed down, and equality between energy 

* E. Opik, Tartu Obs. Pub., 30, No. 4, 1938. 


+ Ibid., also 31, No. 1, 1943 and recent material, in press or unpublished. 

t T. G. Cowling, M.N., 96, 42, 1935. 

§ Indirectly the boundary condition is also extremely sensitive to the conditions in the interior 
in so far as these affect the extent of the peripheral convective region; if, however, the latter is kept 
constant (by a suitable choice of the energy flux), the boundary condition varies little. 
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generation and radiation to space is approached again. ‘This process of self- 
regulation of the energy balance is favoured by the circumstance that a decrease 
in the internal energy flux produces automatically a shrinking of the peripheral 
convective region, as an inherent property of the solution of the equations of 
radiative equilibrium; thus the boundary condition meets half-way the self- 
regulation of the energy sources.* 

The disturbance of the energy balance can be imagined as a suddenly initiated 
convective mixing in the central regions of the star. The writer has calculated 
modelst+ in which the hydrogen content increases outwards, starting from the 
boundary of the core, and which after continued transformation of hydrogen 
into helium are apt to develop a new convective region outside the core. The 
latter leads to a discontinuous change of structure, the local polytropic index 
changing by a finite amount from a value of m=1-2—1-3 to the adiabatic value 
nearm=1'5. The discontinuity is caused by the circumstance that the gradient 
of molecular weight which has installed itself during the radiative equilibrium 
stage is wiped out by the newly established convective equilibrium. Because 
of the discontinuity, the convective disturbance spreads out with finite speed, 
producing mixing of the material over a wide but limited range in the interior; 
a merging of the new convective zone and the original core is possible. In any 
case the result of such a disturbance will be a transport of hydrogen toward the 
exhausted central regions and thus a possible increase of the energy output. 
The detailed study of the processes involved is made elsewhere (unpublished). 
A surprising result is that the time intervals involved in the propagation of the 
convective disturbance are not so short, up to 107 years for solar dimensions, 
with major effects showing up after the lapse of some 10° years since the start. 
However, for the general way of treatment applied in the present paper the 
details are irrelevant. It suffices to admit the possibility of a ‘‘sudden’’} increase 
in the intrinsic energy output. ‘The consequences of this increase can be studied 
in a differential way, by comparing the disturbed model with an undisturbed 
one, both models being built upon exactly the same set of differential equations 
and physical constants but for the changed distribution of the energy sources. 
For the differential study almost any type of a radiative equilibrium model is 
well adapted, but for the sake of simplicity we will use for this purpose models of 
uniform composition (constant molecular weight throughout), namely the 
writer’s series of models numbered from I to 7§ which are representative of 
the static composite models of the entire range of boundary condition from 
M,/M=M,/M (No. 1) to M,/M=1 (No. 7)/|, the last case differing from 
Cowling’s 4] but slightly (in taking into account radiation pressure and electron 
scattering). 

3. The Luminosity of Slowly Expanding or Contracting Models.—The 
integration of the equations of radiative equilibrium of a static model can be 

* For boundary condition cf. also H. N. Russell, M.N., 97, 127, 1937, and L. Biermann (jointly 
with T. G. Cowling), A.N., 267, 131, 1938. 

+ Proc. Roy. Irish Acad. A, §3, 1, 1949; Armagh Obs. Contr., No. 3, 1950 (in press). 

} Sudden as compared with a time-scale of the order of 10* years. 

§ Loc. cit. 

| Outside these extreme values of M,/M the boundary condition may still vary, leading to 
“* super-radiative ” photospheres at M,/M=1 (early-type stars), or to “‘ super-convective ” structures 
with M,/M=M,/M (perhaps red dwarfs) in which the convective transport of heat does nowhere 


vanish ; for stars around the solar mass these cases hardly do occur. 
q Loc. cit. 
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performed by assuming an arbitrary value for the intrinsic energy generation, 
E (core-source); each value of E leads to one definite boundary condition and 
vice versa. Actually the range of E (for a given mass of the model) which leads 
to a more or less reasonable boundary condition is rather narrow. In nature 
the boundary condition (although unknown to us) is fixed by the behaviour of 
photospheric opacity and, therefore, the value of E is also precisely fixed, E = Ey 
being the static value leading to the prescribed boundary condition. In the 
static model (core-source) the radiative flux outside the core is constant and 
equal to Q=£); this is also identical with the luminosity Ly= Ep, representing 
the energy radiated to space. A deviation of E from the equilibrium value Ep, 
at the imposed boundary condition, is incompatible with a static solution; it 
requires an expanding (E> E,) or contracting (E<E,) model. For the non- 
static model let QO, be the variable radiative flux in the radiative equilibrium 
zone, H,—the gravitational energy generated inside r per unit of time; evidently 


O,=E+H,. (1) 


Let Q,, H, correspond to the boundary of the convective core, Q,, Hz to the top 
of the radiative equilibrium zone, and let H be the total gravitational energy 
generated by the model, so that its luminosity is given by 


L=E+H. (2) 


Consider a small deviation E — E, from the static value; the principle of continuity 
and average quantities require that, at the same boundary condition, in the 
non-static model a certain average value of the radiative flux intermediate between 
Q, and Q, be equal to the static value Ep, 


QO, < Fy <Q2. 


Ey=kQ, +(1—k)Q3, (3) 
where & is a positive fraction. In the solution of the equations of radiative 
equilibrium a variation of the radiative flux affects the boundary condition the 
more the deeper it is situated; therefore the relative weight of Q, in (3) exceeds 
the weight of Q,, k>1—k or k>4; the semi-empirical value k = § fits well actual 
stellar models as shown below.* 

Substituting Q, and Q, according to (I) into (3), setting 4, = H,/H, h,=H,/H 
and with 


This may be written also as 


a=kh, +(1—k)hy, (4) 
it is easy to obtain 
E-—E,=-—-aH; (5) 
also, substituting H from (5) into (2) we get the deviation of the luminosity of 
the non-static model from the static value as 


L—Ly=L-E,= - @—9 (£- By) (6) 


As h, and A, by their definition are reasonably positive fractions, the characteristic 
factor a is by (4) also a positive fraction. Hence (5) tells us that an increase of E 
requires a negative value of °7, which means expansion, and vice versa. From (6) 
we infer that the change of luminosity is always of opposite sign to the change 
of the nuclear energy generation. 

* Cf. Addendum. 
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Assuming with Eddington that the release of gravitational energy per unit of 
mass and time is proportional to the temperature (actually to 7/u for non- 
homogeneous composition), an assumption strictly corresponding to homologous 
changes of structure but which turns out to be an excellent approximation for 
the actual non-homologous changes likely to occur, we find that for Models Nos. 7, 
6, 5, representing possible configurations of a solar mass with a relative mass of 
the core about 0-14, about 22-23 per cent of the gravitational energy is released 
inside the core. We assume the round value A, =0-25, further 4, =1 and k=§; 
this yields a=} and from (5) and (6) 


H=-— 2(E - E,), (5’) 
L~ Ey =(8=<E). (6’) 


The expected change in luminosity is opposite in sign and equal in amount to 
the change in the energy generation, and the release of gravitational energy is 
the double of that. The value of a may vary with the type of the model ; however, 
for widely different models the above equations are changing but insignificantly. 
These general considerations are substantiated by actual calculations of 
expanding models. To the usual equations of hydrostatic and_ radiative 
equilibrium the differential equation of the release of gravitational energy is 

added: 
dH,=CTdM,; - (7) 


the constant C, negative in the case of expansion, is chosen arbitrarily; E is then 
determined from an integration over the whole convective core according to 
equation (I): 


1, 
Q,=E+C| TaM,. (8) 


0 
Q,, the radiative flux at the boundary of the core, is given independently by the 
assumed structure of the core (polytrope m = 1-5) and the law of opacity ; continuity 
of the radiative flux at the boundary of the core is assumed. 

The comparison between a certain set of non-static and static cases for the 
same boundary condition is given in Table I. The non-static models of this 


TABLE I 
Expanding and Static Models at Same Boundary Condition and Constant Central Temperature 
M,/M 0°297. -0°386)=—soggr 0682 0824 0929 @8©0°969 1°000 
Exp. 0°297 0°297 0°295 0289 0'284 0'279 390 0"277 0'277 
MMs 3. ‘ 
Stat. 0°224 0188 0°165 oO'rsi O°144 o'140 0°139 0138 


Exp. I'0O0O)=—s«dI'000—'—«édT‘000—s«édar'014~—s«1048 105 1°185 1°423 
Stat. 1*000 1*000 1009 1‘037 1086 -186 1°286 1°576 


P Exp. 6-00 6:03 6-10 6-60 7°58 9°19 11°49 19°96 
PciPm*) Stat. 6:00 6-06 6°50 7°69 9°5 13'0 16°9 31°3 


E 
—2°5 log —0'36 —o'59 —0o'79 —0'88 —0'95 —0°95 —0'95 
0 


L 
—2°5 log - +o° +o-78 +091 +0°97 +0°99 +1°01 +1°02 
0 
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table have been built all on the same core, corresponding to zs =1-5 (Emden’s 
reduced radius), all the range of M,/M being obtained by varying the parameter C 
in (7); the results were reduced to constant mass and central temperature with 
the aid of homologous transformations (using mean values of the radiation 
pressure); the values for the static models are interpolated according to M,/M 
from the actually calculated set of Nos. 1-7. In the table, M,/M is the relative 
mass of the core, R the relative radius of the model, p,/p,, the ratio of central to 
mean density. The last two lines of the table give the change of the energy 
output and of the luminosity, both expressed in stellar magnitudes; from a com- 
parison of these two lines it appears that, for the present finite and rather large 
changes, the simple relation (6’) deduced in the present section is of course 
no more valid, but that a similar relation referring to the logarithmic changes is 
approximately fulfilled: the logarithm of the luminosity varies by about the 
same absolute amount but in the opposite direction as compared with the logarithm 
of the energy production. For the first two cases of small changes the 
infinitesimal relation (6’) is, however, almost exactly fulfilled : 


M,/M =M,/M 0-386 0-491 
(E-E,)/E, +0°20 +0°39 +0°72 
(L—L,)/E, —O-21 — 0:39 —0O°51 


4. Changes of Structure.—Any change of stellar structure can be considered 
as consisting of two components: a general homologous transformation in which 
the potential energy changes inversely as the radius; and a non-homologous 
transition at constant potential energy which we will call further a “‘ pure change 
of structure’’. 

In usual notations the potential energy of a spherically symmetrical body is 
given by * 

_M 
Q=G | =<, (9) 


~0 


7M 
Q=3| PVdM, (10) 


~0 


with 
RT 
PV = — =(y-1)c,T, II 
5 (11) 
where y =c,,/c, refers to the mixture of gas plus radiation. When the molecular 
weight is constant and when 1 —f is small, an average value of 8 over the whole 
model can be used and (10) may be represented as 
2K — 
Q= TM, (12) 
Bp 


. TdM 
P= tt. (13) 





* Cf. A. S. Eddington, Internal Constitution of the Stars, Cambridge, 1926, and various 
other sources. 
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For a true polytrope of index m (9) is reduced to 


1 oe 
ad ac" te abe (14) 


Table II represents a sequence of pure changes of structure along true 


polytropes; for these B=const. or B=1 is to be assumed somewhat artificially ; 
for n=3°25 the data are based on Chandrasekhar’s* calculations, the rest on 
Emden’s.t According to (12), Q=const. implies T=const., whence the 
relative central temperature is determined from 
7 R=" as 
T M'(n+1) 
and the relative central density is 


po~(5—n)9.P ; 
Pm 

R’ and M’ are Emden’s reduced radius and reduced mass (Eddington’s 
notations). | 

The table shows that a pure change of structure with increasing m requires 
an increasing central density and a decreasing central temperature. Both 
these changes mean a release of heat at the centre. As the potential and the total 
energy of the model remain constant, an equal amount of heat must be absorbed 
in the outer regions. Thus a pure change of structure requires a certain ‘“‘ energy 
of transition’’, w, to be transferred from the central regions to the outer regions 
(increasing n, positive w) or in the opposite direction (decreasing n, negative w). 


TaBLe II 


Polyiropes at Constant Mass and Potential Energy 


n R Ab oe Pp. 
relative relative relative 
“0000 1°884 1*0000 “000 
-1667 1‘806 09582 *196 
*4000 1°758 0°9332 462 
"7500 I°710 0‘9080 ‘690 
"0000 1°693 08989 1°836 
*5000 1°661 08815 2°423 
35°0000 1°655 0°8784 3°632 


>RWWNDND 
Sonoudun 


This remains true for all kinds of pure change of structure, whether polytropic 
ornot. Asshown below, the energy of transition for the composite non-polytropic 
models of the above-described type is considerably smaller than for true polytropes 
and amounts to a small fraction of the potential energy, practically of the 
order of 10-%, 10~? at the utmost. Therefore, instead of using the laborious 
method of numerical integration in each particular case, a differential formula 
can be derived. The formula yields the energy of transition in units of the potential 
energy with sufficient accuracy (to quantities of second order) for transitions 
between two arbitrary models when the variations of T., p, at constant potential 
energy are given. 
* S. Chandrasekhar, Ap. ¥., 89, 116, 1939. 


+ R. Emden, Gaskugeln, Leipzig, 1907. 
t A. S. Eddington, op. cit. 
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The energy released per unit mass in a change of state is 
dqg= —(PdV +c, dT). (16) 


Assuming 6 =1, P=p, (gas pressure), y= $ (monatomic gas), an approxi- 
mation admissible for stars of about the solar mass and even larger, with the 
aid of (11) and (12) it is easy to obtain with p=1/V, when for dq, dp, dT the small 
finite quantities g, Ap, AT are substituted 


QT (2 Ap ) 
=} ~{- — —- —}, I 
q FG; T (17) 


which gives the heat released per unit mass in the finite but small change of state. 
For homologous transitions Ap/p and AT/T are constant, whence qT; for 
non-homologous changes this is no more true, but the absolute value of g remains 
a quantity of the order of 7; we may assume 
fp &F 


q a — HF wa bx, 


where x= M,/ M is the relative mass and 


Be. Se, (18) 


On a broad line u = T/T, does not very essentially differ from 1 —x, as shown by 
the run of these quantities for the composite model No. 7: 


I—x 1-000 0-862 0°624 0-407 0°295 O-I2I 0-000 
u I'000 = -0'793-—Ss(0'606—-0-473. Ss «000-0277 —: 0000 


This is typical for all models which may come into question. ‘Therefore (17) 
may be represented fairly well as a quadratic wave 

Q T, 

=}— =(1-—x)(a—bx). 

I= 39 REI ) 


For pure changes of structure 


~0 


-M ‘1 
| g@M=M | qdx 
~0 


must vanish, which yields 6 = 3a. 
The energy of transition represents the maximum absolute value of 


“z 
M)| qdx 
“0 
which is attained evidently at 
q=9, 2 == > =}; 


b 


% separates the region of positive g from that of negative g. Performing the 


integration, further assuming as a broad average 7\,/T =1-7 (cf. Table II, around 
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n= 3-3°25), and substituting (18) for a, we get an expression for the energy of 


transition 
w= 9(24% AT) 


8\3p 7. 


which can be used irrespective of the type of model. ‘Table III contains collected 
data for the seven typical composite models.* 

The effective value of the polytropic index in the sixth line of the table is 
chosen to satisfy the potential energy in equation (14); as compared with the 
true polytropes of Table II, at constant potential energy the range in radius 
corresponds more or less to n (effective); the range in 7, and p,, however, is 
much smaller than for true polytropes of an equal range in m; moreover, there 
is an “inversion” in the trend of T, between No. 6 and No.7. As aconsequence, 
the energies of transition between these core-source composite models are 
considerably smaller than for true polytropes. 


(19) 


TaBLe III 
Potential Energy and Energy of Transition of Composite Models 


[4=const. } 
Model No. 2 3 4 
M,/M 0°358 0478 0684 


M,/M 0189 0167 = or150 


Pel Pm : 6:06 = 642 7°70 


M=Mo 
R=Ro 


n, effective 1°50 1°51 1°55 1°68 1°84 2°18 


M=Mo 
48 J 
Q, 10 ere ( ¢~T™ ) 


—2°5 log L+const. —o'84 —0°83 —068 —o40 -—o'21 —0'04 


3°245 3°283 3°408 4°006 


Q, 10“ erg ( 


3°235 3°22) 3S 33218 3 3205 3190-33173 


At Constant Mass and Potential Energy 
T., relative I*0000 §=60'9970 «09922: «009838 = 09804 «= 00'9779 = -00'9795 
p, relative 170000 ©60°9998_ =—-1'0235 —sdr'0980)=—s-d1'1631 = 1°2257 = 1°25 14 
R relative 1°000—s: 1003 1‘O15 1'053 I'10g9-—s«a1'238~—Ssd 609 
Transition I>2 2>3 3>4 4-5 5-6 6>7 
w/Q, 10> units +35 +257 +694 +525 +471 +153 


In Table III, the seventh line gives the potential energy on the assumption 
that the average release of nuclear energy per unit of mass and time is proportional 
to the 20th power of the central temperature, with the relative luminosities as in 
the eighth line (stellar magnitudes); this case, which may be near to the actual 
conditions for dwarf stars, yields an almost constant potential energy for the 
different composite models. 


“3. Opik, Tartu Obs. Pub., 30, No. 4, 1938. 
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5. Decay of the Disturbance through Homologous Transitions.—For a first 
approximation we neglect the pure changes of structure and changes in the 
boundary condition which actually should accompany every disturbance; the 
type of structure and the boundary condition we regard as rigidly imposed on 
the model during all the changes it is likely to undergo; similarly, we assume 
the equilibrium luminosity Ly = E, to be independent of the radius, disregarding 
thus factors of the order of R-®® which would follow from Kramers’ law of opacity. 
We remain thus within the frame of Section 3 of the present paper, assuming 
the response of the rate of expansion and of luminosity to a disturbance in the 
energy generation to take place instantaneously according to the simple formulae 
(5’) and (6’). 

We assume a law of energy generation E =C(1 +A)p,7% and assume that the 
disturbance consists in A changing suddenly from zero to a certain finite value 
and remaining constant after that. The equilibrium luminosity is thus Ey = Cpy)T%, 
Po To, Ro and Q, being central density, central temperature, radius and potential 
energy of the initial undisturbed model. By homology 7,~R™", p,~R® and it 
is easy to verify that after the disturbance has started 


E=E,(1 +2) (3) (20) 


The output of gravitational energy is a fraction 


y-5 4 


won 
of the variation of the potential energy when y= 3; hence 


H=35 =—2(E-E,) 
by (5°); ss 
Ry 
- 2? 


we obtain the differential equation (t = time) 


Ry dR 4E, R,\**3 
R? dt “= ak +) (7) -r]. (21) 


The expansion (A>0) or contraction (A<o) proceeds according to this formula, 
approaching asymptotically the new equilibrium state E, = E, which, according 
to (20), corresponds to an equilibrium radius 


R,=R,(1 +A)**s. (22) 


Q=0 


For small changes (21) may be simplified, by setting x=R/R,—1 and by 
retaining only members of the first order with respect to x and A; this gives 


dx 4E, 
= q, A-& +3) ], 


or, setting t=0 at x =0, we get 
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Qo 
t,= ———_—— 2 
° 4E(s +3) (24) 


represents the time of relaxation of the disturbance. Assuming a solar luminosity 
E, =3°8 x 10* erg/sec. and a potential energy of (2,=4 x 108 (Table III, Model 
No. 6), with s=15 we get tj = 4-4 x 10° years as a probable figure for the Sun. 
The disturbance of the luminosity decays with the same period, in agreement 
with (6) and (20), 
L—E,= —AE ge". (25) 


6. Complex Decay of Disturbance.—The above assumption of an instantaneous 
response of the luminosity to a disturbance, with homologous changes, is rather 
schematical. The expanding model does not undergo exactly homologous 
changes (cf. Table I), wherefore a certain energy of transition appears as a 
delaying factor; also, the boundary condition will be violated to a certain extent. 
Here an attempt is made at a second approximation in which the changes of 
structure and of the photospheric temperature are taken into account, although 
still rather schematically. As the initial type of static model we assume No. 7; 
according to Section 3, expansion requires an increase of M,/M, the relative 
mass of the core; from Table III we infer that this would correspond to a shift 
of structure from No. 7 toward the earlier numbers; we therefore assume 
conventionally that the structure may undergo a continuous displacement along 
the composite models of ‘Table III. As a continuous coordinate of structure we 
assume the relative energy of transition y=w/Q, starting from No. 7; this is 
found as follows from the last line of Table III: 


TasBLe [V 
Model No. 7 6 5 4 3 2 I 


y X 10° ° —153 —624 —1149 — 1843 —2100 —2135 
T(pp=const.) 5800 4880 4650 4557 4465 4396 4350 


The last line gives T,, the photospheric temperature which follows from the 
numerical solution of the models, assuming a constant value of the photospheric 
density, p,, chosen so as to give T,=5800 deg. for Model No.7. The external 
radiation is prescribed by 7, (and radius). Instead of the actual figures, the 
luminosity as function of the structural parameter y is represented by an empirical 
interpolation formula 


L=E{t—(—y + 10-19}, (26) 
satisfactorily corresponding to T, of Table IV; Ey is the equilibrium luminosity 
of Model No. 7 (y=0). 

The equilibrium luminosity £5 of different structures (yo) differs, however, 
from E,. The actual deviations log Ej/E, are identical with the figures of the 
eighth line of Table III. According to formula (3), E,;=Q is identical with a 
certain average value of the radiative flux inside the radiative equilibrium zone 
for the given set of models. This is well represented by the empirical formula 


E; =Q=E,[1 +197(—y)*%}. (27) 


The actual luminosity of the expanding model differs from Q and is given by (26). 
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The radiative flux has to take care of the following components: the energy 
from nuclear sources, E (positive); the gravitational energy, H, (negative for 
expansion), of the homologous component in the change of structure; and 
the flow of the energy of transition, the maximum value of which is given by 

dw dy 
J=— =Q= 28 
dt dt’ (28) 
a negative quantity for the transitions from No. 7 toward the earlier numbers. 
The radiative flux at the first interface will be 


O,=E+0:25H +J, 


with respect to E and H in agreement with Section 3 (H, =0:25H), whereas for J 
the coefficient should really be somewhat less than 1 (Section 4). 
The flux at the second interface will be 


0,=L=E+H*; 
according to formula (3), k = %, we obtain 

O=E+}H+ J. (29) 
This expression, being comparable with (27), is of fundamental importance 


for the solution of our problem. 
In agreement with formulae (5’) and (6’) we assume 
H 


E=E,-—, (a) 


Laka. 
2 


Ditferentiation of (a) yields 
aE __ dH m= 
:"  's ‘ 


On the other hand, assuming for the homologous component of the change of 
structure the formulae of Section 5 for small changes, (6) with (25) yields 


H =2(L — E,) = — 2AE,e**, 
whence by differentiation 
dH H 
Substituting (d) into (c), we get the last of our differential equations required for 
the solution 
dE _ cel (30) 
7 > at, ’ 


The above system of equations is conveniently reduced to relative variables: 
tT=t/ty; y=w/Q as before; z=E/E,; h=H/E,; j=J/Ey. Setting 
_Lotp si 


Q  4(s +3)’ 


ractically correct for transitions between Nos. 7 and 6 which alone come into question here, 
as shown below. 





*p 
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we obtain: 
=aj; 
dr y 


dz 
Pat al 


=+4h+%j=1+197(—y)*; 


10 


| 
| 
s+h=Ylr—(—y +10], | 


Table V represents the result of a numerical integration, for s=15, a=4; 


for r=0 the relative energy generation is assumed to increase suddenly from 
% =I to z=1°5(A=0'5). 

The first column gives the time in years, assuming for a solar mass the time of 
homologous relaxation equal to half a million years; the third column gives 
L/E , the relative luminosity; the fourth column gives 2—2z=2-—E/E, which 
would be the march of luminosity in the case of instantaneous response, according 
to formula (6’); the last column gives the effective radiative flux inside, which 
is seen to change very little. Remarkable is the great stability in the type of 
structure as revealed by y; the maximum change (at t=1-5 x 10°) corresponds 
to only about one-sixth of the interval No. 7~No. 6. On the other hand, the 
variation of the luminosity is very much influenced by the ‘inertia of structure 
and boundary condition’’, as follows from a comparison of columns 3 and 4 and 
as is shown in Fig. 1. 


TABLE V 
Decay of Disturbance by Simultaneous Expansion and Change of Structure 





t 
years y L/E, 
(to=5 X 10°) 


© 
® 





<o 1*000 1000 ° ° 
° 1000 0°500 0*500 
0°005 - ‘2 xX1071! 0995 *500 — *§05 370 
0°05 ‘8 x 1071 ‘979 "500 "521 "359 
0°5 ‘8 x107~° 947 *500 ess *335 
5 *§ X107* "907 *500 — °§93 *305 
50 ‘70 X 1077 ‘857 “500 643 
500 *s x107~° *796 *500 — "704 
5000 7.x 197° "722 "504 "774 +163 
5 x10! - x 107% 643 *540 — 817 073 
s X 10° x r9o~* 627 581 "792 — +030 
1*5 X 10° x 1075 622 “620 758 + *003 
2x 10° x z07* 627 657 “716 + +027 
2°5 X 105 ° ‘s. Xs0~° 638 692 670 + +044 
3 x 10° ‘9 x1075 658 "724 618 + +051 
+ 
+ 


SSS8888888 


88888 


3°5 x 108 oe x10" “690 "754 — +556 -048 
x 10° ‘2 x1075 0°753 0°780 —0'467 0°020 
4 


mt 





The complex decay (A) yields a very broad minimum, the lowest luminosity 
of 0-622 being attained 1-5 x 10° years after the start; the first approximation 
(B=2-—2) starts much deeper, at 0-5, then rises steadily. For t>4 x10 both 
curves approach the curve of exponential decay and merge into each other. 
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As to the response of luminosity at the beginning of the disturbance in case A, 
although not instantaneous it is nevertheless catastrophically rapid, noticeable 
already during the first year and exceeding one-half of the largest effect after 
500 years. 
L/E, 
1-0 





i 1 1 1 


l 2 3 4 t 
. 1.—Variation of luminosity caused by a 1°5-fold sudden increase in the energy output. 
Abscissae : time in 10° years, for a solar mass. 
Ordinates : relative luminosity. 
A: complex decay (second approximation). 
B: instantaneous response ( first approximation). 
Broken line: original (equilibrium) luminosity. 





7. The Ice Ages.—The problem can only be very briefly reviewed here. It 
is gratifying to notice that the standpoint assumed earlier by the writer with 
respect to paleoclimatic changes* coincides essentially with the views expressed 
by many leading modern geologists + : namely, that the normal average temperature 
of the Earth is higher than it is at present; that we are still inside a period of 
general deterioration of the terrestrial climate, a period already lasting for about 
a million years and connected with the phenomenon of Quaternary glaciation ; 
and that the theories of continental drift and of the migration of the poles cannot 
be made responsible for this, as well as for earlier disturbances of the terrestrial 
climate. The writer pointed out that a general increase of solar radiation, due 
to the decreasing hydrogen content of the Sun, would be responsible for a gradual 
‘secular’ increase of the mean temperature of the Earth which apparently may 
be traced between the Cambrian and the Tertiary. This normal trend is inter- 
rupted by more or less short periods of glaciation, spaced about as follows } :— 


rea ‘ , i Silur- Carbo- 
Glaciation Algonkian Cambrian Ordovician § Staten’ Deane Quaternary 


Time, 10° years 
_ 6 —4°8 nme ~ —2: 
from our era 7 4 3°8 ia 2°3 ° 





* E. Opik, Tartu Obs. Pub., 30, No. 3, 1938. 
t J. H. F. Umbgrove, The Pulse of The Earth, The Hague, 1947 (2nd edition). 


} The list of the main glaciations is the same as quoted by the writer in 1938, but the epochs are 
assumed according to Umbgrove. 


§ “ Mountain glaciation ” according to J. H. F. Umbgrove, loc. cit. 
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It is proposed to account for these climatic anomalies by certain repeated 
disturbances in the energy balance of the Sun as outlined in Sections 1 and 2 
of the present paper. The theory of a single disturbance is expounded in Section 6. 
The phenomena of the Quaternary glaciation are, of course, rather complicated. 
The smooth curve of solar radiation of Fig. 1 cannot explain the repeated advance 
and retreat of the glaciers known to have happened on a large scale at least four 
times during the Quaternary. However, even the warmer interglacial periods 
(to which our era belongs, too) are definitely colder than the preceding Tertiary 
or Cretaceous; for the latter Urey * finds, from the relative abundances of oxygen 
isotopes in deposited lime, a temperature at formation from +17°5 deg. to 
+26°5 deg. C., the average for the six studied specimens being +21-6 deg. C. in 
Hampshire, England. Even assuming this to refer to the summer temperature, 
it differs significantly from the present summer average of the sea temperature 
about +15-5 deg. C. near the surface. The fluctuations of Quaternary glaciation 
are thus secondary waves superposed upon a major general depression. For the 
major depression the above picture of the complex decay of a disturbance in the 
solar interior may be advocated. ‘The fluctuations may be produced by compli- 
cations from the outer shell, the inner disturbance affecting the extent of the 
peripheral convective region.t Repeated, or continuous disturbances (long- 
enduring propagation of the disturbance in the interior) may explain the slight 
deterioration of climate which is characteristic of the end of the Tertiary 
preceding the start of the first ‘acute’ glaciation. 

Terrestrial causes seem to be generally inadequate to produce overall changes 
of climate to sucha scale. The apparent (but loose) connection of past glaciations 
with the periods of mountain building (Holmes) may well be accidental, in view 
of the small number of coincidences; especially when we consider th4t mountain 
glaciation is possible only when mountains are present. ‘The amount of carbon 
dioxide in the air as a climatic factor cannot have any influence because absorption 
by water vapour covers more effectively all the regions of the spectrum where 
carbon dioxide may exert its ‘“‘greenhouse”’ effect. Moreover, vegetation is so 
efficient in removing CO, from the air} that, since the time when plants came 
into existence, the CO, content in the terrestrial atmosphere could hardly ever 
have exceeded the present amount, which actually approaches the minimum at 
which plants can exist (on Mars, according to Kuiper, the pressure of CO, is 
about the sameas on the Earth, a very significant fact suggestive again of vegetation). 
A different distribution of the areas covered by land and sea may be more efficient 
in producing climatic changes all over the Earth’s surface. The effect of an increased 
water surface would be a general lowering of the temperature, because of greater 
cloudiness, a smoothing out of seasonal changes and a smaller difference between 
pole and equator. However, the two hemispheres of the Earth, representing a 
difference in the land and sea distribution which is certainly larger than the 
possible fluctuations of the mean conditions on the Earth in the geological past 
have been, set an upper limit to the expected effect: the southern hemisphere 
is indeed cooler, but only by about two degrees (Hann); seasonal differences 
of temperature are largely smoothed out; but the maritime softening of the 

*H. Urey, Science, 108, 489, 1948. 
t+ F. Hoyle has recently advocated an external cause for such fluctuations which perhaps is not 
quite so open to objections which are adduced at the end of this paper against meteoric effects. 


} The present amount can be removed by the photosynthesis of plants within but a few hundred 
years. 
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climatic difference in latitude is hardly perceptible. There is hardly more 
of an ice age in the southern hemisphere than in the northern. From this it 
appears that systematic transgressions or regressions of the sea on a world-wide 
scale are still inadequate to account for the large climatic changes of the past. 
Besides, the warm period of the Tertiary differs from the present not only by an 
elevated average temperature, but also by a smaller climatic difference, the warming 
up being the most pronounced in the polar regions; this is hardly compatible 
with a more continental climate, whereas an increase in solar radiation alone 
explains it: at the elevated average temperature the moisture content, thus also 
the heat content of the air, is considerably increased, and the atmospheric circu- 
lation is more efficient in the convective transport of heat; as a result the 
temperature differences between the climatic zones are reduced. 

It seems difficult to escape the conclusion that the main cause of these climatic 
changes is to be sought outside the terrestrial globe. 

There have been several attempts to advocate changes in the eccentricity 
of the Earth’s orbit and in the obliquity of the ecliptic as a cause of the ice ages; 
since the older work of Adhemar (1842) and James Croll (1864-1889), this idea 
has been considered repeatedly ; the most detailed treatment is by Milankovitch * ; 
a calculation of the perturbations of the Earth’s orbit and obliquity has been made 
for a similar purpose by Pilgrim ¢ for the past 800,000 years. All these considera- 
tions have the same weak point in common, namely, in ascribing to small causes 
large consequences out of any proportion. The obliquity of the ecliptic influences 
the annual distribution of heat received at different latitudes, greater obliquity 
corresponding to increased insolation at latitudes above 45°; the annual effect 
according to Milankovitch amounts to +0-78 per cent at latitude 60° per degree 
of obliquity and, as the extreme range of obliquity is from 22°-2 to 24°-5 
according to Pilgrim, a maximum range of insolation of +0-8 per cent results 
from this cause; the range in the black-body radiative temperature corresponding 
to this is +0-6deg.C.; the average temperature of the Earth is not influenced 
by this cause, a rise in high latitudes being counterbalanced by a decline within the 
tropics and vice versa. Nothing of the size of the actual climatic changes can 
be explained by such a minute effect. The eccentricity of the Earth’s orbit 
affects the intensity of insolation in different seasons, by small amounts about 
as the above, ‘but the total quantity of heat received during the year remains 
unaffected by it; only a second-order effect can be expected from that, largely 
smoothed out by the Earth’s predominantly maritime climate: the result is 
doubtless nil. 

Even if we admit with Milankovitch and others that the advancing and 
retreating pulse of Quaternary glaciation reflects the combined fluctuations in 
the eccentricity { (e) and obliquity (e) of the Earth’s orbit, the absence of 
glaciation in the Tertiary and earlier remains still unexplained; the fact that 
Pilgrim did not dare to push his calculations back farther than 800,000 years 
before our era does not mean that the fluctuations did not exist before that; 
the variations of e and e depend upon planetary perturbations and are continuing 
certainly all the time the Earth and the solar system exists. Obviously something 


*M. Milankovitch, Théorie mathématique des phénoménes thermiques, etc., Paris, 1920; Kalorische 
Fahreszeiten und deren Anwendung im paldoklimalen Problem, Beograd Akad., 1923. 

+t L. Pilgrim, Versuch einer rechnerischen Behandlung der Eiszeiten, Stuttgart, 1904. 

} Combined with the precession of the equinoxes and the motion of the perihelion. 
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of a more fundamental change has occurred at the end of the Tertiary, and unless 
an explanation is given for the mere start of the first glaciation it may seem 
superfluous to argue about its subsequent course. 

The most elaborate, and probably the most futile, attempt to devise a theory 
of the ancient climatic changes on the Earth has been made by W. Képpen and 
A. Wegener.* Of the 256 pages of this treatise 99 deal with the Quaternary. The 
authors accept Milankovitch’s and Pilgrim’s astronomical data for a basis ; however, 
they are confronted with the insufficiency of the astronomical cause. They do 
not realize that the true fluctuations in temperature will always be smaller than 
those calculated from a black-body radiative equilibrium formula, because for 
a constant annual amount of heat received by the whole Earth local or seasonal 
anomalies will be counterbalanced by the store of heat in the oceans and by 
atmospheric circulation having an equalizing effect; the theoretical amplitude 
of +0-6 deg. C. should be reduced to from one-half to one-third of that amount. 
Képpen and Wegener take an opposite course, applying to the small figures a 
certain numerical amplification. Moreover, they restrict themselves to the 
summer solstice when the fluctuations are the largest, and instead of using 
directly the data for insolation, they convert them into effective shifts of latitudes. 
Because the curve of insolation runs rather flat with latitude at the solstice t, 
small changes of insolation are represented as big shifts in latitude, and the 
shifts in latitude so obtained are converted into temperature changes by using 
mean temperatures for the corresponding latitudes. Such a method yields final 
changes which are 10-20 times larger than the actual variation in the radiative 
equilibrium temperature, thus up to 50 times exceeding the probable variations 
expected. But this amplification is still insufficient, and Wegener’s theory of 
continental drift is introduced: the northern continents have to be displaced by 
about 1000 km. in half a million years. This is not enough, and so the pole of 
rotation of the Earth is to be moved by 20 degrees of latitude during 600,000 
years. And now one might expect that after having adjusted arbitrarily so many 
factors everything fits well into the scheme, but this is not so. Numerous facts 
which do not fit are either dismissed or ad hoc explanations are produced. As to 
the displacement of the pole of rotation, the apparent changes in the barycentre 
hitherto observed during half a century are entirely explained by errors in the 
stellar proper motions used and by changes in the observing programme. f 
But assuming the changes as they stand, they imply a motion of the pole of only 
0”-06 per century, or 600” in one million years, 200 times less than assumed by 
Képpen and Wegener. Also all hitherto announced directly observed continental 
shifts are traced to systematic errors.§ Leading paleontologists strongly oppose 
the conception of continental drift as proposed by Wegener and it has been stressed 
‘that the continents were essentially stable throughout the whole time involved 
in mammalian history (i.e. since the Triassic)”’.|| All the elaborate edifice of 
Képpen and Wegener breaks down, as not one of their arguments has a real 
foundation. Since the beginning of the Tertiary at least, and very likely for the 


* W. Képpen and A. Wegener, Die Klimate der geologischen Vorzeit, Berlin, 1924. 

+ The curve of geometrical insolation even rises toward the pole; however, Milankovitch 
introduces a correction for atmospheric absorption (rather primitively) and obtains a slow decline 
toward the pole. 

t T. Nicolini, Draft Reports of I.A.U., Ziirich, 1948, Commission 19, p. 36. 

§ Ibid., p. 39, and Wm. Markowitz, Trans. Amer. Geophys. Union, 26, 197, 1945. 

|| G. G. Simpson, Amer. F. Sct., 241, 1, 1943; also Bull. Geol. Soc. Amer., 58, 613, 1947. 
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earlier periods, we may assume that the climatic changes recorded in the terrestrial 
rocks reflect chiefly the changes in the intensity of solar radiation, only partially 
modified by local factors. 

The radiative equilibrium of the Earth may be influenced by other cosmic 
factors also. For example, meteoric dust has been many times advocated as a 
possible cause of the ice ages. There is a fundamental misunderstanding about 
this point. Meteoric particles by their impact will produce extra heat on the 
surfaces of the Sun and the Earth, but any absorption of solar radiation will be 
more than compensated by the re-radiation of solar heat from the dust cloud 
inside and outside the Earth’s orbit. A thorough discussion by V. Riives* settles 
definitely this point in the sense that from such a cloud only a warming effect 
can be expected. If the problem is inverted and it is assumed that the “ normal” 
state corresponds to a meteoric medium sufficiently dense to increase the solar 
radiation by 10-20 per cent, an impossibly high space-density of galactic matter 
must be postulated. 

It seems therefore that disturbances in the energy balance and internal 
structure of the Sun, perhaps complicated by secondary disturbances in the 
hydrogen (or hydrogen+helium) content of the subphotospheric convective 
region and leading to temporary disturbances in solar radiation, offer the most 
plausible explanation for the phenomenon of the recurrent ice ages on the Earth. 


Armagh Observatory, 
Northern Ireland: 
1949 March 1. 


ADDENDUM 


The behaviour of the luminosity of a non-static model, considered in 
Section 3, may be commented upon also in the following way. Ina static model 
with a core-source of energy the radiative flux between r, and 7, (Fig. 2) is constant, 
being a definite function of the parameters and boundary conditions involved: 


O=/(M, To Pe T ., Pp» Bs Kqy «+ -) = Ep. (1) 
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This function cannot, of course, be expressed explicitly, but it is a definite 
outcome in the process of numerical integration and may be called the “ radiative” 
equilibrium function”. In Fig. 2 the straight line AB represents the behaviour 
of 0, =Q=Ar,=const. over the radiative equilibrium zone. 

* V. Riives, Tartu Obs. Pub., 30, No. 7, 1940. 
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Let us now abandon the condition of the constancy of the radiative flux. 
This may be done by interposing arbitrary positive and negative energy sources, 
and Q, may be represented by the dotted line in Fig. 2. Under unchanged 
conditions (parameters, etc.), and when the deviations are infinitesimal, the average 
value of the radiative flux will equal the static value, O0,=Q=E,y. Now, the mean 
value theorem requires that Q,,.; >>O,>Qmin OF QOmax >Eyp>Qmin- This is true 
also for finite small changes and, as calculations indicate (Table I), remains valid 
even for large finite changes. 

Let us now consider the case of an expanding model. Negative gravitational 
energy sources are present and the radiative flux is steadily decreasing outwards 
(9,>Qs2, Fig. 3); QOmax=Q1» Omin =Q2; thus, according to the above, 


0,>E,>Q2. (2) 


This is the principal inequality used. 


Q 4 








Ly 
R 





Fic. 3. 


In an actual model the parameters of (1) are changing with the change of 
structure, and therefore the value of Q in (1) will somewhat change also and 


instead of (2), 
0,>90>9, (2') 


should be written. Actual integrations, however, indicate that the deviation of 
QO from E, is so small that (2) remains always valid. 

Let us consider the transition from a static case 0, = E)=const. toward the 
non-static case 0,0,Q, to be obtained by the cumulative effect of single pertur- 
bations 80 = QO, — Ey, each covering the range dr (Fig. 3); a single perturbation 
5Qdr will, in the process of outward integration, produce in the boundary 
condition a variation F(r).5Qdr. As the radiative equilibrium solution in 
outward integration is divergent, the character of the divergence being an 
exponential one*, and as the region between 7, and R is irrelevant as governed 
by a fixed adiabatic equation of state, F(r) will represent a rapidly increasing 
function of r,—r, such that F(r,): F(r,) may be of the order of 10“ and more. 
The total perturbation in the boundary condition will be 


5B = [ "89. Fin dr. (3) 





* E. Opik, cf. Tartu Obs. Pub., 30, No. 3, 1938. 
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Setting 5B =0, which is justified because the physical boundary condition (7, p,) 
is extremely stable, being limited by the laws of photospheric opacity and radiation*, 
we have acondition 


f "30. F(r)dr=o0 (4) 
which, by the theorem of mean value, implies at once 


QO, >Ey>Q2. 
In the practical application EZ, has been assumed a weighted mean of the extreme 
values, 


E,=kQ, +(1—)Q2, 


which is formula (3) of the paper. The actual shape of the function Q, as the 
sum of a constant core-energy and the negative energy of expansion proportional 


7M 
to | : T dM, is about as pictured in Fig. 3, convex toward the r-axis, and the 


0 
weight F(r) of the perturbation 5Q in (4) increases strongly toward r,; under 


such circumstances the fraction k exceeds } in any case, and for the core-source 
models considered here it can be shown that k fluctuates within a narrow range 
around %. The actual integrations for models of uniform molecular weight and 
of moderate p,/p,, lead to k= % very closely. 


* This differs from a mathematical boundary condition which may be varied at will. 
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Summary 


The work of a previous paper (1) by the author on the effect of steady large 
scale convection currents on a magnetic field is extended to include turbulent 
convection. It is shown that turbulence reduces the effective electrical 
conductivity of a gas, in some cases to a small fraction of its normal value. 





1. Introduction.—It has already been suggested by Cowling (2) that the decay 
of the Sun’s general magnetic field might be accelerated by the turbulence in the 
convective core breaking up the field into small pieces which could decay inde- 
pendently. This, however, cannot be the actual mechanism of the decay, for 
however large a number of pieces were broken from a given set of tubes of force 
in a given region, the tubes of force themselves must still thread the region as a 
whole. They would in general merely be displaced from their original positions 
and compressed, as the volume they occupied was diminished by the volume 
occupied by the broken pieces. In general a field can only decay as the complete 
closed tubes of force constituting it disappear by contracting to a point. It will 
be shown that the effect of turbulence is to facilitate this contraction. 

A representation of the field in a conducting medium in motion is established 
in Section 2, which makes it possible to see qualitatively why turbulence should 
reduce the effective conductivity. The mathematical treatment is then given 
in Section 3. 

2. A physical representation of a magnetic field in a conducting meditum.—Alfvén 
(3) has given an excellent picture of the lines of force of a magnetic field in a 
perfectly conducting medium, as elastic strings embedded in the medium. This 
picture can be extended to allow for the finiteness of the conductivity. In this 
case it is possible for material to flow across the field, thus introducing a definite 
velocity of the medium relative to the lines of force. It is necessary therefore to 
make a clear definition of the velocity of the lines of force. 

Let K and B be a solenoidal vector field and its vector potential, respectively ; 
the stream-lines of K are merely the mathematical family of lines given by 
Lagrange’s differential equation 


K. ss & (2.1) 


z vy 


Any particular line can be labelled at a given instant by stating the point in which 
it intersects some reference surface S, intersected once by all the stream-lines in 
the region considered. The physical conception of a velocity can be attached 
to a stream-line only by assigning a velocity to its point of intersection with Sp. 
No physical significance can be attached to velocity along the line; the velocity, 


REY RN a ae, 
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denoted by u, will therefore be taken as perpendicular to the line. It can be 
determined in terms of K and of its values on So, as follows :— 

The ordinary theory of differential equations (4) shows that, in general, two 
independent families of surfaces can be found to satisfy (2.1). Let ¢(xyzt) =const. 
and ys(xyzt) =const. be two such independent families of surfaces containing the 
stream-lines. Now K can be expressed in the form 


K=F grad ¢dagrady, (2.2) 
where F is some function of position and time. Since 


div K=o, 


grad F .(graddagradw) vanishes, and F is therefore a function of ¢ and #. 
The vector potential can therefore be written 


B= i) ” Fdbgrady +grad W, (2.3) 


where W is an arbitrary function. 

Consider now a stream-line defined by 6=A, 4 =p, where A, are constants. 
Adopting Lagrange’s treatment, select any point on a stream-line (A, 1) at time 
t =o, and let (x,y,z) be the position of this point at a subsequent time ¢t. Hence 
x=x(A, p, v,t), etc., where v is a parameter determining the distance of the original 
point along the line (A,). We can therefore write 


0¢/dt +u . grad 6=dA/dt =o, (2.4) 
where u=(0x/0t, Oy/Ot, 0z/dt) is the velocity of the stream-line. Similarly 
0/Ot +u. grad =dyu/dt =o. (2.5) 
Thus 
ua K= Faa(grad¢a grad) 
= F(u. grad) grad ¢ — F(u. grad ¢) grad & 
= F0¢/dt grad p — Foys/dt grad ¢, (2.6) 


on substituting from (2.4) and (2.5). Finally, on taking the curl of both sides, 
while noting that F= F(¢, #), this equation reduces to 


curl (uaK) =0K/2t, (2.7) 


which is the required relationship determining uw in terms of K and the boundary 
conditions on Sp». In order to show that the boundary conditions on Sy are 
sufficient to determine a unique solution, suppose it were possible for two distinct 
solutions u,, U, to exist, taking the same values‘on Sy. The vector w=u,— a, 
would vanish on So, and would satisfy the equation 


curl (wa K)=0 (2.8) 


everywhere. On noting that w is perpendicular to K, it will be shown by the 
following lemma, that w vanishes everywhere, thus proving the uniqueness of 
the solution. 

Lemma.—lIf w is a vector satisfying curl (w,~ K)=o in a given region, and is 
such that wa K vanishes on a surface S, intersected by all the stream-lines of K 
passing through the region, then w~ K vanishes throughout the whole region. 

Since curl(w~ K)=0, there must be a potential function W such that 


w»aK= —grad W. 
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Both w and K lie in the surfaces W=const., hence W=const. on a stream-line 
of K. But gradW=o0 on Sy, hence W=const. throughout the whole region 
considered, which proves the lemma. 

The boundary values of u over S, cannot, however, be assigned arbitrarily, 
as the variations of u over S, must always satisfy (2.7), or strictly, the component 
of (2.7) normal to So, since it is only this component which involves exclusively 
the derivatives of ua in the surface Sy. Now the general solution of (2.7) can be 
written 


ua K=ua,a K-—grad W, (2.9) 


where dy is a particular solution, and W is a function which can be assigned 
arbitrarily on Sy. W is then uniquely determined everywhere because, in virtue 
of grad W being perpendicular to K, Wis constant on the stream-lines. On taking 
the vector product of both sides of (2.9) with K, and noting again that w has been 
defined as being perpendicular to K, the explicit general solution is thereby seen 
to be 

u=u,+grad Wa K/K?, (2.10) 


where K is the magnitude of K. This also shows that the boundary conditions 
for u can always be expressed in terms of a particular allowable set of values 
u, on Sp, together with an arbitrary function W, given on Sp. , 

It can be shown, by a generalization of Kelvin’s theorem concerning the 
circulation of fluid round any closed curve moving with the fluid, that 


curl(va K) =0K/at (2.11) 


is the necessary and sufficient condition for the flux of K through any closed curve 
moving according to a velocity field v, to be constant in time. Equation (2.7) 
shows that, in particular, the velocity ua of the stream-lines of K satisfies this 
condition. The flux of K through any curve moving with the stream-lines of 
K is therefore constant. This justifies the physical picture of tubes of magnetic 
induction which can only vanish by contracting to a point, or by converging on some 
limiting line and disappearing into it in the manner of a line sink, and vice versa. 
In a decaying magnetic field, for example, the tubes of magnetic force can only 
disappear by contracting to a point on a null line of the field, while if, on the other 
hand, the field is being built up by, say, a current flowing down the length of a 
cylindrical shell, tubes of magnetic force will appear from the shell, already 
formed. 

By taking K to be a magnetic field H, (2.11) is the same as the electromagnetic 
equation for a perfectly conducting material of unit permeability, moving with 
velocity v and carrying no externally impressed E.M.F.’s, provided v does not 
vary rapidly enough to produce an appreciable electromagnetic radiation. The 
above-mentioned extension of Kelvin’s theorem, applied to magnetic fields, 
therefore gives the well-known result that, for a perfectly conducting medium in 
motion, the flux of magnetic force through any closed curve moving with the 
medium is constant in time. This property gives rise to the usual picture of a 
field frozen into the material carrying it. 

This picture can be developed to include the effect of the finiteness of the 
conductivity. Consider further the equation curl(v, H)=0H/dt for a perfect 
conductor. Equation (2.7), which holds for any solenoidal vector, shows that if 
u is the velocity of the lines of magnetic force, then curl(a, H)=0H/dt, since 
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His solenoidal. (It is assumed here and throughout this paper that the material 
carrying the field has unit permeability.) Thus if v-u=w, the velocity of the 
material relative to the lines of force, then curl(waH)=0. The lemma then shows 
that, by defining u as the component of the material velocity normal to H on the 
reference surface Sy, wu is equal to this component throughout the whole region. 
The lines of force therefore move with the material, which completely justifies 
the picture of a frozen field. 
In the case of finite conductivity, the electromagnetic equation is 


curl H/470 =va H +cE™ —c grad V —0A/dt, (2.12) 
where o is the electrical conductivity expressed in e.m.u., E** the impressed 


E.M.F., due for example to thermal effects, A the vector potential, and V 
an arbitrary function. Now (2.7) shows that 0A/dt can be written as 


0A/ot=urnH +grad U, (2.13) 


where U is an arbitrary function of position and time. On substituting for 
0A/dt, (2.12) then becomes 

curl H/47o0 = waH +cE™ —c grad V, (2.14) 
where U has been absorbed into the function V. 

In cases where curl H and E™* are perpendicular to H, such as, for example, 
the Sun’s general magnetic field or sunspot fields neglecting currents circulating 
in meridian planes, u can be defined so that waH=curl H/470 —cE™' on the 
reference surface Sy. The work of the earlier part of the section then shows that 


curl H/470 =waH +cE™ (2.15) 
everywhere, while the component of w normal to the lines of force can be written 
explicitly as 

w,, = Ha(curl H/470 —cE™')/H?. (2.16) 
The first term of this expression admits of a semi-geometrical interpretation, as 
follows :— 
Put H= Ht, where t is the unit tangent to the line of force. Hence 


Hacurl H= H*tacurlt + Hta(grad Hat). 
By ordinary vector theory it can be shown that 
tacurl t=} grad t? —(t. grad)t = —(t. grad)t, 
since t?=1. But (t.grad)t= —«n, where « is the curvature of the line of force, 


and n the unit normal drawn outwards from the convex side of the line. We 
therefore have 


w,, =«n/4r0 +ta(grad Hat)/4n0H +cE™H/H?, (2.17) 


which gives a clearer interpretation for w,, than the expression in (2.16). 
Consider now the effect of the magnetic field on the motion of the material 
carrying it. ‘The equations of motion are 


Dv 
Pp, = —Stadp +pg +", 


where F™ is the mechanical force exerted by the magnetic field on the material. 
Now F”= — Hacurl H/4z, which can be substituted from (2.16) to give 


F" = —ocH*w, +ocE™,H. (2.18) 
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The equations of motion can therefore be put in the form 


Pr = —grad p +pg —oH*w,, +ocE™, H. (2.19) 


In the absence of any external E.M.F. E**‘, the magnetic field can therefore be 
regarded as a porous elastic material allowing free movement parallel to the field, 
but resistant to motion across it, with a coefficient of resistance oH?/cm.3/unit 
relative velocity. 

With regard to the decay of a field in a stationary medium, (2.17) shows that 
the velocity of the lines of force contracting through the medium decreases as the 
conductivity of the medium increases. This is another aspect of the well-known 
property that the rate of decay of a field decreases with increasing conductivity. 
The conductivity may therefore be regarded as a measure of viscosity in the medium 
in resisting motion of the linesofforcethroughit. Ifthe lines of force are distorted 
by turbulence of the medium, thereby increasing their curvature, (2.17) shows 
also that their general mobility may possibly be increased, with a consequent 
lowering of the effective conductivity. This cannot, of course, be taken as a 
proof that such a reduction actually occurs because, although the numerical 
value of the curvature at any point is increased by the turbulence, the effective 
contribution to the mean relative velocity w,, has different signs at different 
points. The theory in this section has been developed principally as an intro- 
duction to the theory of turbulent conductivity given in Section 3. 

3. The theory of turbulent conductivity.—Consider a field such as the Sun’s 
general magnetic field, symmetrical about an axis, and distorted by a meridian- 
plane circulation velocity. Using cylindrical polar coordinates as in a previous 
paper (1) by the author, the meridian-plane components of the field Hy which 
would be produced in the absence of turbulence are given by 


where B,=B,(za, 2, t) is a stream function constant on the lines of force. If the 
lines of force are represented by the family of curves uo(w, z, t) =A, then By = B,(A, t) 
and H)=G,/h,,m, where h, AA is the distance between adjacent curves uy=A, 
Uy =A+ AA, and G,=dB,/dA. 

The effect of turbulence is to distort the lines of force about their mean 
positions 7=const. Assuming that it does not affect the essence of the problem, 
the turbulence will be taken for simplicity to occur as tight folds in surfaces of 
revolution about the axis, thereby admitting a two-dimensional treatment. The 
distorted lines of force will be denoted by u(a, z, t)=const., so that any line 7=A 
is distorted into the line u=A. The distorted field H is given by 


_10B_ Gau 10B_ , Gau 


A. ane 


ado wow’ 
sothatH=G/h,w. Here G=0B/0A, where B = B(A, t) but is not the same function 
as By). The mean field H,, is therefore given by 

Gau 


a ~ waz’ 


G 
A, Hy.= + — ’ 


so that H,,=G/hzw. 


re 
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B is determined by the ¢-component of the electromagnetic equation in the 
form developed in Section 2, viz., 


Fag CU = wall +cE™—cgrad V;, (3.1) 


where E**' is as usual the impressed E.M.F., V is a function which is uniquely 
determined once the velocities of the lines of force are assigned, as in Section 2, 
and w again is the velocity of the medium relative to the lines of force. All the 
functions are independent of ¢, hence B is oe : 


@ (1 aB)  @B Est, (3 
- of = Bo +53 |4rom=—|w, = S +e [w+ (3.2) 


On multiplying both sides of this equation by pw?, where p is the density of the 
medium, and integrating over a finite narrow strip AS in the (w, z)-plane between 
two adjacent lines of force u=A, u=A+AA, it follows that 


-2f $4 = (= 0B + | 
dB dB ener 
iyo ae E*tdS. (3.3) 


The left-hand side of this equation and the first term on its right-hand side may 
be transformed into line integrals by Stokes’ theorem. After some reduction and 
rearrangement of terms, while noting that H is normal to, and w parallel to the 
short ends of the strip, this results in the equation 


{ca+anf ds — Go| as | 4 


2 if =| 5 -ag(=) + = 2 feeble 
ene ee Ow Om \ o + og ‘€z\ o 
= ~{B(A+Aa) - BQ)} | paresis 
u=A+i 
—|{ (BA)—B)wdiv pwdS +c | | pwE*'dS, (3.4) 
JJ as AS 
where w=|w|, and ds is an element of length of a distorted line of force. On 


dividing by AA and proceeding to the limit AA-o, while noting that dS =h,, Ad ds, 
this equation can be written 


1a(af If x1fdB a (pa) aB a (pa? 
xanlef_2 sal rds) + te sabe () + be-ag( oe) [ues 


=— G| pww ds +c | pwEh, ds. (3.5) 
u=A A 


0B 2 [pw?\ aB 2 
tos) + ae (3-0) 
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where 0/dn denotes differentiation along the normal to the curve u=A in the 
direction of increasing A. Moreover, since the strip is long compared with the 
dimensions of the turbulence, then 


Jie "S)** |e) * Gay" 


where 0/07 denotes differentiation along the normal to the mean curve # =A in the 
direction of increasing A. Hence 


" 1 [0B @ (pw*\  aB 2 (pa? ai pal ere 
| sel dma) +5-a5(SS) |beee--ef 25/( a 


‘ [___{Hanerad(®)} teas (3.8) 


Some of the remaining integrals can be expressed in terms of certain mean values 
characterizing the turbulence, but independent of the strength of the field. Put 


j,_,a/[,_,as=« and fads [[  heds=nir, (3.9) 


where the integrals are taken along sections which are long compared with the 
dimensions of the turbulence but short compared with the distance of variation 
of the mean field. Here « is the ratio of the length of a distorted line of force 
to the corresponding section of the corresponding mean line; thus y=A;/h,, 
where h,, is the mean value of h, in the region concerned. + is thus a measure of 
the constriction of the tubes of force consequent on the distortion; hence y>1I 
in a highly turbulent region. It should be noted also that h;/h, =H/H,,, hence 
y =H /qH,,, where H is the mean value of the scalar strength of the distorted 
field, while g= H[1/H], and is thus a mean function of the turbulence of order 
unity. 

In order to examine « further, consider a strip of area AS, between two mean 
curves i=Ay, 4=A,+AA, and small compared with the area of variation of the 
mean field while large compared with the area of a turbulent eddy. Now 


ee, wae 


Ag+ BA 
[| hgdsda=% { | h,,ds d= ¥ AS, (3.10) 
=A KJ 4, u=m2 K 


hence x =y. 
Consider now the mass in the region between the moving lines u=A, #=A. 
The rate of flow of mass out of the region across a strip of the surface of revolution 


formed by rotating the line u=A about the axis of symmetry is 27 i] yn gh twds, 
while the rate of flow into the region across the corresponding strip of the surface of 
revolution of the mean line 7=A is 27 f _pwwt ds, where @ is the normal component 
of the mean mass velocity of the material across this mean line. But the total 
effective area of the region vanishes, since areas cut off by u=A on opposite sides 


* By Stokes’ theorem we have 


SJslee{eaa(F)} + 5a(F)] s- Jose) 4-J..ca(F)« 


where S is the region between the curves i=A, “=A. The left-hand side vanishes since the mean 
area between the curves vanishes. (3.7) then follows. 
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of the mean line make contributions of opposite sign to the total rate of increase 
of mass in the region, hence 


pww ds = put ds. (3.11) 
/u=A u=A 
By taking the length of the mean strip to be small compared with the typical 
distance of variation of the mean field, but still large compared with the dimensions 
of the gpm (3.5) can now be written 


2 if i 
zalJ._,— : 7H, 4s) +i {Hn agrad (= = hig ds 


=| porwraHjphzds+c| pwtEShzds, (3.12) 
J gaa Jai=a 


where it has been noted that « =y and H=qyH,,. On making a retransformation 
by Stokes’ theorem, and noting that H,,, is normal to the curves 7=A, 7=A+AA 


bounding the mean a AS at its short ends, we can write 


X(\_.Sertet)- Sli, (Eorm)} a 


Equation (3.12) can then be written 


all, [ {eur (oH, )f qS+- 20 —|_ { H,,ngrad (= “)t h,, ds 


a [_ pee eer pwrEsth, ds. (3.14) 


“USA 


Finally, on dividing by A;AS, where AS is the length of the mean strip, while noting 
that dS =h;Ad ds, (3.14) can be written 


z 2 + pw — 2 Frext 
= {eu (= - qy7H,, )f + {Hm asrad ( : )} =pw*{waH,,}, +cpw ES. 


47 
(3-15) 
Now w=0, —@, where 6, is the component of the mean mass velocity 0 of the 
material normal to H,,, and @ is the velocity of the mean line of force. Hence on 
substituting {GsH,,},=0A,,,/0t, where A,,, is the ¢-component of the mean 
vector potential A,,, (3.15) becomes 


2 2 
a {eur (= H,,)} ieee {Hy agrad (=) 
47 oT 47 0 J), 


=pwOrH,}, +cpw*ES — pw0A,,,,/dt, (3.16) 

where 
oy =0/g/". (3-17) 
(3.16) is therefore the electromagnetic equation determining the mean field H,,, 


and is best considered by reducing it to a more normal form by multiplying by 
470;/pa*, The equation is now 


pos {oot( Ba)), + par rere ("S) 


= 4707{0AH,,}, +470,cES — 470,0A,,,/0t. 
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In highly turbulent regions y>1; consequently o;<o, and (3.18) can be 
written grin as the ¢-component of 


curl H,, +s ot “T, grad (= =) H,,=470,0 aH, + 470,cE™ — 470,0A,,/0t. (3.19) 


Except for Be second term on the left-hand side this is the same equation, in every 
respect, as for the field in a medium of conductivity o;. o; is therefore the 
reduced effective conductivity of the turbulent medium, at least with respect to 
the meridian-plane component of the field. 

The theory just given does not apply to the meridian-plane components of 
(3-19), but it is tempting to extend the principle of the reduction of conductivity 
to these components also. Such an extension would enormously increase the 
difficulty of the analysis, and it will not be attempted in this work. A statistical 
method might be more tractable, although it is questionable whether such an 
approach is fundamentally possible. 

The effect of the mean mass velocity DB of the meridian-plane circulation can 
be examined by the method given by the author in a previous paper (1). It was 
shown there that a mean circulation velocity in meridian planes merely tends to 
distort the field to follow the lines of flow, without seriously affecting the strength 
of the field. 

If, as just mentioned, (3.19) could be taken to hold with respect to its meridian- 
plane components, then the effect, on the field, of differential angular velocity 
could be examined as in the author’s previous paper (1). The lowering of the 
conductivity would allow a greater measure of differential angular velocity along 
the lines of force of the meridian-plane component. ‘The results of this previous 
paper show, however, that in the reversing layer the transverse component H, 
which would result from a variation in angular velocity of the observed order, 
along the the lines of force, would be of the order of 10® times the strength of the 
meridian-plane component. The effective conductivity would therefore need 
to be reduced by a similar factor before any appreciable angular velocity variation 
would be possible along the lines of force of the meridian-plane component. 
The actual extent of the reduction of conductivity due to turbulence is discussed 
in the next section. 

Illustrative example.—Since the results of this section represent a great 
departure from the classical concept of conductivity, an illustrative example will 
be taken in which the problem is reduced to its bare essentials. 

Consider the magnetic field in a long circular cylinder of uniformly conducting 
incompressible material, produced by a uniform E.M.F. E impressed along the 
length of the cylinder. If the material of the cylinder were not in motion the lines 
of force of the resulting field would be concentric circles about the axis of the 
cylinder, with a field strength H, at distance r from the axis given by 


Hy =2trocE. (3.20) 


The effect of turbulence is to distort the lines of force into crinkled curves following 
mean concentric circles. The turbulence considered takes the form of steady 
motion in small cells where, although the velocity varies rapidly from point to point 
over a right cross-section, it is constant in time at a given point and is constant 
along the length of the cylinder. The electromagnetic equation is therefore 


curl H =470vunH + 4r0cE, (3.21) 
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where v is the velocity of the gas and E is a vector of strength E parallel to the axis 
of the cylinder. Only the axial component of this equation is non-vanishing, 
the vector form being retained for clarity. If the integral of this equation is taken 
over the area S enclosed by a distorted line of force C we have 


|| culH.dS=4ro || vnH.dS+4noc|| E.dS. (3.22) 
Since His solenoidal we can write 
[{ vaH.dS= | | v.grad BdS, (3.23) 
~/S Ss 


where B is a stream function for H, constant on the lines of force. On expressing 
v. grad B =div(Bv) — B divv, transforming the div(Bv) term into a line integral, 
noting that B=const. (By say) on the bounding line of force, and subsequently 
retransforming into a surface integral, we have 


i v.grad BdS = | I, (B, — B) div ods. (3.24) 


But divu =o everywhere since the material is incompressible, hence 


| | eaH-dS =o. (3.25) 


After transforming its left-hand side by Stokes’ theorem, while noting that E is 
uniform, (3.22) can now be written 


| Hds=4nocE || dS. (3.26) 
Cc J/JS§ 


For a line of force whose mean circle radius is 7, (3.26) becomes 


2nrcH = 4n0cEnr®, (3.27) 


where, as in the preceding theory, « is the ratio of the length of a distorted line 
of force to the corresponding section of the mean line, while H is the mean value of 
H along a section of the line which is long compared with the distance of variation 
of the velocity. The mean field H,, can then be written 


H,, =20,cEr, (3.28) 


where o7 =0q(H,,/H)*, while g= H[1/H] as before. 

It will be noticed that o; is not a constant of the material but depends also on 
the strength of the mean field. Equation (3.28) provides a simple example of 
how the whole character of the electromagnetic equation is changed by this 
circumstance. On substituting from (3.20), (3.28) can be written 


qH,,H, = H?, (3.29) 


where H, is the field which would be produced in the absence of turbulence. 
Now since H,, <H, (3.29) shows that H,,<H<H,, so that in cases where the 
source of the field is in a turbulent region even the mean scalar distorted field 
strength H is reduced to below the normal field strength Hy. Equation (3.29) 
alone is not sufficient to determine H,,. However, in anticipation of a result 
obtained in the next section we can quote an additional relation 


H=)'4H,,, (3.30) 





No. I, 1950 The effect of turbulence on a magnetic field 79 


where A= O(4z0vd), v and d being the typical velocity and diameter, respectively, 
of the turbulent eddies. In the Sun’s convective core for example we have 
o~7 x 10-*e.m.u., while Cowling has estimated (7) o~2 x 10%cm./sec. and 
d~10°cm., giving A~2 x 10. Now, H,, =qH,/A, H =qH,/A"?, hence the mean 
field resulting from a given E.M.F. in the convective core would be reduced by a 
factor of order 10! by turbulence. 

This effect must not be confused with the distortion of a field originating in 
some remote non-turbulent region, such as the distortion, in the hydrogen convec- 
tion zone, of a given field emanating from the radiative envelope. Here H can 
be very much greater than, while H,,, is of the same order of magnitude as, the 
field in the radiative envelope. For intense turbulence the conductivity in the 
turbulent layer would be negligible, so that H,, is simply the irrotational completion 
of the field at the boundary of the radiative envelope, and as such is easily calculated 
without further knowledge of the turbulence. 

4. The turbulent conductivity in the convective core and in the surface layers.— 
The problem is to find the order of magnitude of the ratio of the mean scalar 
field strength to the strength of the mean field, in these regions. It will be assumed 
that turbulence takes the form of a combination of two types of motion, pure 
rotatory motion, or eddying, in small cells, and rapid translatior. of small portions 
of the gas which on moving from point to point isolate themselves from their 
initial surroundings. This distinguishes turbulence from the type of motion 
in a system of quasi-steady convection currents in the form of long filaments, where 
the material at any point is not isolated from the material on either side of it along 
the filament, although the filament itself may be constantly changing its surround- 
ings. Even in this case, of course, there may be turbulence at the surface of the 
filaments. 


























Fic. 1.—The effect of an eddy in turning through one revolution. 


The effect of eddying on the field is to wind up the lines of force like clock 
springs, as illustrated diagrammatically in Fig. 1. 

The effect of turbulent translation is as indicated in Fig. 2. As the cell 
isolates itself from the material initially surrounding it, so also does the field it 
carries isolate itself from the surrounding field. As the cell moves through the 
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field it wraps lines of force round itself which break away in succession from the 
parent lines as the surrounding material closes in behind the cell. Fig. 2 shows 
the correspondence between the lines of force of the isolated field inside the cell 
and their parent lines in the general field; the limiting line of force about to be 
broken away is shown by the heavy line. The point B on this line where the two 
fields are attached is a null point of the field, as is also the point A within the cell. 

A finite solenoidal vector such as a magnetic vector can only vanish on closed 
lines, although in special cases these lines may coalesce into surfaces. In Fig. 2 
the null line involved is a line through B perpendicular to the plane of the figure. 
This is the line of attachment of the broken piece to the main field, and it joins up 
with a similar line through A running inside the cell, thereby completing a closed 
curve as required. 












































4 
Initial position of cell 
Fic. 2.—The effect of turbulent translation. 


It will be seen therefore that, on this model, no line of force is displaced by 
much more than the typical diameter of the turbulent cells, although pieces 
broken from the field may travel farther. The broken pieces do not of course 
contribute to the mean field. These considerations justify that step in the previous 
section where it is tacitly assumed that it is possible to select strips which are long 
compared with the dimensions of the turbulence, although short compared with 
the distance of variation of the mean field. It should be remarked here that the 
presence of isolated pieces of the field must modify the detailed analysis of the 
preceding section, particularly as regards equations (3.10) and (3.11). It is the 
writer’s opinion that the final result as given by (3.19) is not affected. The result 
may not be very accurate, however, where the turbulence is reduced to a system 
of convection currents as is apparently the case in the hydrogen convection zone, 
and possibly also in the convective core and in the chromosphere. Hoyle’s 
theory (5) of the chromosphere, for example, is based on a study of long convection 
columns driven by the accretion of interstellar material. However, even in long 
convection filaments, where material is constantly being interchanged between 
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filaments, it may still happen that no one portion of material actually travels very 
far along a filament before crossing to an adjacent filament, thereby becoming 
isolated from its surroundings. 

The intensity of distortion of the lines of force is limited by one of two factors. 
It has been shown that the mean scalar field strength is increased by the turbulence, 
as compared with the strength of the mean field. This scalar intensity cannot 
increase far beyond the point where the field begins to affect the dynamics of the 
turbulence. This occurs when the magnetic energy density H?/87 of the distorted 
field is comparable with the kinetic energy of the turbulence. Alfvén, for example, 
has suggested (6) that a scalar intensity of between 160 and 500 gauss may be 
present in the photosphere, owing to convection there. The second possible 
limiting factor depends on the normal conductivity. Bearing in mind the picture 
developed in Section 2, together with equation (2.17), it can be seen that, as the 
curvature consequent on the distortion increases, the lines of force, in attempting 
to straighten themselves, slip back faster through the gas distorting them. 
A point will be reached when the lines straighten themselves as fast as the gas 
eddies tend to distort them and no further distortion takes place. The intensity, 
when limited in this way, would be a very large multiple of the. initial mean 
intensity. (2.17) shows that the maximum mean curvature of the lines is given, 
in order of magnitude, by 


K=470W, (4.1) 


where « is the mean scalar curvature. Since any small bulge in a line of force is 
straightening as fast as the turbulence is distorting it, w must therefore be the 
velocity of the gas in that bulge relative to the velocity in adjacent bulges. Now 
the turbulent velocity v of the gas varies over a typical distance of the order of the 
diameter d of the eddy cells, hence 


w~vil/d, (4.2) 


where / is the typical diameter of the bulges in the lines of force. But J~1/k, 
hence on eliminating w between (4.1) and (4.2) we have 


I~d)X*”, (4.3) 


where A=4zo0vd. Hence in any distance L along a mean line of force the corre- 
sponding distorted line of force doubles back on itself A”*L/dtimes. It has been 
shown that the amplitude of the distortion of the lines as a whole is of order d. 
Consider then a region of diameter d as traversed by a certain tube of force. The 
contribution to the mean field from this tube is simply the strength of the tube. 
But in a distance d the tube doubles back on itself A”? times, hence the total contri- 
bution of the tube to the mean scalar strength of the distorted field is A’? x strength 
of tube. Hence 


H~™?H,,. (4.4) 


As was mentioned in the previous section A”*~105 in the convective core. 
The largest field yet accountable for theoretically is of order 1o~* gauss in the core, 
as given in a theory (2) by Cowling. Thus, neglecting factors such as turbulence 
tending to reduce it, the field could not be less than this. The limit which would 
be imposed by the conductivity is therefore of order 10 gauss. This, as is shown 
later, is less than the limit imposed by the available energy of the turbulence. 
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In the surface layers the factor A"? would not exceed 10° owing to the lower 
conductivity. Cowling’s field is of order 10~* gauss in strength at the surface; 
the distorted scalar strength of this field would therefore be limited to about 
10-* gauss by the conductivity, which again is less than the limit imposed by the 
energy of the turbulence. However, since it is generally thought that a mean 
field of order 30 gauss, of unknown origin, is present in the surface layers, it will 
be assumed that the scalar distorted intensity is limited in the surface layers and 
in the core by the energy of the turbulence. 

The turbulent conductivity is therefore best expressed as 


o7=0q (=). (4-5) 


where H is a function of the turbulence only. a, thus depends on the mean field, 
which greatly complicates the study of the electromagnetic equation. It should 
be noted, however, that the expression in (4.5) only holds where the strength of 
the mean field is less than the scalar intensity available from the turbulence. 
Larger mean fields than this tend to prevent turbulence in any case, allowing 
motion only parallel to the field. Sunspot fields in particular exceed the level 
of the turbulent intensity in the photosphere as assessed by Alfvén (6) and are 
therefore not appreciably affected by it. In spot-free areas of the photosphere 
where it is assumed that the mean field is of order 30 gauss and the turbulent 
intensity of the order of 300 gauss, (4.5) shows that the effective conductivity is 
reduced by about 100: 1, toa value of order 10-e.m.u. A reduction of this order 
of magnitude is of little importance; it would allow a drift velocity of gas across 
the field of the order of 100 times the amount accepted hitherto, but still smaller 
than the drifts observed. 

If, however, the mean field were found to be considerably less than the 
classical value of 30 gauss, say about I gauss, as recently measured by Thiessen, 
the conductivity would be reduced by about 105:1. Latitude drifts of up to 
10°cm./sec. would then be possible without appreciably affecting the field. 
The results of the author’s previous paper (1), however, show that a reduction of 
10’ or more in the effective conductivity would be necessary to permit an appreciable 
variation in angular velocity along the mean field. 

In the convective core, also, the turbulent intensity is probably much greater 
than the mean field. Following Cowling’s theory (7) of the convective core, the 
kinetic energy $pv* of the turbulence is of the order of 2 x 10’ergs/cm.°. This 
would produce a scalar intensity of the order of 70,000 gauss; thus if the mean 
field were taken as 10‘ gauss, the effective conductivity would be reduced by about 
50:1. Such a reduction has two main consequences; first, in theories seeking to 
explain the Sun’s field by E.M.F.s in the core, the E.M.F. to be searched for must 
be greater than previously supposed, as was shown at the end of Section 3, and 
secondly, the decay time of the field, if decaying from a field built up at some remote 
epoch, must be somewhat less than that derived from calculations based on the 
normal value of the conductivity. This calculation is made difficult by the 
non-uniformity of the conductivity throughout the Sun. Using Rayleigh’s 
principle, Cowling (2) calculated this decay time to be approximately 107° years, 
and the eigenstate distribution corresponding to an intensity of 25 gauss at the 
equator at the surface, to give 2000 gauss at the centre. With this value of the 
central field the conductivity in the core would be reduced by the turbulence, 





No. I, 1950 The effect of turbulence on a magnetic field 83 


from its normal value of about 7 x 10~*e.m.u. to approximately 6 x 10-7 e.m.u. 
These rough considerations show that the decay time depends, not on the core, 
but on the inner part of the radiative envelope surrounding it, and that the 
field in the core must decay relatively quickly to an irrotational completion of the 
field in the radiative envelope. The actual calculation of the decay time of the 
whole field will not be attempted in this paper. 

5. Conclusion.—(i) Turbulence reduces the effective conductivity in the core 
and in the outer layers of the Sun. Sunspot fields are not affected however. 

(ii) The decay time of the Sun’s general magnetic field is somewhat less than 
the figure of 10! years given by Cowling (2), although not by a factor of any great 
consequence, while the mean field in the core, if the general field is decaying from 
some initial state, is irrotational. 

(iii) No system of meridian-plane convection currents in the Sun can provide 
a general amplification of a field produced by a given E.M.F., while turbulence 
in fact reduces the field. 


I would like to thank Professor T. G. Cowling for his invaluable criticism of 


the paper. I am also indebted to him for suggesting the short proof given of the 
lemma in Section 2. 
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Glasgow : 
1949 December 5. 
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P. 426, line 7, for wn read w/n. 

P. 427, equation (22), for the first term on the right-hand side read 'f(x) s—}4. 
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Summary 


Equations due to McCrea are applied to distinguishing between (a) the 
amoun tof radial separation allowed by viscous drag among circumstellar gases 
under different radiation pressures, and (b) the more common interpretation 
of ‘‘ stratified atmosphere ”’ as a radial variation in stages of optical excitation 
due to screening and attenuation of radiation without separation of materials. 
A range in variation of radiation pressure for McCrea’s theory is obtained by 
adapting the equations of Gerasimovic. Asymmetric ejections of gases are 
classified where they can present spectra of high and low excitation together 
and thus simulate a stratified shell. Using Grotrian’s time for an absorbing 
envelope to adjust itself to change in radiation, as well as McCrea’s time for 
gases to adjust their velocities, some limits are found to the delay in spectral 
changes imposed by rate of uncovering of radiation in a binary system, and by 
the rate of penetration of a screen. Simple application is made to some 
features of Z Andr, ¢ Auri, 8 Lyra, p Cass and y Cass, as elucidated by 
Struve and Swings, Merrill, Kuiper, Kopal, Baldwin and others, preliminary 
to studying the more violent motions in W binaries and novae. 





1. Movements in circumstellar gases contributing to rates of spectral change.— 


Among the problems of non-steady stars, some of the most complex arise from 
the presentation of spectra implying high and low excitation simultaneously or 
in rapid succession. ‘The term “stratification’’, frequently used in discussing 
such problems, does not always leave certain whether it means (a) that differential 
gas motions have enforced a radial distribution in composition, or (5) that dilution 
or screening has provided radial differences in the stages of excitation without 
necessarily any chemical separation. But this distinction is itself incomplete if 
the simultaneous spectra arise not merely at different radial distances from a star 
but from different portions of the area presented; so the alternatives (a) and (b) 
involve the additional possibility (c) of asymmetry due to equatorial or local 
extrusions in a single star or to interference between the atmospheres of a binary 
pair. The degree of such symmetry can range from the prominence localized in all 
dimensions in Kopal’s suggestion for { Auri, to include the spiralling stream 
from Kuiper’s model for 8 Lyra, and then the separated rings around RY Pers 
or the optical structure in the equator of y Cass; even a shell of P Cygn type, 
when as small as that said to surround ZAndr, must exhibit some slighter 
asymmetry if its source is a component of a binary. 

Any rate of change of the spectrum may either denote development and decay 
of the conditions underlying (a), (6) or (c), or, if the object is binary, may be 
controlled mainly by the orbital motion giving atmospheric eclipses. Single 
instances of minima in the light-curve may be due to random drift of interstellar 
cloud, suggested as one of the alternatives in discussion of N Herc (1934) but an 
unlikely coincidence for frequent repetition; asymmetric ejecta allow more 
chances of repeated interposition of a screen from rotating stars and binaries. 
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Interposition of a short-lived prominence may never recur, whereas inequality 
in recurrence of abnormal opacity in 6 Lyra has been attributed to the gradual 
development of a ring from a jet with a maximum density in the second convolution. 
Rotation of the inner portion of a shell around Z Andr or « Scor, relative to its 
now detached nebulosity, offers tiuctuations even more unequal in repetition if 
the object is a binary or a single source of asymmetric ejecta. On the other hand, 
ascription of fluctuating opacity to local association and dissociation in situ is 
attempted in ter Haar’s model of N Herc (1934), also for p Cass in 1947, and in 
explanations of changing shell transparency for Z Andr. 

A complete discrimination between these possibilities must ultimately account 
for duration and rate of extinction and repetition of any spectral change, in terms 
of geometrical factors of rotation and orbit and physical factors of local change 
in composition, density and excitation. But that task awaits two types of decision : 
(i) Can gases acquire such differential impulses, e. g. by selective radiation pressure, 
as to maintain separation as shells without degenerating into common motion by 
viscous drag? (ii) What time elapses (A) during establishment or decay of 
such separation, (B) during adjustment to changes in radiation, (C) during the 
emergence into line of sight of one or other portions of a rotating and asymmetri- 
cally erupting atmosphere or one of the components of an eclipsing binary ? 

Preliminary exploration into some aspects of this complex problem seems 
possible by applying McCrea’s (x) criteria of gas separation, and Grotrian’s (2) 
“time of relaxation of a system exposed to changing radiation’’, to some of the 
published evidence on complex variable stars; from these we assemble a set of 
examples. It is hoped that atentative decision on some of the above ambiguities 
in these well-known stars may enable progress to be made towards the even more 
complicated simultaneous spectra of novae, e.g. where Stratton and Gaposchkin 
postulate jet models of N Herc (1934), and in the W binaries of Beals, Wilson, 
and Gaposchkin, where the greatest elaboration is to be expected in the penetra- 
tion of one star’s materials by the ejecta from its companion. 

2. Data.—The groups of facts, and inferences by their observers, to which 
we shall relate our calculations are as follows. 

(i) ZAndr. A sequence of researches, especially by Merrill (3) and by 
Swings and Struve (4), on the shells ejected in 1939-40 and 1946, discusses the 
separated components of lines at different speeds. In their interpretation, 
Swings and Struve postulate a small P Cygn shell, of the order of 10" cm. com- 
pared with Beals’ mean 200 shells in more typical P Cygn stars: the sub-dwarf 
companion to R Agar also has been associated with a small shell radius. On 
the other hand, these writers call for a far wider extent of the subsequent nebulosity 
here as in « Scor, 10! cm. or even 10!® where the forbidden lines [Fe II] arise. 
We shall compare the Grotrian calculation for the shell with Struve’s own use 
of the same method for recovery of the nebulosity when the shielding shell was 
penetrated in 1940 August. 

(ii) Kuiper’s (5) model of 8 Lyratraces a spiral stream from the Bg component 
flowing across the F component to a diminishing density discussed by Struve, 
and the location of the point at which the asymmetric ejection occurs is significant. 
The transient opacity preceding eclipse of the B component (r=2-80) of ¢ Auri 
in 1939 is ascribed by Kopal (6) to a giant prominence at 1-7 x 10'! cm. above 
the K component’s photosphere (r=2000). McLaughlin (7) finds recent 
evidence for some permanent asymmetry of the K component. 
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(iii) Baldwin’s (8) model of y Cass, rotating at variable speed over a pulsation 
cycle, enables our treatment to compare the meanings of stratification when there 
is a difference of 70 km. per sec. between H and He, in contrast to the behaviour 
of other ringed systems: many other Be stars show cyclical swing between 
double components of a single line, commonly H, more rarely Fe*+ and Ca*. 

(iv) Greenstein’s (g) treatment of the brief changes in radial velocity and 
opacity of p Cass in 1947 November, as explainable in terms of the collapse of a 
shell upon lower strata, when a burst of radiation ionized the distant hydrogen, 
making it nearly transparent and removing its radiative support, offers ground 
for exhibiting McCrea’s and Grotrian’s theories of adjustment, in a case resem- 
bling more closely Grotrian’s original model. 

3. Gas separations limited by viscous drag.—A well-known attempt to refer 
spectrophotometric changes to motions in gaseous ejecta is ter Haar’s (10) model 
of molecular clouds in N Herc (1934), based on suggestions by Chandrasekhar 
and Stratton. The single depression in p Cass has attracted somewhat similar 
ideas, which merge into the older “‘veil”’ theory of long-period variables and 
may also be relevant to the irregular dimming of R Cor B stars. ‘The meeting, 
mixing and separating of gases postulated in all such theories, and more readily 
examined in the simpler stars of the present paper, would occur when stellar 
radiation fluctuates or even where steady radiation becomes differentially absorbed 
by mixed gases; the motions would certainly exemplify our first meaning 
(Section 1) of stratification, if initial impulses due to radiation pressure can be 
maintained until shells develop. Similar considerations apply along radii when 
the asymmetric streams we have mentioned can only provide segments of a shell. 

Optical changes will accompany such altering composition of a mixture, of 
partial concentrations m,, mg, .... per unit volume, when dn,/dR, dn./dR .... 
reverse in sign at various values of R, permitting maxima to develop in the partial 
pressures of constituents with differing opacities. To decide the possibility and 
the time for establishment of such segregations, we proceed to apply McCrea’s 
criteria; these state where the interchange of momentum at atomic or ionic 
collision can permit initially differing accelerations to maintain a difference in 
radial velocity. The equations needed from McCrea’s paper are as follows. 

If g,, Zo are the effective values of radial acceleration from a star for two 
substances, or for one material in two states of excitation and therefore subject 
to different radiation pressure, and p is the density often determined mainly by 
one constituent, then McCrea’s velocity U of separation permitted by viscous 
drag is given by 

dU/dt =g, —g2.—pyU*. 
We shall assess below (Section 5) the range of (g,—g.) including both gravity 
and radiation pressure. McCrea distinguishes two approximations, U<v when 
€=1 and y=a, and USv when «=2 and y=, where v is thermal velocity. 
The quantities « (temperature-dependent because of the intrusion of v) and 8 
(independent of temperature because v is relatively smaller) can be evaluated 
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The o’s and m’s are collision diameters and masses of the two kinds of atom. 
McCrea justifies in considerable discussion the approximation of gas-kinetic 
diameters for stellar ions. In the first case a limiting velocity of separation U, 
is attained, and in the second case Uj. A defined “time of relaxation” T, 79, 
governs rise or decay in transition, for a given fraction or multiple of each limiting 
value :— 


Up =81 —82/p%, Ty =1/pa, 


Ma 81 —8> _— SSE EES 

eb A/ pp’ 7 V eB(gi—82) 

For the purpose of our applications we compute « and £ for a set of gases relevant 
to our stars, Fig. 1 indicating that the possibility of stratification by actual 
separation between H and any other material increases tenfold between He 
and Fe, and that the greatest sensitivity to temperature in the temperature- 
dependent approximation occurs in the lightest pairs. We find also that for 
non-hydrogenic mixtures the dominant variable is the product m,m,, 8 altering 
by only about 10 per cent for the He mixtures while « falls to half of its (hydrogen) 
value, which means a velocity of separation rising by 100 per cent. 
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It is apparent that besides the variable p explored by McCrea in his applications, 
e.g. to solar prominences, the difference (g,—g.,) will become decisive in the 
larger mass-motions of erupting stars and of Struve’s non-static shells. The 
possible range of (g, —g.), treated below in Section 5, will be satisfied if we take 
extremes from 1/100 to 10 times the solar value of g(grav.). This will allow 
inclusion of velocities of separation as function of concentration for pairs of gases in 
both giant and dwarf atmospheric extrusions at high and low temperature, with 
extremes of 60,000 deg. and 6000 deg., with instances when one gas is dominated 
by gravity or rotation and the other gas subject to either a balancing or a reducing 
or a dominating radiation pressure. Accordingly, from Fig. 1 we plot Fig. 2, 
covering a sufficient range to read off the separation velocities for pairs of gases 
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with extremely differing absorptive properties. Both approximations Uy and Uj 
are plotted, as in critical cases the permission to separate may differ at different 
thermal velocities, which are themselves not homogeneous around an erupting 
star. Arrows on the plot indicate what displacement from the graph of H/He 
will cover other likely mixtures with hydrogen, and we recollect that with non- 
hydrogenic mixtures U is somewhat increased. 

Fig. 2 shows sufficiently how velocity of separation can rise very high at 
nebular concentrations (v<10> per cm.'), but at high chromospheric and 
prominence concentrations (v>10*) can only become appreciable at the largest 
reversed values of g conferred by selective radiation pressure and leading to 
motions actually observed in eruption. The separations are considerably 
reduced by temperature. At v>10! separation from hydrogen is inappreciable 
for He even under the highest radiation pressure, and is only incipient for 
heavier atoms. 


Log,, U(H/He) 
in km. per sec. 




















U, Logi 
Fic. 2. 


4. Time of adjustment of shell or segment to fluctuating radiation —In 
spectrophotometric changes following gas motions, adjustment of the system 
depends not only upon momentum exchanged at collision as in McCrea’s theory, 
but also upon the fact the: some of the relevant values of g are due to radiation 
pressure: for an extens:ve circumstellar atmosphere will require finite time 
in reacting to even an instantaneous change of radiation. Merrill’s lag of 1/5 
period at some of the fluctuations of very large shells may well be such a finite 
travel of the locus of self-adjustment. It becomes relevant to apply Grotrian’s 
problem of a hydrogen sphere at whose centre a Planck spectrum at enhanced 
temperature is established within an interval brief compared with that needed 
for the gas to come into equilibrium. It is important to note that Grotrian’s 
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approach has not the immediate generality of McCrea’s; he considers for one 
band of frequencies in one element, the time ‘‘7” for ionizing absorptions and. 
recombinations to raise to within a given large fraction of its maximum the fluor- 
escence stimulated by a Lyman continuum as studied by Zanstra. Only to rough 
order of magnitude does this give information about the over-complex problem of 
unspecified gas mixtures exposed to fluctuations in other frequencies; but 
hydrogen is a striking feature of most stellar eruptions, and the equations we 
shall quote indicate where other ionization potentials and other temperatures 
both rising together enforce upon other materials a comparable behaviour. 
We point out some differences in Section 7. Grotrian applied his calculation 
to the recovery of N Herc (1934) after its minimum of 1935, and Struve and Swings 
have also applied it to the recovery of Z Andr after its P Cygn shell had ceased to 
impede radiation from the star in 1940. Adaptation to any segment of the sphere 
is legitimate except at the bounding edge of asymmetric distribution such as a 
jet, since R is the only dimension involved. But the serious limitation to which 
we refer in considering generalization in Section 7 is that of total absorption 
as end-point, demanding a sufficient depth of gas; for thinner atmospheres, 
Grotrian’s 7 is an upper limit only. 

If « is a constant controlling the rate of recombination and N the total number 
of ionizing quanta from the star absorbed per sec. in the surrounding nebulosity, 


r=1/2VaN, «=a'/V, 
where V is the volume contained by the gas shell and 
a’ =2/n'?; 

n’ is the number of ions, also of free electrons in this approximation, per unit 
volume, and z is the number of recombinations per unit volume to the mth state 
in Cillié’s (11) formula, 

Va75 10 3/2 

= Gate aie (*) M(m, T,)n"*. 
In Zanstra’s (12) discussion of fluorescence in circumstellar gases, 
87? R2x3 


re x2 
N= — za T* | —— dx 
ch J, e@-1 


x=hv/kT and in this case v is the frequency at the Lyman limit. We utilize 
Cillié’s tabulation of the function M(m,T,) for m=2, ..., 10, and Zanstra’s 
tabulation of the integral, and we compute 7 for 15,000 deg. and 20,000 deg.. 


TABLE I 
T in seconds 
r=10! 108 
2°4 X 108 | 7°5 x 108 
6°3 x 10? 2°0 X 108 
2°4 X 10? 7°5 X 10° 
6-3 x 10! 2°0 X 10° 


7°5 X10? 
2°0 X 10? 








go Martin Johnson, Adjustments within shells and asymmetric Vol. 110 


for the cases of R (of star) 10!°, ro", 101", and 7 (of shell or nebulosity) 1o!", 101%, 
to, 1015, The heavy line marks the conditions in R and r, to the right of which 
the time of adjustment reaches the order of a day and rises to the order of 
three years, a range covering the stellar distances we here consider. 

5. Differences in acceleration contributing to incipient gas separation.—We now 
consider grounds on which the range in (g, — g2) was selected for applying McCrea’s 
theory in Section 3, and the conditions bringing our particular applications 
within that range. The effective radial acceleration of a particle outside a star 
at distance R is 


g(eff.) =[GM/R? —w?R] —g(rad.). 


The following will be relevant sources of variation in this quantity :—{i) Since 
Struve’s shells have radii several times the stellar radius, the latter itself greatly 
exceeding solar in some cases, the first term on the right-hand side can alter 
enormously between photosphere and outer envelope but to the same extent 
for all atoms; in an extreme case, ZAndr can exhibit spectra at values of R 
perhaps 10% times solar. (ii) In rotating stars the second term rises at the 
equator, until after expansion to large R, w is diminished by viscous drag; 
Baldwin’s discussion of yCass, where line widths fluctuate over a pulsation 
cycle between 260 and 480 km. per sec., compares the maximum velocities 
permissible under rotational stability, 300 in Jeans’ liquid star, with 500 needed 
to make the second term exceed the first at the M and p of a Be star. Elvey’s 
classical measurements trace a single line from 50 to 250 km. per sec. in a series 
of stars. (iii) The third term is the only one to vary from element to element 
at a fixed R, both by differential absorption and by monochromatic screening 
and backward re-radiation, and by Milne’s mechanism of escape from a screened 
frequency through Doppler displacement. (iv) During temperature changes, 
e.g. in yCass, ZAndr and the flash of p Cass, the third term will also alter in 
integrated amount as well as among different elements. (v) During binary 
eclipse, e.g. of B Lyra and of {Auri, the third term can alter suddenly on 
emergence of one star from the shadow of a companion which is not radiating 
symmetrically over its whole surface. 

The analysis of Gerasimovic (13), deriving the absolute value of g(rad.) in 
the particular case of absorption within the Lyman continuum, and showing it 
2 or 20 times the gravitational term in a Be star even without rotation, can be 
adapted as follows for tracing the variation among other elements :— 


g(rad.)= G=oP" { «(v)I(v) dv, 


mc 


sd 


(1 —x) and f are the fractions, neutral and in a given state, appropriate to ionizing 
absorptions in the continuum beyond », r is residual intensity; corresponding 
quantities take their place at higher continua for second ionization, etc., in other 
elements than hydrogen. The rise of g(rad.) from 10 to 10 cm. per sec.? in H 
between 7000 deg. and 15,000 deg. is checked, soon after exceeding gravity, by 
the fall in the neutral fraction and, for this case of a transparent ion, Gerasimovic 
finds a dependence upon electron pressure making the main expansions begin at a 
low level almost photospheric. The high velocities observed in P Cygn shells 
up to 300 km. per sec., 1000 in W stars, 2000 in novae, involve not only the drag 
by the hydrogen but the fact that no other materials are confined to neutral 
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absorption only. To obtain some insight into likely differential effects contribut- 
ing to McCrea’s (g,—g.), we plot in Fig. 3 the most variable of Gerasimovic’s 
factors, log J(v)/m, where m is atomic mass as multiple of that of hydrogen. ‘The 
other important variations will be in “dilution” and “screening” effects upon r, 
and the time spent in a state capable of the particular absorption, represented 
here by (t—x). Gerasimovic derives his rise and fall of radiation pressure on H 
by combining the temperature coefficients of (I — x) and J(v), but for the present 
purpose of comparing elements these are better exhibited by graphing the actual 
intensity at the series limit, and indicating separately where the control by (1 — x) 
makes one or another portion of the several graphs significant. Accordingly 
we bias our graphs by thickening the line over the range where fraction in the 
relevant neutral or singly ionized state has not fallen below 10 per cent: the 
required information we derive from an auxiliary graph of logmt/n (14) 
according to ionization potential and temperature. The actual range over 
which Gerasimovic found H liable to radiative extrusion from a star is dotted 
on our graph, and extends to n+/n->10%, showing that for comparable J(v)/m we 
are underestimating rather than overestimating the liability of other elements to 
motion under radiation pressure. 


Log 10) 
m 

















10 6x 108 
TEMPERATURE 
Fic. 3. 

Fig. 3 indicates that in the temperature range where Be stars can extrude 
their hydrogen, metals in their: exceedingly brief neutral states would be com- 
parably accelerated, but in their much more prevalent ionized states would suffer 
accelerations 1/10 or 1/100 that of H, rivalling the latter only when it becomes 
transparent through ionization. Neutral gases of ionization potential comparable 
with H fall to below 1/10 of the latter’s acceleration, but experience line absorptions 
pending second ionization. At Be temperatures the graph indicates factors up to 
108 or 10 between the J(v)/m of different atoms, each having considerable abundance 
in the same excitation stage; so that although at Wolf-Rayet temperatures all 
elements may be moving outward and at solar temperature only rare local 
prominences, B stars will be particularly rich in large values of the critical (g, —g.) 
between pairs of gases. Of any pair, both experience the same g(grav.), large 
near the photosphere and small in a distant shell or modified at a rotating equator ; 
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but one of the pair may show an effective g reversed compared to the other, and 
a negative multiple of g(grav.) when its ionization continuum is favourably 
situated in the Planck spectrum of that temperature. Such cases will be covered 
by our selection of (g, —g,) between a pair of elements up to 10 times and down 
to 1/100 of the solar g(grav.) in Fig. 2. 

6. Applicsations.—In proceeding to find where a few particular stars may 
exemplify the features obtained from our graphs and table of McCrea’s and 
Grotrian’s arguments, it will suffice at this stage to work to orders of magnitude 
only. 

(1) Shell and nebulosity of ZAndr.—({a) Viscous drag conferring common 
shell velocity: if the estimate which we quoted from Swings and Struve for the 
P Cygn shell of Z Andr is correct, Fig. 2 indicates a high degree of drag. For 
strong lines in so small a shell, concentrations will probably exceed 10" per cm.®, 
and we see that even with (g,—g,) at its maximum of I0 g, concentrations 
exceeding 10!” per cm.* will prevent separation between He and H or H and O 
from reaching 10 km. per sec. and only allow H/ Fe slightly to exceed that separation 
under the most favourable conditions of temperature. Since PCygn shells 
are not themselves an indication of localized jet extrusion, and we classed Z Andr 
among the lesser degrees of asymmetry, it would seem that the differential line 
velocities in the 1939 and 1946 extrusions must represent stratirication of different 
optical stages, not of different materials separated, the alternative (6) in our first 
paragraph. ‘These stages must interchange absorption and emission predomin- 
ance within quite short distances, the dilution of radiation being due to steep 
local absorption gradient rather than here to any inverse-square attenuation. 
It is arguable that the phase lag between line velocities in those portions of the 
Z Andr system which appear to pulsate must also denote strata of excitation only, 
not of material, perhaps complicated by the lag in propagation of response which 
we discuss below: whether similar argument can include the multiphase pulsation 
of Mira stars will depend on finding whether all their visible gases can also be 
held to exceed 10! per cm.* concentrations. 

(6) Gas separation in nebulosity: the above conclusions here become 
reversed. The [Fe II] region at 10! cm. is perhaps large enough for asymmetry 
to have become smoothed out, but concentrations may well have fallen to a 
nebular 10* per cm.* and cannot exceed the very high nebular 10° per cm.? of the 
“stellar planetary’? IC 4997. At that maximum density Fig. 2 indicates that 
separation between any gases absorbing as selectively as in Fig. 3 can readily 
exceed 100 km. per sec., and indeed there is nothing in the viscosity to prevent 
separations up to 10* km. per sec. being reached, exceeding Wolf-Rayet and 
nova ejection velocity, for He/H and Fe/H if (g,—g.) becomes sufficiently great. 
The net effect will only remain small when the radiation pressure is a mere 
correction to a predominant gravitation. 

(c) Time of adjustment to radiation: Struve and Swings have already 
applied Grotrian’s equations to estimating three or four weeks for the hypothetical 
“old nebulosity”’ to react to light restoration when the P Cygn shell became 
transparent in 1940 August. We add to this, from Table I, the very different 
situation in the shell itself; at the estimate quoted, where shell radius is five 
times the small stellar, the time of relaxation of a completely absorbing hydrogen 
shell following alteration of the Lyman radiation will be certainly much less 
than a day, and for many values of the quantities included in Grotrian’s equations, 
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only afew minutes. This ensures the premiss omitted in the argument of Struve 
and Swings concerning the outer nebulosity, that time for the alteration of the 
screen must be negligible compared with time for the more distant target to 
adjust itself. 

(d) Time for adjustment of gas velocities: in comparison, in the [Fe II] 
region where alone gases can separate, stabilization of the separations in McCrea’s 
theory can be shown from our Section 3 to occupy the following times, at a 
concentration of 105 per cm.*, at 10,000 deg. The times are maxima, taking 
the smallest (g, —g,), of the order of 10*, capable of providing large separations. 


T,(H/He) =5 hours, T,(H/ Fe) =84 hours, 
To(H/ He) = 5% hours, T((H/ Fe) =19 hours. 


(II) Adjustments among the ejecta from binaries.—{a) Viscous drag:—({i) In 
RY Pers, Hiltner (15) has distinguished rotational velocity range of 200 km. 
per sec. in He compared with 110 km. in H, both rotating round the F5/B4 pair. 
RZ Scut and several other binaries also present slowly rotating H shells or rings 
around faster spinning He with velocity ranges up to 250 km. If this is a 
separation by McCrea’s mechanism, Fig. 2 shows that concentrations must be 
far below 10! per cm.® where any large (g, —g.) becomes capable of enforcing a 
noticeable U. But the process leading to stratification is here susceptible of 
a different solution from that expressed in Section 5; it may not have developed 
from the differential radiative accelerations within an initially mixed gas, if 
Struve’s suggestion is correct that such ringed binaries with sharp lines are 
exhibiting the debris of their original formation, not contemporary ejection. 
(ii) The close binary 8 Lyra and the more open ¢ Auri are not instances of spherical 
shell stratification at all, and indeed Fig. 3 suggests that radiative extrusion is 
unlikely from the cooler components of either except at particular points where 
gravitational attractions from two components can compete. But consideration 
of those points allows McCrea’s theory to become relevant in the novel question 
of lateral drag between adjacent streams side by side, as follows. Kuiper 
distinguishes an “A” type of ejection from contact binaries, emerging at a point 
between the two stars and wrapping round the hotter with a large apsidal 
rotation which causes intersection of streams. At distances not large compared 
with the stellar dimensions, the ejecta will therefore possess all the requirements 
for exchange of momentum to enforce common velocity. ‘The binary models of 
AX Pers and ZAndr are assignable to this “‘A” form, but 8 Lyra belongs to 
Kuiper’s ‘‘ B”’ form, where ejection occurs instead at the distant edge of the smaller 
star. Either model allows adjacent regions of excess or defect density, and the 
meeting or divergence of material from the two stars; in general the mixing will 
occur at densities above the minima needed in Fig. 2 for functioning of McCrea’s 
drag, but configurations of a binary system are possible in which a “ pocket”’ of 
abnormally low “‘g”’ between the stars allows radiation pressure to become 
effective at unusually low temperature: separation of gases would then be locally 
possible. Kuiper’s “ B”’ mechanism produces a ring only, whose optical develop- 
ment in the second spiral convolution is a measure of the collisions exchanging 
momentum; the slightness of contribution from the F component in f Lyra to 
its B5 ring denotes the shallowness of the former’s atmosphere, but such depth 
as the Bg stream penetrates into the F atmosphere will give a very complete drag 
of the one stream into the other, owing to atomic concentration at such small 
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distance from the target star lying well within the limiting region in Fig. 2. 
Struve’s (16) observation of a retardation from 300 to 200 km. per sec. on return 
of the ejecta, and the observed turbulent broadening of lines, are further indications 
that the gas filaments could not maintain separation for long even if they had 
initially differing accelerations. This instance of the carrying of adjacent gas 
at lateral contact, instead of solely when strata overtake one another, is applicable 
to the spread of jets into cones in some models of novae. (iii) In Kopal’s model 
of an obscuring prominence above the K component of ¢ Auri, gas so shallow and 
yet opaque enough for the observed dimming must be considered to have a 
concentration well within the “ drag” region of Fig. 2, and indeed Kopal estimated 
on other grounds 10" per cm.*. This drag would account for the obscuration 
being general, not monochromatic, even the extrusion of a single substance 
carrying the remaining gases with it—an extension of McCrea’s own application 
of his theory to the Sun. 

(6) Time lag in adjustment: although McCrea’s time of adjustment would 
be very brief in 8 Lyra, it is important to know whether Grotrian’s time for 
absorption of radiation is small or great compared with the time taken for a new 
target area to be exposed by orbital travel and rotation. Using Kuiper’s, Kopal’s 
and Struve’s estimates for 8 Lyra, R(star) =3-5 x 10!" cm., r(interposed ejecta) 
=I-I x 10!8 cm. as averaged over the first convolution, so that our Table I gives 
time for absorption of Lyman radiation less than 10° sec. Since exposure of the 
ejected material from F5 to radiation from Bg was estimated by those writers to 
occupy the fourth to fifth day after ejection, the radiative adjustment occurs in 
a time brief compared to even the relevant fraction of the short orbital period of 
12-9 days. On the other hand, in ¢ Auri, our Grotrian calculation for the linear 
dimensions R(K)=200Ro, R(B)=2-8Ro provides times of adjustment ranging 
from three hours to two days compared with the period of 973 days, according to 
assumptions chosen as to orbital geometry and asymmetry of ejection. Since 
the observed veiling of B occurs in 24 hours, when the new absorption lines in 
the B spectrum have already appeared several days previous to eclipse, Struve 
and Swings’ assumption of instantaneous exposure in using Grotrian’s equations 
for Z Andr is no longer valid in these cases : in Z Andr the fastest completion 
was the unveiling of the source of light and the slowest the reaction of the distant 
absorbing target, but in § Lyra the reaction is more rapid than the unveiling, 
and in ¢ Auri the processes are comparable in speed. McCrea’s times are again 
briefer, for  Auri Ty being less than one second and 7; not greatly exceeding one 
minute. 

(III) p Casstopeiae.—In Greenstein’s treatment, which we quoted, the sug- 
gestion is that the shell collapsed in 1947 November with an inward velocity reach- 
ing 31 km. per sec. after about a day. At the estimate R (star) = 1600 =10%cm. 
with a much larger shell, say, r=10'* cm., Table I gives times for adjustment 
2-4 x 108 sec. at 15,000 deg. and 1-5 x 10 at 10,000 deg., Grotrian’s model being 
almost exactly the process of hydrogen ionization postulated by Greenstein. Not 
only this response to radiation but also the adjustment of any differential gas 
velocity is rapid compared with the time of general collapse, since if we insert a 
minimum concentration 10° per cm. appropriate to this large shell, in our 
treatment of McCrea’s theory 


T,(H/He) = 1g sec., T,(H/Fe)=5 min., 
T((H/He)=10} min., 7 (H/Fe)=36 min. 
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At the density implied by the dimensions, real stratification in the sense of 
separation of materials is not prohibited if there were sufficient radiation pressure 
to cause it; but in so far as the spectrum suggests a K type this would scarcely 
be sufficient. In the “normal”’ star differences in velocity up to 7 km. per sec. 
were observed. When the “collapse” occurred Greenstein’s change from 
—60 to —12 km. in a day was drawn from 70 lines, so that if his explanation of 
loss of hydrogen support is correct, either the H was dragging other elements in 
the lower and denser portion of its path of fall, or the radiation had failed in its 
support of other absorbers too and all were experiencing together the common 
“‘g”? which Greenstein evaluates. Explanation of how all the materials had 
previously been extruded from the shell is not so forthcoming: radiation pressure 
is not the only agency, as has been proved in the Sun, but it may be relevant that 
before 1946 the star showed F character, and even when seemingly K or M its 
excited levels in Fe, Mg, Si, N, O exceed those expected at the contemporary 
temperature. So the envelope may be a relic, on small scale like a nebulosity 
surviving from an eruptive past. The evidence of Keenan (17) and of 
Thackeray (18) is against an erupting companion of the Z Andr or R Agar type. 

(LV) yCassiopeiae.—In Baldwin’s account of this as a single star, he lists 
stratification from upper to lower atmosphere in a large number of successive 
emission and absorption lines, but suggests that this need not imply separation 
of materials if the ‘‘ shell” is merely the optical effect of atoms streaming through 
a stationary absorbing region, visible only when each shows itself suitably excited. 
It is useful to check our treatment by asking whether it gives reason for such 
restriction to “excitation strata’’ instead of material separation. ‘The case is 
of special importance in view of the likelihood that here the equatorial ring 
extrusion of Be type is accompanied by pulsation, 7 and R varying with a period 
long compared to that of rotation; this confers periodic velocity change along 
radii but varying with latitude, so that if radiation pressure passes a peak at phase of 
high temperature, differential motions are bound to be initiated. Furthermore, 
the range quoted by Baldwin, 8000 deg. to 13,000 deg., does cover the critical rise 
to radiative extrusion in the Be stars of Gerasimovic. Baldwin’s stellar radius 
pulsates from 10 to 18 solar radii and his shell radius from 20 to 100; relative 
therefore to a 10! cm. nebular radius of Z Andr, where concentrations may have 
fallen to 10* or 10% per cm.°, the shell of y Cass or the ring, at one-hundredth of 
that radius implies concentrations of 10! or 10! per cm.%. At such densities, 
Baldwin’s separations up to 100 km. per sec. between H and He cannot be a 
separation of materials, Fig. 2 indicating that this would only occur below a 
maximum of 10* per cm.* at very high radiation pressure. Hence in any use of 
the term “stratification”, stages of excitation only must be involved as Baldwin 
suggested, unless future analysis of the periodic velocity in each line were able 
to prove that simultaneous radiation from photosphere and ring was capable of 
simulating a radial distribution. Ifthe former is true, the explanation we classified 
as (b) instead of (c) in Section 1, the lack of correlation between Baldwin’s radii 
and ionization potentials shows that no general “ dilution” of the inverse square 
type familiar in nebulae is to be expected; the local gradients in line absorption 
and population of levels for emission must be very steep as in the inner shell 
of Z Andr. 

7. Conclusion.—We have distinguished, with the aid of McCrea’s equations, 
some of the rare cases where gases can be driven apart under differential radiation 
pressure, from the commoner cases where stratification can only imply a radial 
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spread of the loci of excitation stages. With the aid of Grotrian’s equations 
added to those of McCrea, the time of adjustment of gas velocity or of response 
to changing radiation has been related to time of transient exposure in some eclipsing 
binaries and stellar extrusions. We have treated only a few aspects of five stars; 
in extension to other binaries with interpenetrating atmospheres, e. g. W binaries, 
and to recurrent novae and peculiar variables, the following are among the points 
needed to be taken into account. Whereas McCrea’s theory is subject to the 
limitations which he himself stated, and Grotrian’s treatment of hydrogen spheres 
is essentially adaptable to other elements and to sectors of a sphere, the whole 
simplification of adjustment to radiation will lose validity for light emitted 
beneath an insufficient depth of gas, since the equations depend upon total 
absorption. ‘This makes our estimate of Grotrian’s “‘r” to be upper limits only. 
It also introduces the important difference between the times of adjustment at 
mere eruption of a hitherto quiescent photosphere into empty space and 
eruption into a nebulosity surviving from a previous outburst as in Z Andr in 
1946. The spectral sequence of “repeater” novae, R Agar, T Cor B, SS Cygn, 


etc., may be dominated by terms in Grotrian’s “‘r’’ which do not arise when the 


outburst is the first which a steady star has ever perpetrated, and it may become 
possible to infer quiet or unquiet previous history from analysis of present rate 
of development. ‘The stratification of W binaries is a related problem, since 
here the ejecta incessantly impinge upon other ejecta, and will demand the greatest 
complexity in employment of the method whose simplification has only been 
possible in the few cases quoted. It is also worthy of note that in atmospheres 
where simultaneous application, e. g. to He and Ca*+ were attempted, in stars with 
such very discordant energy levels, the time of adjustment in fluorescent response 
in Grotrian’s theory will be so different as itself to constitute a source of separation 
of the stages in excitation—an optical stratification altering with time. When to 
a single change in radiation one fluorescent response takes days and another 
takes weeks to establish, a fresh set of interpretations will be necessary to some 
of the fluctuations in spectrophotometry. 


The University, 
Birmingham : 
1949 December 29. 
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SPECTROPHOTOMETRIC MEASUREMENTS OF EARLY-TYPE 
STARS. 1. METHODS OF OBSERVATION AND 
RESULTS FOR Oe5 STARS* 

E. A. Baker 


(Communicated by the Astronomer Royal for Scotland) 


Summary 


This paper gives an account of a technique of measurement and reduction 
of line intensities which has been developed at Edinburgh. The photo- 
graphed spectra are measured with a non-recording microphotometer at over 
a thousand fixed wave-lengths between A 4100 and A 6800, and the measures 
are combined for five or more negatives of each star, thereby reducing grain 
and other accidental errors without smoothing out real detail. The results 
so obtained for five of the Oes stars, selected as having narrow Balmer lines, 
have been further combined; and a mean continuum derived from these 
is used as a basis for the reduction of the results for individual stars. It is 
found possible by this procedure to free the results from the effects of 
continuum curvature due to the sensitivity curve of the emulsion. 

Results are given for seventeen lines in each of seven stars of Draper 
type Oes5. They are in good agreement with those obtained by E. G. 
Williams at Mount Wilson and are systematically lower than those found 
by Petrie at Victoria. In addition the broad “ interstellar’’ absorption 
centred near A 4430 has been studied and found to cover a wider wave-length 
range than previously supposed. For this limited sample of seven stars it 
is noticed that the central absorption of A 4430 shows a remarkably high 
correlation with colour or gradient excess, whereas the correlation of any of 
these with the intensity of the interstellar D lines is much weaker. Discussion 
of these results is postponed until further material is reduced and a larger 
sample of stars thereby rendered available. 


* The full text of this paper appears in Publications of the Royal Observatory, Edinburgh. 1. 
No. 2, 1949. 
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THE VARIABLES IN THE SELECTED AREAS AT 75° AND 
60° NORTH DECLINATION. 2. WORKING METHODS 
AND RESULTS FOR AREAS 5-7* 


E. A. Baker and R. W. Wrigley 
(Communicated by the Astronomer Royal for Scotland) 


Summary 


In continuation of a previous publication, this paper describes the methods 
used in a search for variables in these Areas, and gives results for Areas 5 to 7. 
Since results for Areas 2 to 4 have already been publishedt this completes 
the search of the Areas centred at declination +75°. The published results 
include only such stars as lie within three degrees of the centres of the Areas, 
become brighter than magnitude 14 at maximum, and show changes of at 
least 0-5 magnitude on more than one occasion. 

From plates taken with a ten-inch triplet eight pairs of each Area are 
“* blinked ” and a list of suspected variables drawn up, each with its own 
comparison stars. The photographic magnitudes of the comparison stars, 
accurate to about 0-1 magnitude, are obtained from various sources, that 
preferred being the Bergedorfer Spektral-Durchmusterung, extended in range 
by a short and long exposure method and in area by examining the field 
corrections on de-centred plates. If after a preliminary test the suspected 
variable promises to fulfil the above criteria, it and its comparison stars are 
compared on all plates by the step method and the results combined to give 
in the first place more self-consistent comparison star magnitudes and 
finally the magnitudes of the variable for the date of each plate. 

Maps with photographic magnitudes of comparison stars, positions correct 
to about one second of arc, and from twenty to thirty observations are given 
for each variable. In Areas 5 to 7 these particulars are given for two known 
and seven new variables. 





* The full text of this paper appears in Publications of the Royal Observatory, Edinburgh, 1, 
No. 3, 1949. 
t M. N., 975 541, 1937; 98, 65, 1937. 
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